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A Method for Integrating a Set of Ordinary 
Differential Equations Subject to a Type of 
Numerical Indeterminacy 


1. Introduction. This paper proposes the numerical integration of a system of 
first order ordinary differential equations by means of a convergent sequence of 
equation systems. The schema suggested may prove useful whenever: (a) the 
equations of the set are coupled in an appropriate fashion; (b) ordinary point by 
point procedures fail due to subtraction which makes the analytic formula for 
the derivative of one of the variables numerically indeterminate. 

The technique is adapted to machine use. On a C.P.C. it has solved a set of 
four simultaneous non-linear equations which describe a hypothetical one-dimen- 
sional free-radical flame [1]; on an IBM 701 it has integrated the seven simultane- 
ous non-linear equations for a hydrogen-bromine flame. 

2. Character of Equation Systems for Which the Method may be Useful. 
Certain physical problems such as free-radical flames are described by sets of 
simultaneous differential equations in which the calculation of one of the deriva- 
| tives from its analytic formula involves serious loss of numerical accuracy due to 
subtraction. The conventional solution of such problems presumes the steady 
| state hypothesis which replaces the differential equation by the algebraic approxi- 
mation that the derivative vanishes identically. Unfortunately, questions such 
as the importance of free-radical diffusion in flames must be studied without the 
| ad hoc assumption of the steady state. (When applied to free-radical flames, this 
approximation states that the derivative, dG;/dt, of the fractional mass rate of 
flow of the free-radical B vanishes identically. Such an assumption must not be 
introduced but should be proved valid if we seek to establish the importance of 
' diffusion.) Furthermore, it is generally important to determine the accuracy of 
such an approximation. 

An outline of the qualitative features of one type of equation system which 
presents this numerical problem will clarify the essential character of the difficulty 
One meets when he applies traditional point-by-point procedures without assum- 
ing the steady state. Consider, therefore, a set of M simultaneous differential 
equations in an independent variable Z, 


(1) du(Z)/dZ = Fi, 1 Sj < M, 


'where in general the F; may be functions of any or all members of the set of 
dependent variables {u;} and, or, the independent variable Z. Suppose that the 
set of functions { F ;} has four qualitative characteristics which will be enumerated 
‘presently. The first two of these four conditions state that the first function, F;, 

is poorly determined numerically. As subsequent discussion explains, the com- 
bination of the four conditions describes a coupling which should cause traditonal 

"point by point procedures of integration to fail and which suggests the particular 
Successive approximation schema outlined in Section III. 
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The conditions are: 
(2) F, = pw Q1,ky F, = > 2,k- 
& a 


(a) F, and F, are each given as a sum of terms. 

(b) F; is smaller than one of the terms of equation (2), say, a1,:, by a factor 
of 10”: Fi/a;,, = 10-”, p a positive integer. 

(c) Calculation of an s digit value for one of the terms of F:, say, d2,:, requires 
that u, be known to s digits. Moreover, d2,; and F, are both of the same numerical 
magnitude. ; 

(d) An s digit value for the term a;,, numerically determines an s digit value 
for ue (e.g., 21,1 might be given analytically as a;,; = u2f where f is a function of 
other variables whose value can be computed accurately to s digits). 

Since the conventional methods of computing the increments 


(3) Au (Z,) - u (Zn) can u (Zn-1) 


in the course of a point by point numerical integration require calculating the 
first derivatives, F;, it follows from (a) and (b) that Fi[{u:}, Z,] and therefore 
Ai%(Z,) can be computed to only (s — p) digits whenever the {u,(Z,,)} are known 
to s digits. (This discussion assumes that an implicit method is being used in the 
integration. If an explicit method such as the Runge-Kutta were used, the se- 
quence of subscripts on Z would be different, although the principles would be 
the same.) Furthermore, the error in A,#;(Z,) appears in u;(Z,) and, therefore, 
by condition (c), in F.{{u;}, Z,]. Unfortunately this introduces error first into 
A,u2(Z,) and thence into u2(Z,). The completion of a vicious cycle through which 
this error reappears in the calculation of F,[{u:}, Z,] follows from conditions 
(b) and (qd). 

If the desired solution of the equation system actually does maintain condi- 
tions (b), (c), and (d) over an appreciable interval in Z, then these conventional 
methods might be expected to fail completely. The feedback of error might be 
expected to destroy quickly the approximate equality 


(4) Fi{{u}, Z] = 0. 


In the case of the flame equations for one hypothetical free-radical system, 
the build-up of error was so rapid that in eight digit calculations the solution was 
lost in as few as six integration steps. (See Appendix A of reference [1] for an 
illustration of the failure of conventional techniques. A heuristic interpretation 
of this failure concludes the appendix.) 

3. Proposed Iterative Solution. 

A. General basis of the method. If the analytic formula for the second derivative 
of uy 


(5) d*u;(Z)/dZ* = dF,/dZ 


did not present the same computational problem, the equation system could pre- 
sumably be solved by the simple expedient of adopting F; as a variable and adding 
the differential equation for F; to the system. However, if the approximation 
F, = 0 holds for any appreciable range, then the calculation of dF,/dZ might 
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be expected to exhibit the same numerical indeterminacy as the calculation of F; 
itself. This is the case in several flame problems. 

The very failure of an attempt to circumvent the difficulty by adding the 
equation for dF,/dZ to the set suggests that perhaps the approximation F,; = 0 
is maintained over a considerable range in Z. In this case the four conditions 
(a)—-(d) of Section 2 suggest that the solution to the complete set of equations 
might be approached by integrating the following sequence of equation systems. 
The kth approximating set of equations will make use of (k — 1)st order approxi- 
mations for u; and du;/dt. For k = 1, the steady state approximation, 


(6) Fi[{u:}, Z] = 0, 


is made. According to conditions (b) and (d), equation (6) gives a numerically 
determinate implicit formula for uz. 

An approximation for ; can be obtained as follows. Recall that the hy- 
pothesis has been made that the analytic expression for dF,/dZ exhibits the same 
numerical indeterminacy as F;. Since, moreover, “2 is supposed to make a nu- 
merically important contribution to a:,:, which in turn has been supposed to be 
one of the most important terms of Fj, it is reasonable to hope that the term 


(7) (0a; ,1/OuU2)(du2/dZ) = (0a1,1/Ou2) Fe 


will make a numerically important contribution to dF,/dZ and that a value for 
dF,/dZ will determine without significant loss in accuracy a value for F,. Now, 
according to condition (c), an s digit value for F: is supposed to determine an s 
digit value for u;. Thus the approximation 


(8) dF,/dZ = 0 


may provide a numerically determinate implicit equation for F: and therefore 
for u,;. These approximations are valid in some free-radical flame systems. 
B. Outline of the steps in the iteration scheme. 
1. Step fork = 1. 
Solve the first set of approximate equations: 


(9a) F,[{u:}, Z] = 0 (an implicit algebraic equation for 42“) 
(9b) (dF,/dZ) = 0 (an implicit algebraic equation for «,) 
(9c) (du(Z)/dZ)= Fi, 3555 M. 


2. Step for (k + 1). 
Assume that the u;“, 2 < 7 < M have been computed. 
(a) Numerically differentiate u2“) to obtain (du2/dZ®. 
(b) Since the u;“, 2 <7 < Mare known, calculate u,“ from the equation 


(10) (duz/dZ) = F.[{u™}, Z). 


(Recall that condition (c) states that equation (10) is a numerically 
determinate formula for 1.) 

(c) Numerically differentiate u;“ to obtain (du,/dZ). 

(d) Use the kth approximations which have just been computed for u,“) 
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and (du;/dZ) to integrate numerically the (k + 1)st set of approxi- 
mate equations: 


(11a) (duy/4Z)™ = Fil{us™, uy}, Z] 
(an implicit algebraic equation for u2“+) 
(11b) (du;/dZ)* = Fi{{u,™, u,**}, Z] 
3<jSM, 2<51< M). 


A detailed illustration of the exact way in which one particular set of equations 
falls into the proposed successive approximation schema is given in Appendix C 
of reference [1]. The convergence is illustrated in the following graph. 


0.10 -- 


0.08-- 
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0.02-- 
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Convergence of #2: The ordinate is the fractional difference between the second and first approxi- 
mations [12 — u‘)]/u,) for the set of four simultaneous equations which describe one hypo- 
thetical free-radical flame. The maximum fractional error is always less than thirty percent and 
reaches 10 percent only after uz has dropped to one tenthousandth its initial value.. The convergence 
of %# is equally good. 





4. Conclusion. If the approximation schema converges to the solution (in the 
case of the two equation systems studied, the schema appears to converge to the 
correct answer), then it possesses two noteworthy advantages: (1) it succeeds 
when other point-by-point procedures fail because of numerical indeterminacy 
of one or more of the derivatives; (2) it can be readily programmed for machine 
calculation. Compared with the relaxation technique which also has been found 
to be stable towards this type of numerical indeterminacy, successive approxima- 
tion requires more calculation. However, the relaxation method requires repeated 
use of judgment which can become difficult if the variables in the equations are 
badly cross-linked. Although by far the greater part of the calculational labor in a 
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relaxation treatment can be programmed for machine, that part which requires 
judgment would be difficult to code. 

The disadvantage of successive approximation lies in increasing the length of 
the calculation. However, as the numerical indeterminacy grows more stringent, 
the initial approximation improves and the number of iterations required de- 
creases. Moreover, this is a comparatively less serious disadvantage for a method 
suitable for machine use. Thus the number of times the calculation is to be re- 
peated for a different set of input data, the character of the equations with re- 
spect to cross-linking, and the comparative labor of constructing the machine 
programs must all be considered when deciding which method will prove most 
convenient for any particular study. 

Epwin S. CAMPBELL 


Department of Chemistry 
New York University 
New York, New York 


1. Epwin S. Campsett, CM-847, University of Wisconsin Naval Research Laboratory, 
July 11, 1955. 


A Smoothest Curve Approximation 


A practical problem which often comes up in numerical work is the fitting of 
a curve to a finite set of known values in order to perform various operations such 
as integration. The usual method of approximation consists of fitting the points 
with one or more polynomials (independent of each other). By letting there be 
more points for each polynomial, with the polynomials being of comparable order, 
the error of approximation becomes asymptotic to a higher power of the interval 
length. 

However, error analyses for such methods usually depend upon the bounded- 
ness of some derivative of a correspondingly high order [1]. But even if the func- 
tion to be approximated is analytic, its correspondingly high order derivative may 
be of sufficient magnitude that for the given interval size, a simpler method would 
give better results. For instance, a fitting with an eighth order polynomial gives 
the following rule: 


. 8 
[sae = 38350 [989 f(0) + 5888 f(1) — 928 f(2) + 10496 f(3) 


— 4540 f(4) + 10496 f(5) — 928f(6) + 5888 f(7) + 989 f(8) ]. 


An application of this formula to a positive function which was everywhere small 
over the range, apart from a sharp peak in the center, would lead to a negative 
result. Try to convince a prospector that there is a negative amount of mineral on 
his land because he finds a rich strike in the middle of it! 

Frequently one contents oneself with a simpler rule which is repeated in blocks 
of so many intervals per block. However, this usually introduces discontinuities in 
the first derivative at the junctions of the blocks. If one were to integrate an in- 
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terval broken up into 1000 equal subintervals, why should the 500th point be given 
twice as big a weighting as the 50ist point (the middle point), as would be the 
case using Simpson’s Rule? While this may not be a great objection for functions 
which do not undergo great changes of nature over several intervals, it may be 
quite unfortunate for those that do. 

The main problem to be considered in this paper is the problem of obtaining 
and examining an integration rule based on integrating a function passing through 
a given set of points such that the function will have a small amount of twisting, 
and such that whatever twisting is necessary will be spread out. To be more speci- 


fic, we will minimize 
a? 2 
eae 


where f ranges over an appropriate set of fitting functions, say those which are 
continuous and have continuous first and second derivatives. 

It would be improper to consider an “‘angle-like’” function, such as f)(t) = 
#(0 <¢ < 1) andf,(t) = 1 (1 < ¢ < 2), and argue that since fy” (#) = 0 al- 
most everywhere, 


Ii Cfo’ (t) Pdt = 0. 


In fact if {f;} (¢ = 1, 2,3, ---) is a sequence of appropriate functions which 
converge to fo, then 


2 
lim Cf’ (t) Fadt = OO, 
ima J 0 


Another application of this problem is that of bending a stick a little bit. A 
bent stick assumes a position which, subject to the bending constraints, will 
minimize its potential energy. Assuming the stick to be originally straight and 
uniform, Hooke’s Law implies that its potential energy will be proportional to the 
integral along its arc length of the square of its curvature. If the problem is two 
dimensional, and if the bending is sufficiently slight that the arc length may be 
considered as being practically proportional to some coordinate axis, then we get 


the problem of minimizing 
a? 2 ‘ 


In considering this stick application, one may also obtain the function given in 
Theorem 1 by mechanically balancing torques against curvatures. 


1. THEOREM. Let to < t; < +++ < ty. For eachi = 0, ---,m — 1, let f bea 
cubic on the interval [1;, ti+1 ]. Also let 
a? v7 vt 
ap (*0) =f" (to) = 0 = f" (tn). 


Fort =1,---, m — 1, letf(ts—) = f(tit), f'i—) = f' it) andf'(ti—) = 
f'' (tit). Let g(r) be any admissible function on [to, ta] such that g(t;) = f(t) for 
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4 = 0,---,n. Then 


[ [e” (t) P dt > [ Cf” (t) P dt 
if and only if g ¥ f. 
Proof: Define n(t) as g(t) — f(t). Then since f’”’ is constant in each interval 


* Ce” (t) Pat — PR Cr’ Fat = ‘ ("Far + 2 ig a (Of" (tat 
Y [. Ca" (Pat + 2S La’ tend" Uten) — a COP" 


—s = | oe L'a! (t)dt = [. "En!" (t) Pat + 2En (te) f" (tn) 


i=0 


— 9 (to)f" (te) — 2° f'" th interval) [a (tiss) — 2(¢0)] 
i=0 


o [owrora +2[0 -0]-2 50 


10 
which is positive for admissible non-zero ’s. 

Note that were ¢, and #, interior to the end points of the interval on which f 
is defined, then for minimizing f [f’ (t) dt, f would be the same as before on 
[to, tn]. Also, f would be linear at the left of tp and the right of ¢,, and the first (as 
well as the second) derivatives would match at ty and fy. 

The function f of Theorem 1 consists of ” cubics, and each cubic contains four 
coefficients. The restriction that f passes through f(t;) at t; gives 2n determining 
equations. The conditions f’(t;—) = f’(t:+) and 

Ff" (ti-—) = f'" (t+) (¢ =1,---,m — 1) 
give m — 1 equations each. The equations f”(t.) = 0 and f’(t,) = 0 give still 
two more, giving a total of 4n. Therefore we have 4n linear equations to deter- 
mine the 4m coefficients of f. 

2. Definition. In Theorem 1, let ¢; = 1 (0 < 4 < m). Then for m an integer 


(0 < m < n), define w(m, nm) as [[roa where f(z) = én, (¢ = 0, --- ,m). 


The following is a table of w’s form < 10. In this table w(m, n) isthe m + 1-st 
entry of line . Also the common denominator is factored out at the left. 
¥ [1; 1] 
¥% [3; 10; 3] 
Mo [4; 11; 11; 4] 
Ms (11; 32; 26; 32; 11] 
gg (15; 43; 37; 37; 43; 15] 
Mos (41; 118; 100; 106; 100; 118; 41 ] 
4 42 (56; 161; 137; 143; 143; 137; 161; 56] 
Mss (153; 440; 374; 392; 386; 392; 374; 440; 153] 
iso (209; 601; 511; 535; 529; 529; 535; 511; 601; 209] 
VY gag (571; 1642; 1396; 1462; 1444; 1450; 1444; 1462; 1396; 1642; 571]. 
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Since f (and therefore f fdt) is linearly dependent on the f(#;), we get the 
integration rule 4 fd = a f(m)w(m, n). For instance, from the third 
m=(0 
line of the table, 


[2 seoae = Fano + 1scy + 119) + 471 


Notice that as max [m; n — m] gets large, | w(m, n) — 1| gets small. 

The following rules which follow from Theorems 6 and 7 are probably the most 
useful for writing down the w’s. Define D(m) as the common denominator for the 
w(m, n)’s. Since 


no =a (RF + S25F 





and 


pon en avi (Myr + (55) 


it follows that D(2n + 2) = D(2n) + 2D(2n + 1) and 
D(2n + 3) = D(2n +1) + D(2n 4+ 2). 
A similar rule also holds for the numerators. For 
m < 2n, D(2n + 2)w(m, 2n + 2) = D(2n)w(m, 2n) 
+ 2D(2n + 1)w(m, 2n + 1), 


and for m < 2n + 1, D(2n + 3)w(m, 2n + 3) = D(2n + 1)w(m, 2n + 1) + 
D(2n + 2)w(m, 2n + 2). For 

D(|\2m—n ) 

0<m <n,w(m,n) = 1+ (-1)™ 


2D (n) 





where D(0) = 4. Also note the obvious fact that w(m, n) = w(n — m,n). 
3. Definition. For n a positive integer, consider a function f on [0, 2] such that 


f isa cubic on [7,4 + 1] (¢ = 0, ---,n — 1); f(t) = 0 (¢ = 0,--- , mn); 
f'(0) #0; f"(n) = 0; 
and fort = 1,--- ,n—1,f’(i—) =f’ (t+) andf”’(i—) = f’ (t+). Then define 


R(n) as f’’ (0) /f’(0). Define I(m) as {> f(t)dt/f’ (0). Also define R(O) = (0) = 0. 


Define x as the number 2 + +/3. Define y as itys = /x. 


4. LemMa. For n a non-negative integer, R(n) = -2/3= rt = , 


oo 0 
Proof: R(0) = 0 = —2+1/3 pe . To use mathematical induction, let us 


show that the validity of this formula for R(”) implies its validity for R(m + 1). 
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Let f be the function of Definition 3 defined on [0, m + 1]. f(0) = f(1) = 0 and 
f'(0) # 0 imply that there exist real numbers a and } + 0 such that f(t) = 
Lt — af? + (a — 1)#*] for 0 < ¢ < 1. Then 


_f'() _ 44-6 





R(n) 





Ff) a@=2 

— _ 6 — 2R(n) 
implies that a = ie Therefore 

. 2 oe E _ Rio)) 

R(n +1) = F(0) = —2a 4 t— Rin) 
which, substituting for R(m), equals 
gett — x (atl) 
—2V3 a eee 
i (—1)* 


5. LEMMA. For n a non-negative integer, I(n) = oh 6a" + x*)* 


9 ORE. > =~ (—1)° : : : 
Proof. I(0) =0 = 12 ~ 6° 4+ x) ° To use mathematical induction, let 


us show that the validity of this formula for J(m) implies its validity for I(m + 1). 
Let f, a, and \ be the same as in the proof of Lemma 4. Then 


In +1) = atx [soa = BT [soar + rare | 


a 1 R(n +1) 1 
m4 = Gt 1 =} - FEE [rem - 3] -21 


| 


- P = ( - 1) n+1 
12 6[«"+ + x (at)) 7 





on substituting for R(m + 1) and I(m). 
6. THEOREM. For n a positive integer 
mtl (-— )-@+) 
w(0,2) = y y . 
a) 2V6Ly* + (—y)™] 
Proof. Let f be the f of Definition 2. Then f(0) = 1, f(1) = Oandf”(0) =0 
imply that there exists a real number a such that 


f® =1—at+ (a — 1)#for0 <t# <1. 


_ f"(1) _ 6(a — 1) _ 3R(n — 1) — 6 
3 — R(m — 1) 
4 — R(n — 1) 








. Using 





R(n) = —4 
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3 + R(n) $3 6 
t+ Rin) andsoag = 3+ Rw) ° 


w(0, n) = [Psa = [sma + y'(1)I(n — 1) 


from Lemma 4, we see that R(m — 1) = 4 











a = Need (—1)=" 
oa 4+ “3 FE aes | 
3). e. -26=~-3  {—8)" 

“3-i2+~ 6 =t+a= 

i 

i 


1 1 2 (—1)” 
~ 2R(n) + E + na |= +x° 


On substituting for R(m) from lemma 4 and after some manipulation, we obtain 


™ nth _ (—y)—@+) 
0") = 97a + (1 





7. THEOREM. For 
y= + (—y)?"— 
ay + (—y)* J 

Proof. Let f be the function of Theorem 1 where t; = i — m (i = 0, --- , nm) 
and f(#;) = dm, Then f(0) = 1,f(—1) = f(+1) = 0, f'(0—) = f’(0+) and 
f’(O—) = f’(0+) imply there exist real numbers a and 6 so that f(#) = 1 + 
at + b#@ — (1+a+5))#for0 <¢<landf(t) =1+at+ d+ (1 -—a-+ d)# 
for —1 < ¢ < 0. Define a = m — m — 1 and B = m — 1. Then 


R(a) = S71) _ 6 + 60 + 4 


fi) 3+20a+5 


0<m<n,w(m,n) =1— (-—1)” 








and 
_f"(=1) _ 6 - 6a + 4% 
A-) > d= Bee 
Solving these last two equations, we obtain 
3R(a) — 3R(B) 


R(6) = 











* ~ 2R(a)R(B) — 7R(a) — 7R(B) + 24 
and 
cehe —6R(a) — 6R(8) + 36 
2R(a)R(8) — 7R(a) — 7R(B) + 24° 
Now, 


ll 


w(m,n) —1 


i a" f(t)dt — 1 


Has’) — 1@s'(—1) + f* sae - 1 


= I(a)[—2a — (b + 3)} + 1(8)[2a — (b + 3)] + i (b + 3). 


Substituting the values for a and 6 + 3 into the last equation and multiplying 
both sides by the common denominator of a and b + 3, we get 


[2R(a)R(8) — 7R(a) — 7R(B) + 24] [w(m, n) — 1] 


= 12I(a)[R(8) — 3] + 127(6)[R(a) — 3] — R(a) — R(8) + 6. 
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Now substitute for R(a), R(8), I(a) and I(8) the expressions of lemmas 4 
and 5 and multiply both sides of the equation by (x* + x~*) (x + x”). 
Then collect like terms to get 


[(48 + 284/3)xe + (48 — 284/3)x-“<*] [w(m, n) — 1] 
= (6 + 4y/3) (—1)*x* + (6 + 40/3) (—1) *x* 
+ (6 — 4/3) (—1)%* + (6 — 473) (—1)e-*. 
Divide both sides of this equation by 2+/3 and substitute x = 2 + +/3 to get 
2[xatee® — x—(at8+2)")] [w(m, n) — 1] 
= (—1)Peott + (—1) ax8H — (—1)8x-Cet — (—1) «x GHD), 
Since x = y*, we may divide both sides of this equation by 
[svete — (—y) Coton] [yaroes + (—y)“Coten] 

which using the binomial expansion, equals 


<< 2;k1B; y i— (—x)-@+) 
jk —jt+i1)mi  nvV6 1+ (-x)™ 











k n n-* (—2m) ‘k! (—x)— 
p> a'(k —2)! 1+ (—x)™ 


_+ 21h! eee a 
- DP a+ (aS | Ae 


aia 
n 











i=l m=0 
1 hit da‘ ; a 
- eee PE + (-1) £ 7(-1)| 
+1 a ao 
[Po — a toad | wiipcibe: |. 


8. Definition. For i — 1 a positive integer, define B; as i factorial times the co- 
efficient of t* in the power series expansion of t/(e* — 1). In other words, let 0 = 
B; = B; = B, = --- , and let Be, By, Bs, Bs, --- be the Bernoulli’s Numbers 


i Fee a fe 
6’ 630’ 42’ 30’ 66’ = 2730’ 6’ 


9. THEOREM. Let P(t) be a finite polynomial on the interval (a, b). Then 


[ P(t)dt = 3° P(a + mh)w(m, n) h 


ht di di 
7 Sem sro + (-1'£ P@)| 


t+1 3 ‘ —m 
[Ban ~ Pisa & (-™)(-#) | 


m=0 








where h is the sub-interval length = °— ©” 
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Proof: If this theorem holds for a given a < 5, then for \ a positive real num- 
ber, this theorem also holds for end points \a and NJ, since new h = } old h and 
d‘ — © ot 
mi? att =a(orb) =X rT P x at ¢ = da (or AD). Also, the three terms of 
this equation are invariant under translations. Therefore, it is sufficient to prove 
this theorem for —a = 6 = 1. Also, since this equation is linear in P, it is sufficient 
to verify it for P(t) = ¢* where k is a non-negative integer. 

If k is odd, then P(t) = —P(—t) implies that all three terms of the equation 
of this theorem are zero. Therefore, we may let k be even. 


Ap P(—1 + mb) w(m, 2) hb — p Pat] 


m=0 


-[E(7@7) om |- a 

















ott yi, - ¢- “2a + (=) ’ \— 1 
[1 - mE ag] + 200 } 
MS geek. 2 ene eet Oo 7 = ‘ 
“HS(0s2Y|i—ceorrsfeg]-14 mn pty 
° 1 <—/n — 2m\ 22 ri \F 
If mis even, + 5 ( - ) -23 (2). 








n B. k (n—1)/ 2 _ k 
Me teedd, - >(* 2m) =; (4+ty. 


m=0 i=0 n 


Using the Euler-Maclaurin sum formula [2], either of these last two expressions 


equals 
1 1 2ikIB; 
ye ets DS AEE fyi ni’ 





Substituting, we get 
1 2 2ik1B; 
| Lew — fre | = 2 ee a 


2w(0, n) 1a(/n—2m\ (-—x)™ 
700.7) - =X ( n or 
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D(2n +1) = VY2[y™ + (—y) Ot], 
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[- Diy] = - ¥3. 


i 
[4 -}0)| = 0. 


Aa 
i<y 
. w 


S| 
le 


11D(1) — D(3)] = = , 


3 
ond 


> 
ome i ™= OE OE 


} os 

sie RIS Rie &) sie Rl 
i 
IS 
os 


i) 
_ 
an 


> 
a 


—66 + ; [26D(2) — Dia)}} = 


I 

5 

— 

an 
eo! — 
Se 


6s( - Bi [-302D(1) + 57D(3) — D(s)] = — >. 

6: = 95 {2416 + ; [-1191D(2) + 120D(4) — Do) } —_ >. 

bs (2) ~ Rs. [15,619D(1) — 4293D(3) + 247D(5) — D(7)] = - tye. 
s (*) = aig {— 156.190 + 5 [88,234D(2) — 14,608D(4) + 502D(6) 


7 
- Dis)}} =5. 
Therefore, the equation 


b n hit 
[ P(t)dt = T(m) + x P(a + mh)w(m, n)h — > G+iI 


dé , di ae: 
|Z Pe) + (-—1) < PCa) | Ber (—1) —z (*)]- 
gives us the series of this theorem. 


10. CoroLLarY. >> w(m, n) = n. 


m=0 


11. THrorem. Let P(t) be a finite polynomial on the interval (a, b). Then 
b n 
I P(t)dt =T(n) + >> P(a + mh)w(m, n)h 
a m=0 


v3 3 ii ii 1 4 iii _— iii 
— “aq LPM) + P8(a)] + 35 ht LPH) — PH(@)] 


V3 ye (Pb) + Peay] — 1 ne [Pr(b) — Pr(a)] 


+ 364 2016 





til mee h" (Pri(b) + Pri(a)] + stew [Prii(b) — Prii(g)] 


_ _1ivV3 
4,354,560 


31 


h® (Priti(b) + Priti(a)] — 9.580.032 


h® [P*(b) — P*(a)] 


+ +++, whereh wm = ng tien (2 + 9/3)"T (2) = 0. 
n>o 


n 
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Proof: For each i, > (—m)*(—x)-"is an alternating series. Therefore, for 


m=(0 


sufficiently large n, 








x (- m) *(— -2)-" - ¥ (- m) *(—x)-™ 


< (n + 1) *¢-(+1), 
Since h*" is proportional to 1/n*, 


lim of [Poa ~ > P(a + mh)w(m, n)b 


n?>o 


ht 
+7 [|Z7m + Zr@| 


| Bas -#*y (—m)«(—x)-]} ="0. 





If we define ¢;(¢ as) > m‘(—t)™, then diui(t) = . o:(t). 
m=0 


—t+f 


Since oo (t) = = »oi(t) = a+ 4. t)* ; 


1+? a 3 O2(t) = 


TF 


—t+4e-—# 
ten 3 ba(t) = 


—t+ lif — 11°44. 
(1 + #)° 





o3(t) = 


—t + 268 — 66t? + 26¢ — t 
(1 + #)§ : 





$5(t) - 


—t + 572 — 302¢ + 302 — 57¢5 + #8 
(1 + ¢)? 





os(t) = 


etc., and keeping in mind that 


1+ 


Sle 


6 6 = 2n in —2n 
V8 4/8, Diam) = 20 + ( »° ] 
to get 


gerttacg) OO 4. Ptt( a) ~iet® 
yt + (—y) eto] 
yr + (—y)*"-" 
2[y" + (—y)] 
12. Examples. Let g(t) be a function on [—1, 1] and let us consider the error 


> g(-1 + 2) 2 (m2) _ p g(t) dt. If g is constant or g(t) = —g(—2), 
m=0 -1 


then the error is zero, as is true of most rules. 





w(m,n) —1= — 


= (-1)" 
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If we let m = 10 and consider g(t) as /, ¢ and ¢* respectively, then the error is 
7.7 X 10-*, 4.5 &K 10-*, and 1.1 X 10 respectively. If we now apply the end 
point terms of Theorem 11, the remaining errors are —2.9 K 10-*, —1.8 x 10-* 
and —4.2 X 10-® respectively. If we let m = 14, then the above quoted errors 
become 2.8 X 10-4, 1.6 X 10-*, and 4.0 X 10-* and also —5.5 X 10-", —3.2 X 
10-" and —8.1 X 10-". The corresponding errors using Simpson’s Rule are 0, 
4.3 X 10-* and 2.1 X 10-* for 10 intervals and 0, 1.1 K 10~ and 5.5 k 10~ for 
14 intervals. 

If we let a = cos (at/4) and m = 12, then the error is —1.7 XK 10-*. If we 
apply the _ h* [g’"(1) — g’’(—1)] correction term, the remaining error is 
—7.1 X 10-’. The corresponding error using Simpson’s Rule is 6.7  10-*. 

13. Semi-Infinite Interval. Let us define w(m, ©) as lim w(m, nm). Then 


w(0, 0) = % 


and 


ae aye eS (—x)-™ for m > 0. 


Given a function g(t) for which the ving’ g(t)dt isconvergent, une may 
0 


approximate this integral by the expression h > *_, w(m, «)g(mh) where N is 
sufficiently large that h }°@_y41 w(m, ©)g(mh) is insignificant. If one also knows 
the first few derivatives of g(¢) att = 0, one may also add the terms obtained from 
Theorem 11. 

If g(t) = e~*, then the series 


1 


h i w(m, oo )g(m, h) " ys h®g"(0) — a0 ’”’ (0) 
v3 h®g “nn 6/ttt 
+ 364 (0) + a h'g""(0) + - 


will converge to / g(t) dt whenever hd < 2r. 
0 


For some functions such as g(t) = e~ # this sequence of derivatives is divergent 
for any h > 0. However, this fact does not imply that using this formula with the 
first few derivative terms would not do a good job of approximation. 

Notice that the weight terms w(m, ) are easy to compute on a machine, 
requiring only one multiplication per term. 
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Numerical Integration over the Spherical Shell 


1. Introduction. Numerical integration over three-dimensional regions has re- 
ceived relatively little attention. The papers listed in the bibliography give the 
work on this subject known to the author. With regards the sphere, Hammer and 
Wymore [6], have given formulas with specific polynomial accuracy. This paper 
gives a formula for arbitrarily high polynomial accuracy and generalizes the re- 
gion to the spherical shell. 

2. Type of formula. The type of formula considered is a weighted sum of 
integrand values 


[ff tee, dnayas =X aes 9. 29, 
i=1 

where the dot over the equals sign is meant to signify that in general the right 
hand side will only be an approximation to the left hand side. The region under 
consideration is the spherical shell of inner radius R and outer radius one. By a 
theorem of Hammer and Wymore [7], a formula over any region which is an affine 
transformation of this spherical shell can be easily obtained. The criterion used 
to determine the weights, a;, and the points, (x;, yi, 2;), is that the formula shall 
be exact whenever the integrand function is a polynomial in x, y, and z of degree 
at most s. 

A numerical integration formula, whether it be in one, two, or more variables, 
is said to have accuracy s provided the formula is exact whenever the integrand 
function is a polynomial of degree s or less, and provided there exists at least one 
polynomial of degree s + 1 for which the formula is not exact. 

3. Numerical integration over the spherical shell. We first transform the 
problem to one in spherical coordinates. Under the transformation 


x = rsin ¢ cos 0 
y =rsin ¢ sin 0 
2=>rcos¢ 


the integral | | f(x, y, 2)dxdydz over the spherical shell with inner radius R 


and outer radius one can be written as 


i i ‘ ; *" 2 sin $F(0, ¢, r)d0dedr, 
R/j0 0 


where F(6, ¢,r) = f(r sin ¢ cos 6, r sin @ sin 8, r cos @). Since any formula for the 
integral in rectangular form is to be exact whenever the integrand function is a 
polynomial of degree at most s in x, y, and z, any formula for the integral in 
spherical form must be exact whenever F(6, ¢, 7) is a polynomial of degree at most 
sinr sin ¢cos 6,7 sin ¢ sin 6, and r cos ¢. The type of formula we attempt to find 
is 


(i) I= i, i i i "sin GF (0, 6, Ndddgdr + DOD Dink Os 61.1) 
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which is a weighted sum of integrand values where the limits on 4, j, and & are to 
be determined. We now state the main result of this paper in the following 
theorem. 

THEOREM. [f it be required that 


(a) formula (1) have accuracy s = 4m + 3, m = 0,1, 2, ---, in r sin ¢ cos 6, 
r sin ¢ sin 6, and r cos ¢, 

(b) the evaluation points be taken at all intersection points of the planes @ = 6;, 
the planes @ = $;, and the spheres r = fr, 

(c) the evaluation points under this configuration be minimum in number, 


then tt ts both necessary and sufficient for the existence of a unique set of real weights 
D; ik that 
(a) 6; sa 2x ;/(s + 1), 7 -” i, 2, 2e2e + he 
(b) the cos $; be the 2m + 2 zeros of the Legendre polynomial Pims2 of degree 
2m + 2, orthogonalized on {|—1, 1], 
(c) ther, bethem + 1 zeros of the polynomial in r* of degree m + 1, Qmss(r?), 
where J 1Qmsi(??)Tn(r*)dr? = 2 fi P’Omsi(?*)Ta(r*)dr = 0, and where 
T,,.(r*) is an arbitrary polynomial in r* of degree m or less. 


The unique weights D; ;, are equal to A,B ;C; where 








2x ‘ 
haw ae = 1,2, ---, 1, 
3 ee ai s+ 
1 ’ Pom42(y) » 
By - yy) ——— dy, = 1,2, ---, 2m + 2, 
7 Pom+2(y j) a y aan Vi y J 
= i * ?7Om+1(1?) ss 
C= ates [Sean = 1,2, m+, 
or 
1 1 
CG. = Vt Omi) ay where t = #2 


~— 20 aa (te) Jas t — be 


In the first expression for C;, the notation Q/,,,(r,2) indicates a derivative with 
respect to 7. 


Proof. Because of the arrangement of the points we have chosen, we write the 
left hand part of (1) as 


1 na 2r 
I= i ef sin +f F(0, , r)d@d¢dr. 


H(, 1) = i ” £6, 6, r)d8, 


Now define 


G(r) = [ " sin 6H(¢, r)de, 
which gives 


rhe [roar 











246 NUMERICAL INTEGRATION OVER SPHERICAL SHELL 


We now look for three formulas, namely 


(2) H(¢, r) a i F(@, ¢, r)dé = YA F(6,, d, r), 
(3) G(r) = I sin oH(¢, r)do # > B ;H(¢;, r), 
(4) hie fi YG(r)dr = Y C,G(r). 


When these three formulas have been obtained, we will combine them to obtain 
formula (1) with D;;, = A<B;C,. 
Formula (2) 


H(¢, r) = ig F(0, ¢, r)dé@ = > A;F(6;, , 17). 


The number ? is to be determined, and this formula must be exact whenever 
F(@, ¢, r) is a polynomial of degree s in r sin ¢ cos @, r sin @ sin 6, and 7 cos ¢. 
Since r and ¢ are parameters here, we can require formula (2) to be exact whenever 
F is a polynomial of degree s in cos @ and sin 6, and hence whenever F = 1, 
cos 6, sin 6, cos 20, sin 20, ---, cos s@, sin s@. This leads to a system of 2s + 1 

equations 


2r 
Art Aste +4, = [ 1d0 = 2x 


II 
o 


2r 
A, cos t@; + Az cos 0, + --- + A, cos t@, = [ cos t6d6 
0 


II 
a) 


2r 
A; sin #0; + Az sin t@, + --- + A, sin £6, -.[ sin t6d0 
0 


for? = 1,2, ---,s. 
One solution to this system (a sufficient condition that formula (2) have ac- 
curacy s with minimum ?) is given by 








p=s+i 
2x * 
| ey r ’ = 1,2,-:--, 1 
A as 1 s+ 
ee oo See ee 
é +1" » 4) ’ . 


The requirement that the A; be real leads to the conclusion that this is the only 

solution up to a rotation. (For the idea which led to this result, I am grateful to 

Professor Fritz Herzog.) The above is therefore a necessary and sufficient condi- 

tion that formula (2) have accuracy s in r sin ¢ cos @, r sin ¢ sin 0, and 7 cos ¢. 
Formula (3) 


Gin) = l "sin 6H(6, do + Y BGs 1) 


The number g has to be determined, and this formula, like (2), must be exact 
whenever F(@, ¢, r) is a polynomial of degree s in r sin ¢ cos 0, 7 sin ¢ sin 0, and 
r cos ¢. Considering the definition of H(¢, r), formula (3) must therefore be exact 
whenever H(¢, r) is 2 polynomial of degree s in cos ¢. This leads to a system of 
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s + 1 equations 
B, cos ‘¢; + Bz cos ‘d2 + --- + B, cos ‘¢, = I’ sin ¢ cos ‘¢d¢d 
for? = 0,1, 2, --- , s. Putting cos ¢ = y, these equations become 
Byyi' + Bayo' + --- + Bay,' = [i v. 


The unique solution to this system for minimum g has been given by Gauss [4] 
and is 








= pe 2m + 2, 
¥i = cos ¢; = the g = 2m + 2 zeros of the Legendre polynomial Po»: of 
degree 2m + 2, orthogonalized on [—1, 1]. 
B;= 1 . Pamtx(¥) gy 


— Pomgelys) Jayoys 7’ 


where P’ means the derivative of P with respect to y. 
Formula (4) 


fa / '26(Ndr = > G(r). 


Again, the number v has to be determined, and this formula, like (2) and (3), must 
be exact whenever F(6, ¢, r) is a polynomial of degree s in 7 sin ¢ cos 6, r sin @ 
sin 6, and 7 cos ¢. From the definition of G(r), formula (4) must be exact when- 
ever G(r) is a polynomial in r* of degree (s — 1)/2 = 2m + 1. Putting g(r*) = 
G(r), formula (4) can be written 


2I = / ‘ eg(P det = > 2Cig(r2). 
R k=l 


The unique solution for the C, and the r,; for minimum 2 is given by Szegé [11] 
and is 


v=m-+ i, 
12 = the m + 1 zeros of the polynomial in r* of degree m + 1, Qn4a(r*), 


where I * 1Omsi(?*)Ta(r?)dr? = 2 / ‘°Qns(?*)Ta(t*)dr = 0, and 


where T,,(r*) is an arbitrary polynomial in r* of degree m or less, 


1 Omer(?4) 
— one I. fom 


where Q’ means the derivative of Q with respect to r*. This completes the proof of 
the theorem. 


Several corollaries are immediate. 

CoroLiary. The weights D;;, = A,B,C, are positive, and their sum is the 
volume of the spherical shell. 

COROLLARY. All evaluation points of the formula lie between the inner and outer 
sphere of the spherical shell. 

CoROLLARY. The number of evaluation points is 8(m +- 1)*. 
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Although 8(m + 1)* is the minimum number of points for the point configura- 
tion we have chosen, it is not necessarily the minimum number for other point 
configurations. For example, the formula developed by Hammer and Wymore [7] 
for 7th degree accuracy utilizes only 27 points, whereas the formula developed 
here uses 64 points for the same accuracy. The formulas presented here, however, 
give formulas of arbitrarily high accuracy. (By removing the restriction that s 
be of the form 4m + 3, which can be done, it is possible in some cases to reduce 
the number of points below the number specified by 8(m + 1)*. 

4. The weights C, and the zeros of Q,,.,:(r?)._ We list below, for m = 0 and for 
m = 1, and for four values of the inner radius R, the polynomials Q,,4,(r*), their 
zeros, and the weights C,. (Decimal numbers are truncated, not rounded.) 














m= 0 
R=0 QF") =7-3 
r2 = .60000 00000 C, = .33333 33333 
341 
R=} = Qe) = 2 - 
r? = .60892 85714 C; = .32812 50000 
93 
~ a es 
R=4 Qy=7- 75 
r? = 66428 5714 C, = .29166 66666 
2343 
~~ age <n See 
aes Ql) = 7 — 3060 
r? = .79155 40540 C, = .19270 83333 
m=1 . 
10 5 
- 7 ae 
R=0 Q.(r?) = r4 9 r+ 1 
r2 = .28994 91979 C, = .13877 79991 
re? = 82116 19131 C, = .19455 53342 
107,605 5,464,615 
= L 2) = : ’ 2 po ide tot ares te Ee 
aie Qa(r*) = 1 — “Or a7a” + 21.161,728 
ry? = .31239 3379 C, = .13890 8126 
re? = 82662 1355 C. = .18921 6873 
5905 63,005 
= 2; = o_o, ae ’ 
R=% Qz(r*) 4500" + 171,696 
ri? = 42940 5421 C, = .13053 6459 
re? = 85456 9355 Cs = 16113 0205 
2,046,079,840 782,901,015 
= 3 2 = — ’ ’ ’ 2 ’ ’ 
R=% Qo(r*) = * — 7302.513,408" + 1,302,513,408 
r? = .65958 19 C, = .09166 95 
re? = 91128 85 C. = .10103 87 
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On a Generalization of the Prime Pair Problem 


1. Introduction. One of the many unsolved problems in the theory of prime 
numbers concerns the celebrated conjecture by A. de Polignac [1] about primes 
which differ by two. Prime pairs (or twin primes) are generally believed to consti- 
tute an infinite set but to date no one has succeeded in proving or disproving this 
conjecture. 

It was established by Brun [2] that the sum of the reciprocals of all prime pairs 
is bounded but whether this boundedness is due to the finiteness of the set or to 
the ‘‘thinness” of an infinite set has not been determined. 

2. Prime n-tuples. In generalizing from prime pairs to prime n-tuples we 
must first define what is meant by a prime n-tuple. 

Definition: a prime n-tuple is a set of m odd primes :, po, --- , such that 
the difference An = p, — p; is a non-trivial minimum. 

The meaning of ‘‘non-trivial” will now be explained. For nm = 3 we can exhibit 
a set of three consecutive odd primes such that A3 = 4 (i.e., 3, 5, 7). We consider 
this a trivial minimum, however, because (with one exception) in every set of 
three consecutive odd integers one of the integers in the set is composite since 
it is divisible by 3. The single exception to this rule is the case where the number 
divisible by 3 is the number 3 itself. Generally, we consider a set of m consecutive 
primes a trivial prime n-tuple if one of the primes of the set is always a factor of 
one of the terms of any sequence of consecutive odd integers having the same value 
for An. 
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3. The generalized problem. We are now ina position to state the generalized 
prime pair problem: 


(1) For a given m, what is the value of An? 
(2) For a given n, is the set of primes which yield An finite or infinite? 


Clearly, the second question remains unanswered for all m > 2. The first 
question, however, can be answered for specific values of m by establishing the 
following conditions: 


(A) Show that for every set of n odd primes pi, po, --- , p, that An > U where 
U is some positive integer. 
(B) Exhibit a set of m odd primes such that An = U + 2. 


4. Anexample. The procedure followed in establishing condition (A) will be 
illustrated for n = 8. Beginning with 1, write down 14 consecutive odd numbers 
(mod 10) 

1387913579135 7. 


Now, assuming the first number represents a prime, the 2nd or 3rd and every 
3rd number thereafter is divisible by 3; one of the 2nd, 3rd, 4th, or 5th and every 
5th number thereafter is divisible by 5, etc., for 7, 11, 13, ---. Taking every 
possible combination of assigning factors there will be one or more combinations 
which leave a maximum number of terms of the sequence which have not been 
assigned factors. Placing the factors which have been assigned underneath the 
terms to which they have been assigned, the maximum combination for this se- 
quence is 

i3s791 
3 3 
5 7 


35791357 
a2. ot 


ee | 


We need not assign any further factors since 11, 13, 17, --- , can always be 
assigned in such a fashion as to involve only terms which have already been as- 
signed factors. If the same procedure is followed for sequences of 14 odd numbers 
beginning with 3, 7, and 9, respectively, it is found that-there is no combination 
that will leave more than 8 unassigned terms and two combinations that leave 8 
unassigned terms, i.e., 


3 91 


nao 


7 3 
3 3 


ee | 


and 


T9LISST9IVISST913 
3 53 $5 3 


If we repeat the same process for 13 consecutive odd numbers we find there is 
no way to assign the factors to leave 8 terms unassigned. Since the difference be- 
tween the first and 13th odd number is 24 we clearly have condition (A) satisfied 
when U = 24. Now, if we can find a set of 8 consecutive odd primes (mod 10) 
equal, respectively, to the 8 unassigned terms in any one of the three sequences 
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determined above we shall have proved that A8 = 26. The following primes com- 
plete the proof. 


11, 13, 17, 19, 23, 29, 31, 37. 


When conditions (A) and (B) have been established, the procedure described 
here gives all the essentially different poissible types of n-tuples for a given n 
and by solving one or more linear indeterminate equations the exact form of each 
type of n-tuple can be obtained. Thus, form = 8, there are, at most, three different 
types of 8-tuples, viz., 


(a) 210x + 11, 13, 17, 19, 23, 29, 31, 37 
(b) 210x + 173, 179, 181, 187, 191, 193, 197, 199 
(c) 30x + 17, 19, 23, 29, 31, 37, 41, 43. 


5. Some further questions. The preceding discussion raises some new questions: 

(3) Are there values of m such that An > U + 2? 

The writer has verified that An = U + 2 for all m < 10. On the other hand 
an exhaustive search of the first two billion numbers failed to disclose a 
single 19-tuple of the types prescribed by the above procedure. 

(4) Are there some possible types of n-tuples for which no n-tuple exists? 

For example, there are, at most, two types of 10-tuples; a 10-tuple for one 
type has been exhibited but not for the other; one possibility is that none exists 
for the other type. 

6. Numerical results. The writer has computed Am for n < 26 and determined 
all the possible types of m-tuples for each m in this range; these are tabulated in 
Table I. It is convenient to describe a given type of n-tuple by specifying the small- 
est prime in the sequence with respect to an appropriate modulus followed by the 
succession of differences between the constituent primes. It is observed that a 
kind of ‘“‘duality” exists between the different types of n-tuples for a given m such 
that if the set of differences for a given n-tuple is known, its dual will have 
these same differences in reverse order. Some n-tuples are self dual and do not 
yield distinct types, such as m = 6. An examination of Table I shows that for 
n = 2, 4, and 6 there is only one type of n-tuple; for » = 8, there are 3 types; 
for nm = 3, 9, 15, and 22 there are 4 types; for m = 13 there are 6 types; for all 
others in the table two types are given, each the dual of the other. In a recent at- 
tempt to discover another -tuple that is self dual the writer ascertained that for 
n = 41 there are 8 distinct types of possible 41-tuples. 

The numerical results obtained in this investigation were performed on an 
IBM 701 electronic calculator. Four separate programs were written. Program I 
was designed to search for 4-tuples; an exhaustive search of all numbers 
<5,073,379 was made and 549 4-tuples were discovered. 

Program II was designed to search for 7-tuples of the types 


(A) 210x + 11, 13, 17, 19, 23, 29, 31 
(B) 210x + 179, 181, 187, 191, 193, 197, 199. 


Seventeen 7-tuples of type (A) and twenty-four of type (B) were discovered 
in an exhaustive search of all numbers < 157,131,419. 








NH HON OH Gog tH 
HAMONANHAN WOM SSS 


>» 2 > > ba ie 





2 OS ee es ee a ee 
oe © » 8. ©- 2 2» ».6.) = > ee 
> © 22> 2 » » SB B_> - ban > i - 
- 2) @ 2 © 2". 2S. 2.2.2 2 2 > Se” 
eo @« © » » © 2 ©» >» Sy 2. 2 2 > 2S SE So oe 


> S 2 8. o 2. 2 SS SSP Sb eS he Be oe oe eS ee ee 


° 2 2 » 2 7 Sf 2 oP DP oP Sue ep SD SS > SS Se Se 


WAPCAHHA SALSA AAHSOCOCCCSCHSNS HH aa HON ON 


NAA LALA Fag ESO ENO ALO AAAAA A HOOMOOONNOSCOHTAOC SSO 

YS ELSES SNANVAEANAOAMAANAHAONONONHTOOCO HOSA ON Gat 

m PEGE SLO IO SLO SOOO. gf HNO SHO OO SHA SH SH SH SHO SHO HATO HOO HHO ONO a HO SO 
a 8 SAA LO AAAS Bf O MO AA AN HO OOOO ONO OOOO WANG HO NTSC G7 God 
<e HG FOO AAO OA OOO VO OOH GALLO ASH OSL CASTS A a Ho it Hoa teeta & 
af Me NE NS ONS HAAN GO HO MH HON ONOGOCOHOAAHTAN SHON CHAS AS 
be & AAA A AALS ALO BOAAAOAHANA HSM OHLVON COC CHO MH HA HOH HOW AO 
S SEAL AAA AS A AAA SO SHA SSH SHO ALON HOA HAN HAT ONO ONO AON ON 
ANAND APA HANAN NODPANHANNANNHNONONHTANHANNNONNHN HOGA HPONNN Hd NOdMdHadanad 
FMT RN MAROON SMM BD MAE TAMAMTATMOREONRMAONRMBARARAGTMRARMRATOR 
Sn oe be ee Bo Da SSE SBS OMOMORMAMANHVSADMMAMHPAMHAAANHAANHDNAN BH BONDAMAHHADCOE 

—_ _ — mBNONANAN CH Ht WH OC BsADAGCAGCAN nN nN i ie) ee ho | OnAw 

S330 SH HSH HH HH HH CMOMOME~M A A © © Onm S33 
= Set et et N N NO N N _ _— =H OO N o* 
SS NN mo mW MAN tN 
N oO 
©Ssssssesgcesesesosessoeooeosooooooososoooooocooeoessesooosesosooocooccocococccco 

= MOMMA MAMABAM AAMAS SSS IMM D DDD SMHS AOD DOD aa a PBR 
ANNAN AN NNN NAAR ARARAAASSSS VT °SSSSSSSSSSSSSSReeewwssos 
NNANNANANANANAAS oO oo cooocoeoo So oooococooco aa 

oF) 6) OF) OF) OF) OF) 0) 0) 0D OD et O00 St St St St et IAD 

in ih MMMM NNMwH oo 

nO“ 

NOOCOCOCBDNANLO SSCOCSSSSSONANMOONNWDDDMDDDODOOOOODSOOOSOCSOOOOOHHOSCSO ~ 

~ AA SRRSSASSSSSSSSIISSSSSRSSSSSSSSSSLLLVLessssSegsggzzesocors 
<1 ee ee ee 
N69 9 09 09 HIN IND OEE DDO AAAASCOAAAANMMMMOMMHHNNNNHNOONEMDDAKDSCOHANANANAMOdMdHMMHOO 

= tt et et et et et SS St St St St St St St HNN NNN NNN NNN AAA 








GENERALIZATION OF PRIME PAIR PROBLEM 253 


TABLE II 
Ist term of n-tuple n Type 


11 1 
1277 
5639 
88789 
1 13143 
1 13147 
1 65701 
2 84723 
6 26609 
8 55709 
10 68701 
11 46773 
25 80647 
65 60993 
75 40439 
85 73429 
119 00501 
157 60091 
178 43459 
185 04371 
190 89599 
207 37877 
210 36131 
240 01709 
256 58441 
394 31921 
429 81929 
435 34019 
450 02591 
678 16361 
691 56533 
733 73537 
742 66249 
761 70527 
792 08399 
804 27029 
841 04549 
868 18211 
879 88709 
936 25991 
1006 58627 
1240 66079 
1247 16071 
1284 69149 
1347 64997 
1362 61241 
1379 43347 
1401 17051 
1442 14319 
1546 35191 
1571 31419 
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Program III was designed to search for 8-tuples of the type 
(C) 30x + 17, 19, 23, 29, 31, 37, 41, 43. 


Ten 8-tuples of this type were obtained in a search of all numbers 
< 137,943,347. 

Program IV was the unsuccessful attempt to find a 19-tuple. As a check on the 
above calculations program II was recently re-run on an IBM 650 magnetic drum 
calculator; there was 100 percent agreement between the results of the two 
machines. 

An examination of Table I shows that all m-tuples in a given range of numbers 
for 7 < m < 15 can be determined by examining only the n-tuples obtained from 
the three types given above. This subsequent calculation yielded three 8-tuples 
and one 10-tuple for type (A), four 8-tuples and three 9-tuples for type (B) and 
three 9-tuples for type (C). These results are tabulated in Table II. Since a 9- 
tuple may also be an 8-tuple, and a 7-tuple, etc., the highest ranking designation 
is given in each case. 

7. Conclusion. The prime pair problem has been shown to be the simplest case 
of a more general problem. Essentially, the problem poses two questions: 

(5) How “clustered” can a given number of primes be? 

(6) How many such “clusters” are there? 

The first question is considered equivalent to the question, what is the smallest 
range of natural numbers that can contain primes? The scheme proposed in 
this paper gives a necessary range; it is suggested but not proved generally that 
this range is sufficient. It is proved to be both necessary and sufficient for all 
n-tuples for m < 10. 

G. H. Hardy [3] and W. W. Ball [4] have made conjectures concerning the 
number of prime pairs and higher -tuples under a given magnitude but the results 
of this investigation are too meager to expect any degree of accuracy from the 
formulas that have been proposed. 

The programs for this investigation were not the most efficient since the 
writer had a limited amount of time to write, debug and run them. With a more 
elegant approach to the programming and by using the more sophisticated logic 
of the IBM 704 these calculations could be accelerated by a factor of 200 or more. 

The writer wishes to express his thanks to the referee for his comments and 
suggestions. 
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Generation of Bessel Functions on High Speed 
Computers 


The generation of values of sets of functions such as the Bessel functions 
J,(x) and Y,(x) is a problem which one encounters frequently in numerical cal- 
culations. The present paper is concerned with describing a method which has 
been found particularly useful on high speed computers. The basic idea is due to 
J. C. P. Miller and is described in detail [as applied to the Bessel functions J,(x)] 
in the introduction to the second volume [1] of Bessel functions published by 
BAAS. The tables in the BA volume were computed by hand as were the tables by 
Fox [2]. However, the method is particularly suited to large scale computers. In 
preparing a general purpose routine to compute J,(x) and Y,(x) one is faced with 
various obstacles. The power series representation of J,(x) is useful for values of 
x which are not too large. In addition to the slowness of convergence the form of 
Y,(x) becomes complicated. The asymptotic expansions of J,(x) and Y,(x) are 
useful when x is large with respect to m. In the regiqn where both m and x are large 
neither of these representations is useful and a third form must be used. It is clear 
that any program which employs these representations and which should be 
generally useful must automatically choose the proper form for computation. The 
method of Miller avoids these difficulties and provides a simple algorithm not only 
for J,(x) but for many other functions with a similar behavior. 

Specifically, the method is as follows (the theory supporting the method is 
given in [1], p. xvii): 

1. Choose & larger than the greater of m or x, assume Ji4, = 0, J, = a ¥ 0, 
where a is some arbitrarily chosen constant. 

2. If s = 1/x, generate the sequence J,_1, Jas, --- , Ji, Jo from the recur- 
rence relation 

J 1 = 2p2J Sears J pt 
starting with » = k. If k is chosen sufficiently large we will obtain J, = cJ, for 
desired values of p ranging from p = 0 to p = n. 
3. To determine the constant c use the functional relation 


J 0 os 2>- J 2m = 1. 
m=1 

The choice of k offers no difficulty since the calculation may be done iteratively. 
Thus, starting with a value of k > max (mn, x) we perform the calculations in (2) 
and (3) once and then repeat with & increased by a fixed amount. The results ob- 
tained are compared in accordance with a preassigned tolerance and the process 
is repeated until the criteria for acceptance are satisfied. In practice, if the values 
of J, are desired for » = 0,1, --- , m the comparison need only be made at 
pb = n. It is advisable to employ floating point arithmetic and choose a as small 
as possible. 

4. To generate the functions Y,(x) we employ the following representation 
for Yo(x), 


Yo(x) = (2/z) [0 log 5 + 7} +2 5 Dees) 
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where 7 is Euler’s constant. We note that the series part here employs the same 
values Je; as occur in the calculation of the normalization constant c. 
5. Compute Y;(x) from the relation 


Ji(x) Yo(x) — Jo(x) Vi(x) = 2/2x. 
6. Generate Y,(x) in a forward direction for p = 1, 2, --- ,m — 1 from the 
recurrence relation 
Vo41 = 2pzY, = Fat 
The method may be employed in a similar manner for other sets of functions. In 
the following we outline the procedure for some of the functions for which it may 
be employed. 
Modified Bessel Functions I,(x) and K,(x). Here the recurrence relation 
(descending order) for J,(x) is 
| par = 2pzl y + I p41 
while the normalization can be obtained from the relation 
In+2> In = &. 
m=1 


While there is an analogous relation for Ko(x) corresponding to that for Yo, the 
possible loss of significant digits makes it unfeasible to use. To compute K(x) the 
integral representation 


K(x) = [ exp (— x cosh u) du 
will be convenient for all values of x > .01. (Alternatively one may use polynomial 
approximations given by Allen [3].) The relation 
Kilo + Koli = 1/x 
will yield K, while the recurrence relation (ascending order) 
Kyi = 2p2K, + Ky1 


may be used to generate as many values of K ,,; as may be desired. 

Spherical Bessel Functions j,(x) and n,(x). The spherical Bessel Functions 
jn(x) = V2/2x Jny3(x) may be generated with the aid of the recurrence relation 
(descending order) 


jo. = (2p + 1)2jp — Jos 
with the normalization constant determined from the fact that jp) = sin x/x. To 
determine the function m,(x) = (—1)"*! ~V/a/2x J_in4i)(x) we have 


sinx cosx 
and m, = — —— — = 
x x 


cos x 
u%=-> => 








and it is convenient to compute these functions directly and compute as many 
further values of , as desired with the previous recurrence relation (ascending 
order) which is also satisfied by 7,. 





ay 


eo 


me 


the 


In 
1ay 


ion 


the 


the 


nial 


ons 
‘ion 


any 
ling 





NUMERICAL INTEGRATION FORMULAS 257 


Modified Bessel Functions of Half-Integral Order. The functions 
in = Vx/2x In4x() 
ken = (—1)"*1 W/x/2x Kay3(x) 
both satisfy the recurrence relation 
(2p + 1)2ty = ipa — toys 


and may be computed in a fashion entirely similar to the computation of the 
spherical Bessel functions. We have also 


sinh x 
dailies 
ko = 
a we-*(x + 1) 
— 2x? 


IRENE A. STEGUN 


MILTON ABRAMOWITZ 
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Washington, D. C 
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Remarks on the Disposition of Points in 
Numerical Integration Formulas 


1. Introduction. Numerical integration formulas of degree 2 (i.e., exact for 
polynomials of at most degree 2) consisting of m + 1 equally weighted points have 
been developed for certain regions in n-dimensional Euclidean space. Thacher 
[1] discusses the equations which a formula of this type must satisfy for regions 
which are invariant under the group of linear transformations which leave the 
n-cube with vertices (+a, +a, --- , +a) invariant; we call these symmetric re- 
gions. Hammer and Stroud [2] give two such formulas for the m-simplex. Hammer 
[3] has shown that a set of 2” equally weighted points lying on the coordinate axes 
form a formula of degree 3 for any symmetric region. 

In section 2 of this note we show that the formulas of degree 2 discussed by 
Thacher can be described geometrically. We also show that there is a similar class 
of formulas for the regular n-simplex. In section 3 we describe geometrically a 
wide class of formulas of degree 3 containing 2n points for symmetrical regions. 
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2. Formulas of degree 2. This section is devoted to the proof of the following: 


THEOREM 1. A necessary and sufficient condition that n + 1 equally weighted 
points form a numerical integration formula of degree 2 for a symmetric region or 
for a regular n-simplex is that these points form the vertices of a regular n-simplex 
whose centroid coincides with the centroid of the region and lie on the surface of a 
sphere of radius r = V/ni 2/Io, where 


1n= | ae b= | s#d0 = tee -| X,2d0. 
R R R 


Proof. The weight for the points is J)/(m + 1). We first prove the theorem 
for symmetrical regions; the proof for the -simplex follows easily. Let R be a 
symmetrical region. We then have 





[, #00 = [ecco = 0 4,j = 1, ,n 4 xj. 

R R 

Now suppose the m + 1 points 

(1) vi = (vin, Vang + ** y Vin) 4=0,1,---,m 


are an integration formula of degree 2 for R with equal weights. Then they satisfy 
the equations 


(2) Vor + mit +--+ + mi = 0 G¢mi,-+- ie 


1)J. - 
(3) vosvog + visvig + +++ + ning = mae i; 4,j =1,---,m. 


To show that the points »; satisfy the conditions of the theorem consider the 
matrix 


é 7 








Vou Vi V21 site Val 
Vo2 Vi2 V22 iia Vn2 
Ae 
Von Vin Von oe Van 
L V I2/Io V I2/Io V I2/Io --: V I2/Io | 


By (2) and (3) it is seen that 


AAT = M@HTEUE | . 


Io 
where J is the identity matrix. Hence 


A'TA =e (n + 1)In 7 
Io 


is equivalent to the equations 


= (n + IIe, 
Io .~ 


This shows that the points lie on the sphere of radius V/nI./Io with centroid 
at the origin. To show that the »; are vertices of a regular m-simplex it suffices to 


I. ee 
(4) vavin + +++ + vinvin + F 4,j =0,1,--:,m. 
0 
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show that they are equidistant. By (4) 
@(v;, vj) = vin oes + Vin + vin ooo > Vin — Wve + --- + vind jn) 


- 2(n _ 1)I- 
= ——— 4 


Reversal of the above argument proves the condition of the theorem is sufficient. 
This establishes the theorem for symmetric regions. 

To show the theorem holds for the n-simplex consider a particular set of 
n + 1 equidistant points which lie on an n-sphere with center at the origin and 
radius r = av/1/(n + 2) . The results in [2] show that there exists a regular 
n-simplex S, for which these points are an integration formula. The vertices of 
S, lie on a sphere of radius a. The above argument shows that any ” + 1 equidis- 
tant points on the sphere of radius 7 is also an integration formula for S,. This 
completes the proof of the theorem. 


TABLE 1 
Region I2/Io 
P, a*/(n + 2) 
S. a*/n(n + 2) 
ce a?/3 
Qn 2a?/(m + 1) (n + 2) 


Values of I2/Ig are given in Table 1 for four regions (i.e., their interiors). 
P,, is the n-sphere of radius a with centroid at the origin. S, is any regular n- 
simplex whose vertices lie on P,. C, is the n-cube with vertices (+a, +a, ---, +a) 
or this cube rotated in any manner. Q, is the region defined by the 2" inequalities 
ax; + % + +--+ + x, < a, or any rotation of this region. The regions S,, C,, 
and Q, are the only regular polytopes for > 5 ((4], p. 120). In the following 
section we use the fact that any vertex of Q, is a distance of ay/2 from 2(n — 1) 
vertices and a distance of 2a from one vertex. 

3. Formulas of degree 3. In this section we prove 


THEOREM 2. A necessary and sufficient condition that 2n points v1, --- , Yn; 
—v1,°**,—0_ form an equally weighted numerical integration formula of degree 3 
for a symmetrical region is that these points form the vertices of a Q, whose centroid 
coincides with the centroid of the region and lie on an n-sphere of radius r = +/ nI2/Io. 


Proof. The weight is 5 Iy. Let the points 


(5) vi = (van, +++ 5 Vin) —vi = (—va, +++ , — Vin) += 1, »n 
satisfy 
(6) ViiVij a eee ~_ VaiPni = nls. t,j = 1,--- " Nn. 


Io 
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(The integrals of the odd degree monomials are zero and are identically satisfied.) 
Setting 





| Vil V21 abnas Var 
Ault ™ 99 es Var 
Vin Vin sis Yan 
the proof follows as in theorem 1. We find 
@(+v;, +v;) = sae xj i,j = 1, ,n 
0 
and 
4nI2 
2 ie: sta 
d (v5, Vi To 


This completes the proof. 

The results of theorems 1 and 2 generalize results given by Thacher for the 
2-cube (square) and the 3- cube. There may exist unequally weighted integration 
formulas containing fewer points than those discussed here; however, the equally 
weighted formulas will probably remain important for applications. 

4. Concluding remarks, One of the properties desirable in an integration for- 
mula is that the points be interior to the region. The 2m point formula given by 
Tyler [5] for C, using points on the coordinate axes has the points outside for 








n > 3. We give here two formulas for C, (with vertices (+1, +1, --- , +1)) of 
degrees 2 and 3 for which the points are interior for all m ([4], 245): 
Let I’, denote the point (y1, yz, --- , Yn) where 
pa 4/2 cos 2thE ai Wf sin 278 
Y2r-1 3 n wre 1 Y2r 3 n re 1 
r= 1, 2,--- , [4n} 
([4] is the greatest integer not exceeding 4m), and if m is odd y, = (—1)*/+/3. 
Then Io, Ti, --- , I, satisfy the conditions of theorem 1, and all are interior to C,. 
Let 2; denote the point (0;, a2, --- , on) where 
4/2 (2r — 1)ke if? (2r — 1)ke 
Or-1 = = cos “= Or = sin <———_—_. 
K n 3 n 
‘ ~ 9 i, 2, ary , [3n]’ 
and if  isodd o, = (—1)*/+/3. Then 21, --- , 22, satisfy the conditions of 


theorem 2 and all are interior to C,. 

For the formulas and regions we have considered, only the location of the 
points relative to each other and not to the region was found to be important. The 
insight given by these results may prove useful in further investigations; one prob- 
lem is that concerning formulas for more general regions. 

A. H. Stroup 


University of Wisconsin 
Madison, Wisconsin 


This work was supported in part by the Office of Ordnance Research, U. S. Army Contract No. 
DA-11-022-ORD-2301, and in part by the Wisconsin Alumni Research Foundation. 
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TECHNICAL NOTES AND SHORT PAPERS 
A Method for the Numerical Evaluation of Certain Infinite Integrals 


The solution of many physical problems often necessitates the numerical evalu- 
ation of infinite real integrals, a common example being that of solutions obtained 
with the aid of integral transforms. The evaluation of such integrals is often a 
laborious task, particularly if the integrand is oscillatory, so that it is usual to 
resort to special methods which give information for certain ranges of values of 
the variables; methods of this type are those involving asymptotic expansions or 
the related techniques of steepest descent and of stationary phase. The purpose of 
the present note is to outline a method in which the value of such integrals is ex- 
pressed in terms of a convergent series obtained by a modification of the corre- 
sponding asymptotic expansion. The development is given below for a special 
case only, namely one which might arise in conjunction with the use of sine trans- 
forms; it will be clear however that these results can be readily generalized to 
other types of integrals which are usually reduced to an asymptotic representa- 
tion. Examples may be found in Erdélyi [1]. The method is thus valid whether 
the integrand is oscillatory or not; in fact, though the special integrand consid- 
ered in detail below does oscillate, inspection of the convergence proofs shows that 
this fact is of little importance to the developments presented. A method which 
holds in the case of oscillatory integrands has been described by I. M. Longman 
[2]. 


Basic expansions. Consider a convergent integral J(a) of the form 
(1) I(a) = | f(x) sin xdx; f(x) — 0 steadily as x — ~. 


By f(x) — 0 steadily, we mean that f(x) > f(x.) > O if x. < x, and limf(x) = 
0; see Whittaker and Watson [3]. N successive integrations by parts may be 
shown to give the following result 


N co 
(2) I(a) = > f(a) cos [a + i(x/2)]+ / f(x) sin [x + N(x/2) dx 
i=0 a 
where f‘® = (d‘f/dx*), provided that f(x) is differentiable the required number 
of times, and that 
(2a) f(x) — 0 steadily asx — ~; 4 = 0, 1, 2, - 


The term in equation (2) containing the summation usually represents an asymp- 
totic representation of J for large values of a, and the infinite series obtained as 
N is increased indefinitely in general does not converge. A convergent expansion 
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for I(a) may now be derived in the foliowing manner. Integration by parts gives 
(3a) (a) = |" f(z) sin ade + f(a) cos ai + [ ; f(x) cos xdx 
and further 
(3b) (a) = f(x) sin xdx + ["* f(x) cos ade + f(a) 008 a1 
— f(a) sin az — eg f(x) sin xdx. 
Repetition of this process finally gives 
(4) T(a) = 3 (oss) 008 [osns + (4/2) 


he > [fox sin [x + (4/2) ]dx 


i=0 “2; 


where one may set 
(4a) Gini > Ai; ao =a. 


It will now be shown that the quantities a; may be chosen in such a manner that 
the two series on the right-hand side of equation (4) converge. 

Convergence of series expansion. The first series on the right-hand side of 
equation (4) will certainly converge if the a,’s are chosen so that the series 


@ 


(5) Si = Df (ain) 

i=0 
converges; and this series will converge (absolutely) if a positive number p in- 
dependent of 7 exists such that 


(5a) 1 > p > | fain) /f (as) | 


for all > 1. It will now be shown that such a choice of a,’s is always possible. 
(The author is indebted to Dr. C. C. Chao for his valuable suggestions concerning 
this proof.) 

Choose the quantity a; > a» arbitrarily; then the value of f (a;) is known and 
a2 must be selected so that 


(Sb) lf a2) | < p| fF (ar) | 


as may always be done because of relation (2a). Now however the value of 
f(a) is known, and so a; can be chosen by a similar procedure. Repetition 
of this process yields values of all a,’s in such a manner that relation (5a) is 
satisfied for all i > 1 and therefore series S,; converges absolutely. It should be 
noted that the choice of a,’s is not unique, and that in fact if such a choice has 
been made (a; = a,’, say) then the values a; = a,;"’ will also insure convergence 
of S; provided only that 


(6) a,’’ — a,’ 


~ 


in view of the steadiness requirement of equation (2a). 
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It will now be shown that the a,’s may be taken in conformity with require- 
ment (6) and, in addition, so that the second series of equation (4) ,namely 


(7) S. = > Ig ikea [Pe sin [x + i(x/2) ]dx 


i=0 


also converges. Note first that it follows from equation (2a) that, for any i, 
a number A ; exists such that 


(7a) | f= (x) | < | f(x) | for all x; > x > Aj. 


Let now the quantities a; be selected (consistently with inequality (6)), so that 


(7b) a; > Aj. 
Because of equation (4a) then the relation 
(7c) | fF (digs) /f (as) | < 1 


holds for all 7. 
Consider now the integrals I;; because of the steadiness requirement in equa- 
tion (2a) the quantity f‘®(x) does not change sign within a; < x < aj4: and 


(8a) |Ii|< = |f (aus) — f(a, | 








[ae 
|FO(as) 1 — CFPC + DFME)I| < 21fP@) |; i #0 


in view of relation (7c). Series S. (with the possible omission of the first term) is 
then term-by-term less than the series 


(8b) 22 | f(a) | =2 ) Lf (asus) | 


which has been shown to converge. Hence S; also converges. 
Example. As an illustration of the procedure indicated above, the special case 
of f(x) = x~* will be considered; thus 


, oo 





(9) I(a) = | x* sin xdx; k> 0. 
Here one may take (as will be shown) 
(10) a;=at+ita 
where a@ is a constant; equation (4) then reduces to 
(10a) I(a) = Si(a) + S2(a) 
where 
nie See SAS — 2) eid ~ aw el) 


‘=1 (a ta) &+i-D) 


(10b) 
2 at+ia 
S:(a) = > (1)(k)(R +1) --- (R4+2- 2) [aren cos [4(4/2) — x jdx. 
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Series S, converges if 


(k+i—1)(a+ a 
[a + (¢ + 1l)a]*® 


= lim Gti-* = ‘) lim (5 + ta ) lim (24 ie _ + to ) 
je a+ia 00 a+ (4 + 1)a i020 a+t+ (4 + l)a 
. 1 é 
“re ( ~ 1+ @/a) + ) —e 





(10c) 1> im { 





or in other words if 
(10d) a> (1/e). 


Series S; may now be considered by expanding the integrals it contains in a 
manner entirely analogous to that of equations (8a) and (8b), and it can thus be 
easily shown that this series also converges if a is chosen as specified in equation 
(10d) ; the latter condition then represents in the present case the only require- 
ment for convergence of expansion (10a). 

An advantageous choice of a, consistent with requirement (10d), isa = 2/2, 
since in this case the sine-term in S is constant; in particular, if a = 0 (ora = 
nx), note that S; = 0. No choice of a is of course possible which will make S, = 0, 
so that numerical evaluations of integrals are still necessary. Series S, converges 
quite rapidly, however; as an example consider in fact the case of a = 0, k = 
and a = 2/2 for which the value of the integral in question is well known. The 
result may be written as 

in /2 


(11) (2/2) [ (1/x) sin xdx = 1 = (2/z) ps (4 — 1)! [ xt 


(i—1) x/ 2 
cos [i(4/2) — x ]dx. 


The integrals in this summation were evaluated by Simpson’s rule with the rela- 
tively coarse interval of (1/8). The value of the summation itself may be ex- 
pressed as the limit of the sequence of the partial sums S; of the first 7 terms of the 
series; the first few terms of this sequence were found to be as follows (to four 
significant figures) : 


(ila) Si = .8727; S, = .9762; Ss =.9951; Sy= 9988; Ss = .9996 


and may therefore be said to converge fairly rapidly. Almost the same results 
were obtained when the coarser interval of (4/4) was used in Simpson's rule. 


Bruno A. BoLey 


Columbia University 
New York, New York 


This work is part of a project supported by the Office of Naval Research. 
1. A. Erp&tyt, Asymptotic Expansions, Dover Publications, Inc., New York, 1956, Chapter III. 
2. I. M. LONGMAN, ‘ ‘Note on a method for computing infinite integrals of oscillatory functions,” 
—— Phil. Soc., Proc., v. 52, 1956, p. 764-768. 
T. Waittaker & G. N. Watson, A Course of Modern Analysis, Macmillan Co., New 
York, 1948, p. 17, footnote. 
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A Formula for the Approximation of Definite Integrals of the 
Normal Distribution Function 


A simple expression has been found which may be used to yield approximate 
numerical values for definite integrals of the normal distribution function. The 
expression is as follows. 


e7x7/2 


a t2/2 - 
[, co X + Be 4x 


ox = 6. 


This approximation is satisfactory for certain applications over a broad range 
of values for the finite integration limit, as may be judged from the table below. 





© wet 

X [ exp {(—f/2)dt X -— = L a5 Difference % Difference 
0 1.253 1.250 — .003 — .2 
2 1.055 1.044 —.011 —1.0 
& .7173 . 764 — .009 —1.2 
1.0 . 398 -395 — .003 — .8 
1.2 .1675 . 1674 — .0001 0 
2.0 .0570 -0573 + .0003 + .5 
2.2 .0156 .0157 .0001 -6 
3.0 .00338 .00343 .00005 r.5 
3.5 -00058 3 .00059 2 .00000 9 4.5 
4.0 .00007 95 .00008 06 -00000 11 1.4 
5.0 .00000 0718 .00000 0730 .00000 0012 me 


Some of the properties of Mill’s ratio, the function here approximated by 
X + .8 e~“*, have been recently described by Sampford [1 ]. 


RoGER G. HART 
University of California 
Berkeley, California 


1. M. R. Samprorp, ‘‘Some inequalities on Mill’s ratio and related functions,” Annals Math. 
Stat., 1953, v. 24, p. 130. 


Some Factorizations of Numbers of the Form 2* + 1 


The author has prepared a factorization routine for use on an IBM 701 com- 
puter. In this note, we describe the routine briefly, and feport on some results ob- 
tained during the period February—April 1957 on the computer at the University 
of California, Berkeley. 

In the basic routine, an arithmetic progression is given in which divisors of a 
number WN are to be sought. Only single word divisors, that is, divisors less than 
2°, are considered, but the number N may be many words. After deleting those 
terms of the progression which are multiples of 2, 3, 5, 7, or 11, the remaining 
terms up to N}, or up to a prescribed bound, are tried as divisors of N. In order 
to delete the multiples of 2, 3, 5,.7jand 11 efficiently, use is made of the fact that 
the differences of the remaining terms of the progression repeat with a period 
$(2-3-5-7-11) = 480. The 480 required differences are computed in advance, 
and used repeatedly. 
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There are variations of the routine which make it possible to use the fact that 
+2 is a quadratic residue mod N, and which expedite handling very large num- 
bers, several related numbers, or sequences of numbers. For example, using the 
last of these, it was easy to find that the primes nearest to 2** on either side are 
2° — 31 and 2° + 53. 

The routine can theoretically factor any number N < 2”. But if nothing is 
known about the form of the factors, the maximum possible time (which increases 
roughly with N*) would be prohibitive for the larger values of N. However, any 
number N < 2** can be factored in less than a minute, and any number N < 248 
can be factored in less than an hour. 

The remainders from the various divisions are summed and the sum is re- 
corded at intervals. If a re-run is made and the remainder sums agree, this indi- 
cates that no machine errors have occurred. 

The first work undertaken was a search for factors of the early Fermat num- 
bers, supplementing work of a somewhat different character done previously 
[4], and extending the results of Selfridge [5]. It is known that any prime factor 
p of the Fermat number F,, = 2?” + 1 must satisfy the congruence p = 1 (mod 
2™+2) if m > 1. All possible divisors p < 2* of any Fermat number were tried, and 
also all possible divisors with p < 2** and p = 1 (mod 2’). The total running 
time was about five hours. No new divisors were found. The computation was 
checked by a re-run. Thus we can say definitely that there are no factors in 
this range beyond those previously found, all of which have been known for 
more than thirty years, except the two found by Selfridge [5]. In particular, 
F; has no factor p < 2; since it is composite, it must be the product of either 2 
or 3 prime factors. Also, we can say that Fi3, whose character is unknown, has no 
factor less than 2**. (After the computation described above was finished, I 
learned that some unpublished work of Selfridge, extending [5], had covered 
most of the same range. However, the work had not all been checked, and there 
were some gaps. Selfridge also tried some divisors between 2° and 2**; I was unable 
to check this part of his work.) . 

Another project was the study of the Mersenne numbers 2" — 1, where n 
itself is a Mersenne prime, that is, 7 = 2™ — 1, where m and ” are prime. Accord- 
ing to a well-known conjecture, any such number is prime. An extremely long 
computation of D. J. Wheeler (see [3], page 844) indicated that 2* — 1 (cor- 
responding to m = 13) is composite, and the conjecture therefore false. However, 
no factor of this number was found. It was decided to make a starch for factors 
p < 2** of such numbers. The factors must satisfy = 1 (mod 2m) and p = 1,7 
(mod 8). The only values of m in question are m = 13, 17, 19, 31. Two factors 
were found, namely 

1768(2"7 — 1) +1|27, 1-1 
and 
120(2® — 1) +1] 277-2 —1. 


These results can easily be verified using a desk calculator. Thus the conjecture is 
certainly false. The time required for the test was about an hour and a half. The 
computation was checked by a re-run. Thus we are sure that there are no other 
factors in the range considered. In particular, 2* — 1 has no factor less than 2". 
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Finally, factorization of various numbers of the forms 2" — 1 and 2" + 1, 
which had not previously been factored, was attempted. (A table of the factors, 
known in 1925, of numbers of these forms with  < 500 appears in Cunningham 
and Woodall [1], and a number of additions to this table were given by Lehmer 
[2]. A few other factorizations have appeared elsewhere.) A somewhat arbitrary 
selection of cases was made, consisting of the following: 

(a) 2" — 1 with m odd. The primitive factors satisfy = 1 (mod 2m) and 
p = 1,7 (mod 8). The cases tried were n = 95, 97, 101, 103, 109, 119, 121, 125, 
129, 131, 133, 137, 139, 149, 157. 

(b) 2" + 1 with m odd. The primitive factors satisfy = 1 (mod 2m) and 
p = 1, 3 (mod 8). The cases tried were m = 71, 101, 103, 107, 109, 113, 115. 

(c) 2" + 1 with m = 0 (mod 4). The primitive factors satisfy p = 1 (mod 4m). 
The cases tried were nm = 104, 112, 116, 124. (For the cases m = 128, 256, etc., see 
the discussion of Fermat numbers above.) 

(d) 2" + 1 with m = 2 (mod 4). The primitive factors satisfy p = 1 (mod 2). 
The identity 


Dsit2 a i= (2261 — 2! os 1) (2241 _ Qt oa 1) 


shows that in this case the factorization of 2" + 1 is reduced to factorization of 
numbers of the form 2?**! + 2! + 1. Four numbers of this type were tried, 
namely 267 — 2% + 1, 267 + 2% + 1, 27 — 2% + 1, and 2% — 2%? + 1, of which 
the first two are the factors of 2 + 1, and the others are factors of 2 + 1 and 
qe 4 4. 

In each of the thirty cases, all possible divisors less than 2*° were tried. The 
time for each run was about half an hour, except that when factors were found, 
this was reduced, in some cases very sharply. The following factorizations were 
found: 


287 — 2% + 1 = 5-269-42 875 177-2 559 066 073, 
287 + 2% + 1 = 15 152 453-9 739 278 030 221, 
27 + 1 = 3-56 409 643-13 952 598 148 481, 


2% — 1 = 31-191-524 287-420 778 751-30 327 152 671, 

21% — 1 = 745 988 807-X, 

21% + 1 = 3-104 124 649-Y, 

2'% + 1 = 449-2689-65 537-183 076 097-358 429 848 460 993, 
21% + 1 = 3-227-48 817-636 190 001-491 003 369 344 660 409, 
2157 — 1 = 852 133 201-Z. 


All the factors written out are less than 2°, and hence are shown to be prime by 
the routine. This fact was checked by an additional run. The three factors X, Y, 
Z exceed 2°, and are of an unknown character. 

Thus we have completely factored the numbers 2" + 1, 2% — 1, 2" + 1, 
2"3 + 1, and 2™ + 1, the last factorization being obtained by multiplying to- 
gether the first two of the above equations. (The entry for 2“ + 1 in Cunning- 
ham and Woodall [1] erroneously gives the product of the two largest prime fac- 
tors of 267 — 2% + 1 asa prime.) We have also found factors of two Mersenne 
numbers, 2! — 1 and 21%” — 1, for which no factor was previously known, and a 
factor of (2? + 1) /3. 
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No factors other than those shown above were found in any of the cases tried, 
except for the algebraic factors and the factors less than 300,000 which appear 
in Cunningham and Woodall [1]. Check runs were not made, but it appears to 
be quite unlikely that any factor less than 2*° was missed. 

Appendix. At the suggestion of the referee, two lists are included which show the 
progress which has been made in factoring numbers of the form 2" + 1. These 
lists have been prepared with the help of D. H. Lehmer and J. L. Selfridge. First 
we have a list of all the cases we could find in which complete factorizations have 
been claimed, where for 2" — 1 only odd values of m are considered. 


2” — 1: = 1-99, 105, 107, 111, 113, 115, 117, 123, 127, 129, 135, 151, 521, 607, 
1279, 2203, 2281 

2" + 1: = 0-102, 105, 106, 108, 110, 111, 112, 113, 114, 118, 120, 122, 123, 126, 
130, 134, 135, 138, 146*, 148, 150, 154, 162, 170, 174, 182, 186*, 190*, 
198, 210, 234*, 270 


Some of these were not known to the author at the time the work described above 
was carried out. It should be mentioned that in several cases there is doubt that 
the factorizations are in fact complete; this is true, in particular, in the cases 
marked with an asterisk. Notice that M. Kraitchik [6] had already given a sup- 
posedly complete factorization of 2° — 1, but that we found above that a further 
decomposition of one of his factors is possible. 

The second list concerns the ‘‘original Mersenne numbers’’ (that is, numbers 
of the form 2? — 1 where is prime and p < 257), and brings up to date a similar 
list by R. C. Archibald in MTAC [7]. 


p Character of 2? — 1 
2,3,5,7,13,17,19,31,61,89,107,127 Prime 
(All other  < 100) ,113,151 Composite and completely factored 
163,173,179,181,223,233,239,251 Two or more prime factors known 
109,131,157 ,167,191,197,211,229 Only one prime factor known 


101,103,137 ,139,149,193,199,227,241,257 Composite but no factor known 


RAPHAEL M. ROBINSON 
University of California 
Berkeley, California 


The operation of this computer is supported in part by the National Science Foundation. 
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On the Solution of “Jury” Problems with Many Degrees of Freedom 


1. Introduction. In a recent numerical investigation using the Differential 
Analyser, it was found necessary to solve differential equations of up to the eighth 
order with two-point boundary conditions, the so-called “Jury’’ problem. Now, 
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for second order equations with only one degree of freedom various processes of 
inverse interpolation may be used to estimate, from trial solutions, the correct 
initial conditions for the required solution. Of these, perhaps the best for this 
purpose is that due to Aitken [1, 2], which is, as it stands, equally suitable for 
direct or inverse interpolation. An example of its application is given in section 2. 

The method given below is an extension to many variables of a similar process. 


If we have m functions f(u, v, --- , 2), g(u,v, --- , 2), --- , m(u,v, --- , 2), of 
the same m independent variables, it is necessary to interpolate between n + 1 
sets of values of f, g, --- , m, to find the values of u, v, --- , z, required to form 


the wanted values of the functions. With more than one variable we are restricted 
to linear interpolates only, as the extension to higher order interpolates turns 
out to be unprofitable and difficult to generalise. It is found to be simpler to use 
the process of linear inverse interpolation iteratively, obtaining from the first 
n + 1 sets of values a further set based on the first linear interpolation, then 
using this and the best » preceding sets to repeat the process. 
2. Aitken’s process for one variable. This may be illustrated by means of an 
example of a “Jury” problem. Consider the differential equation 
d’x 
za 


where the solution is required satisfying the boundary conditions 


dx 2. = 
+7 +* = 0, 


x(u) = a,x(v) = b. 
We may then integrate forward from ¢ = u with two arbitrary initial values of 
x’, say x1’ and x’, to reach ¢ = v, when x = c and @ respectively. Defining 
d = c — b, we may interpolate d = 0 to find, in the usual notation, 

x3’ = x’(0; di, de) 

as a starting value of x’ for a new run. This produces d;, and further interpola- 
tion yields both x'(0; di, ds) and x’(0; d), do, ds). The sequence of operations may 
be shown as: 
X1 C1 a 
X2 Ce dz x’ (0; dy, dz) = x3" 
x3" C3 ds x’ (0; di, ds) x’ (0; di, de, ds) = x4’ 
x4! C4 d, x’ (0; di, ds) x’(0; di, de, ds) x’ (0; di, de, ds, dg) = x5! 


and so on. 


It should be noted in passing that convergence is not guaranteed, since a solution 
of the type we want may not exist. Results for the above equation with the bound- 
ary conditions x(0) = 1, x(1) = 0 are shown in Table 1. 


TABLE 1 
run x’ initially x finally parts 
1 —1 +0.1875 +1875 
2 —1.2 +0.0811 —1.3524 + 811 
3 —1.3524 —0.0015 —1.3496 —1.3497 — 15 
4 —1.3497 —0.0003 —1.3491 —1.3491 —1.3490 -— 3 
5 —1.3490 0 — — _ 
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3. Extension to many variables. Let us consider the case of three variables 
(three degrees of freedom) for simplicity in writing. The extension to the general 
case from this is immediate and obvious. 

Then we have three parameters u, v, w, whose values at the start of the range 
have to be estimated so as to make three functions f, g, h of these parameters 
zero at the end of the range. 

Nowf = f(u,v, w), g = g(u,v, w), h = h(u, v, w) may be inverted to give 
u = u(f,g,h),v = v(f, g,h), w = w(f, g, h), provided that the Jacobian of the 
transformation does not vanish at any of the points in question. This we may 
assume for any given set of numerical functions. 
ou ow 
on a 
U, V, W, U;, --- , Wi, respectively, then we obtain from Taylor’s series truncated 
after its first-order terms 


Denoting the values of u, v, w, taken when f = g = h = Oby 


u= U+fU; + gU, + hU i. 
From four such equations (derived from the approximations 1, te, %3, %4) 
fiU; + 2U, +mUn+ U = mM, 
f2U; + g2U, + mU, + ;U = mw, 
fsUs + gsU, + hsU, + U = us, 
faUz3 + gsUg + 4Un + U = Uy, 


we may eliminate U;, U,, U; to obtain the required interpolate U. This set of 
equations may be solved by the Gauss elimination method [3] to give U directly, 
with no back-substitution. 

Further, the sets of equations for V and W differ from the above set for U only 
in that ~; to uw, are replaced by v7 to v4 and w, to w, respectively. Hence the whole 
process may be represented in matrix form as 





fi: gs om 1) [ Uy Vy W, 7] [m1 mh wm | 
fe g2 Ie 1 U, V, WwW, Uz Ve We 
fa gs hs 1 Ur Va Ws 7 te ts ws | 
fe ge he 12 Le Pony | Ue Ve We | 

















and the elimination process may be carried through to give a result of the form 
[FU FV FWl)=[a 6 @w], 


whence U = i/F, V = i/F, W = @/F. This process may then be repeated with 
four sets of values comprising (U, V, W) and the three best approximations from 
(uy, 01, Wi), *** , (Ma, V4, Wa). There is a corresponding reduction in labour if the 
asymmetrical Cholesky process [3] is used. 
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Formal solution of the matrices gives 
_ | 4 fe gs he | 
|1fegsha|’ 


which is seen to be a (formal) extension of the linear interpolation formula in the 
Aitken process, since the usual formula 


flusfi, fe) = (ufo — uefi)/(fe — fi) 


can be rewritten 

| m1 fe | 

f(u; fu fe) = 

| 1 fe | 

4. Example. As an example of this process, consider the steady-state solution 
of van der Pol’s equation subject to a forcing function, 

d*x 
dt? 








dx 
— — 2 — = 
(1 x") Zi + x = cost. 


If x(0) = a, x’(0) = 5, we may be required to estimate both a and 6b so that 
x(2x) = a, x’(2x) = b, that is, the solution has the same period as the forcing 
function. For a given (a, b) let x(24r) = c, x’(24r) = d, and define f = a — ¢, 

= b — d. Then we may integrate forward with three independent sets of (a, b) 
and then interpolate inversely to estimate values 4, 6 corresponding to f = g =0. 
The results are shown in Table 2, and exhibit quite a rapid convergence. 


TABLE 2 
Computed 
Zz using 
run a b c d f 2g a b runs 
+1 0 +1.227 —1.411 -—0.227 +41.4i1 — _ 
+1 —1 +1.140 -—1.499 -0.140 +0.499 — _ 


+1.200 —1.600 +1.196 -—1.443 +0.004 -—0.157 +1.2271 —1.4106 1,2,3 
+1.227 -—1.411 +1.197 —1.441 +0.030 +0.030 +1.1903 -—1.4482 2,3,4 
+1.190 —1.448 +1.190 —1.448 ~— — —_ aa <— 


Or whe 


5. Conclusion. The method described above has been used on various differen- 
tial equations and has shown quite good convergence. When a series of solutions 
with varying parameters have to be found, it is usually possible to ‘“‘predict”’ a 
set of values lying near the required ones. The convergence is then more rapid. 

This note is published by permission of the Director of the National Physical 
Laboratory. 

F. J. WARNER 


National Physical Laboratory 
Teddington, Middlesex 
England 


5 i. L. M. Mitne-TuHomson, The Calculus of Finite Differences, Macmillan, London, 1933, 
3.81. 

2. M. C. K. TWEEDIE, “A modification of the Aitken-Neville Linear iterative procedures for 
polynomial interpolation,” MTAC, v. 8, 1954, p. 13-16. 

3. L. Fox, “Practical solution of linear equations and inversion of matrices,” NBS Applied 
Mathematics Series, v. 39, 1954, p. 1-54. 
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106(C, D].—Juan Garcia, “Nuevas Tablas de Logaritmos,” Las Ciencias, 
Madrid, Spain, v. 19, 1954, p. 567-592. 


A new, five-place “triple entry’’ table of 5D logarithms of numbers and of 
circular functions is described with the aid of representative excerpts from the 
table. Instructions for its use are in the form of varied and illustrative examples. 

With regard to the logarithmic tabulation of the angles in the intervals 
(0°, 3°) and (87°, 90°), the table provides a direct and simple interpolation scheme 
for such angles; the arrangement is “triple entry.” 

Uniformity in the methods of tabulation and interpolation for the angles 
throughout the interval (0°, 90°) is attained. The table is compact, and its entries 
can be read simply and efficiently. 


GENEVEVO L6PEz 
University of California 
Los Angeles, California 


107[F].—D. N. Leumer, Factor Table for the First Ten Millions Containing 
the Smallest Factor of Every Number not Divisible by 2,3, 5, or 7 between the 
Limits 0 and 10017000, Hafner Publishing Co., New York, 1956, 9 p. in double 
col. + 476 p., 43 cm. X 30 cm., oblong. Price $22.50. Originally issued as 
Carnegie Institution of Washington Publication No. 105, 1909. 


[F].—D. N. Lenmer, List of Prime Numbers from 1 to 10006721, Hafner 
Publishing Co., New York, 1956, 8 p. in double cols. + 133 p., 43 cm X 30 
cm., oblong. Price $15.00. Originally issued as Carnegie Institution of Wash- 
ington Publication No. 165, 1914. 5 


The availability of these two monumental tables is a cause for rejoicing among 
all devotees of the theory of numbers—and only a very few of us are not. The 
author begins his preparatory note to the first of these volumes, ‘‘Factors,’’ with 
the sentence: “The value of a factor table depends chiefly on its freedom from 
errors.”’ Time has certainly proved the value of this book for no errors have been 
reported in it; only two serious errors have been found in the companion volume, 
“Primes.” They are listed here. 


p. 11, col. 13, line 1, for 8151 read 8051. 
p. 14, col. 30, line 55, for 51 read 47. : 
p- 99, col. 20, heading, for 224 read 724. 
p- 119, col. 25, heading, for 83 read 883. 


In the review copy the first of these errors has been corrected while the second 
has been marked but not corrected. The first three were noted by D. H. Lehmer 
[1]. The last was pointed out by E. G. H. Comfort (Ripon College, Wisconsin). 

Professor Comfort has also pointed out that the description of the Kulik 
tables in the ‘‘Primes’”’ need to be revised in the light of the note by S. A. Joffe 
[2]. He also noted that he received a copy which had pages 68 and 69 replaced by 
28 and 29 and pages 72 and 73 replaced by 48 and 49. The review copy appears 
faultless in this respect. 
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Both books have been reproduced by photographic methods from copies of the 
original edition—which were themselves reproduced photographically from type- 
script. The original editions were some of the earliest major tables to be produced 
in this manner. 

For those who are not familiar with the tables we give brief descriptions. 
“Factors” gives the least prime factor of each » < 10017000, provided n is not 
divisible by 2, 3, 5, 7. The most convenient way of using it is to find g, r such that 

= 210g + 7r,0 <r < 210 (using a desk calculator, for instance). Then, turning 
to the page given by the first three digits of g and the line given by the last two 
digits, we read in the column with heading r, the least prime factor of m, or a dash 
indicating that m is prime. If r is not among the column headings, a factor 2, 3, 5, 
or 7 is present and must be removed before using the table. Repeated entry into 
the table gives the complete factorization of ». By means of this a much greater 
range can be covered than the other modern tables BAAS [3] which give complete 
factorization. 

In “Primes” we find a listing, 5000 to a page, of all the primes up to 10006721. 

Both volumes contain careful introductions which describe the construction 
of the tables and their checking and printing. In ‘‘Primes’’ there are tables giving 
comparisons of r(x), the number of primes <x, with various approximate for- 
mulae. In ‘‘Factors’’ there is a list of errors discovered in earlier tables. 


Joun Topp 


California Institute of Technology 
Pasadena, California 


05, D. H. LesmeEr, Guide to Tables in the Theory of Numbers, Bull. National Research Council 
105, 1941, p. 161. 

S. A. JorFe, Editorial Note in connection with KuLix’s Magnus Canon Divisorum . 
um 48, MTAC, v. 2, 1946, p. 139-140. 

3. J. Peters, A. Lonse, E. J. TeErNoutn, & E. Girrorp, Factor Table giving the Complete 
Decomposition of all Numbers less than 100,000. (British Association for the Advancement of 
Science, Mathematical Tables, v. 5.) London, BAAS, 1935. 

4. J. KavAn, Factor Tables giving the Complete Decomposition into Prime Factors of All Numbers 
up to 256,000 ..., Macmillan, London, 1937. [RMT 196, MTAC, v. 1, 1945, p. 420-421.] 


108[F, L].—D. H. Leumer, “Extended computation of the Riemann zeta- 
function,’’ Mathematika, v. 3, 1956, p. 102-108. 


This paper describes the methods used by the author in the computation of the 
first 25,000 zeros of the Riemann zeta-function {(s) which was carried out on 
SWAC. All these zeros have Res = %. 

A table of the coefficients of 4 polynomials entering in the calculation is in- 
cluded, and there is a table of the number of failures of Gram’s law divided into 
3 types and tabulated in 15 sets of 1000 roots each, beginning with root number 
9892, which overlaps the previous run [1 ]. 


Morris NEWMAN 


National Bureau of Standards 
Washington, D. C 


D. H. Lramer, “On the roots of the Riemann zeta-function,’’ Acta Mathematica, v. 95, 
1956 p. 291-298 [MTAC, Review 52, v. 11, 1957, p. 107-8]. 
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109[F].—A.Bert GLopEN, Table de factorisation des nombres N* + 1 dans 
Pintervalle 3000 < N < 6000, published by the author, rue Jean Jaurés, 11, 
Luxembourg, 25 p., 30 cm., mimeographed. Price 125 francs belges. A copy 
deposited in the UMT Fre. 


This table is a revision of UMT 108 [1] with all unknown factors > 8.10°. 
The author intends to extend the tables [2] of solutions of x* + 1 = 0 (mod p) 
to 10°. The UMT files also contain an earlier revision of UMT 108 published in 
1952. 


JoHN SELFRIDGE 
University of California 
Los Angeles, California 


1. A. GLopEN, Table de factorisation des nombres N* + 1 dans l’intervalle 3001-6000. [UMT 
108, MTAC, v. 4, 1950, p. 224.] 

2. S. Hoprenot, Table des Solutions de la Congruence xt = —1 (mod N) pour 100000 < N < 
200000, Brussels, Librairie du ‘‘Sphinx,’’ 1935. [RMT 48, MTAC, v. 1, 1943, p. 6.] 

3. A. GLopEn, ‘Table des solutions de la congruence X‘ + 1 =0 (mod p) pour 2.108 < p< 
3.105,” Mathematica (Rumania), v. 21, 1945. [RMT 280, MTAC, v. 2, 1946, p. 71-72.] 

4. ALBERT DELFIELD, ‘‘Table des solutions de la congruence X‘ + 1 = 0 (mod p) pour 300000 
< p < 350000,” Institut Grand-ducal Luxembourg, Section des Sciences, Archives, v. 16, 1946. 
[RMT 346, MTAC, v. 2, 1947, p. 210-211.] 

5. ALBERT GLODEN, Table des solutions de la congruence x* + 1 =0 (mod p) pour 350.000 < 
p < 500.000, Luxembourg, author, 11 rue Jean Jaurés, and Paris, Centre de Documentation Uni- 
versitaire, 1946. [RMT 410, MTAC, v. 2, 1947, p. 300-301.] 

6. ALBERT GLODEN, Table des Solutions de la Congruence x* + 1 =0 (mod p) pour 500 000 > 
p< . 000, Luxembourg, author, rue Jean Jaurés 11, 1947. [RMT 491, MTAC, v. 3, 1948, 
p. 96. 


7. A. GLODEN, Solutions of x* + 1 = 0 (mod p) for 600000 < p < 800000 [RMT 1169, 
MTAC, v. 8, 1954, p. 77.] 


110[F].—Joun Leecu, ‘‘Table of groups of 4, 5, 6, 7 primes from 50 to 
100 17000,” 1 sheet (photostat), 34 X 23 cm., deposited in UMT Fie. 


The table gives 15” where four primes between 50 and 100 17000 are generated 
by adding and subtracting 2 and 4 from 15”. Those values for which additional 
primes are generated by addition and subtraction of 8 and of these the values 
which give additional primes by addition and subtraction of 14 and then of 16 
are suitably marked. Double fours, primes of the form 210” + (11, 13,17, 19) are 
also marked. 

From author’s remarks 


111[F].—R. M. Rosinson, ‘Table of factors of numbers one unit larger than 
small multiples of powers of two,’’ Los Angeles, 1957, 312 pages listed on an 
IBM 402 from SWAC punched card output. The table and the punched cards 
have been deposited in the UMT FILE. 


The main table tells whether any number of the form N = k-2" + 1 is prime 
or composite for odd k < 100 and all m < 512. If the least factor of N is less than 
10‘ it is listed. For k = 1, 3, 5, 7 the bounds on # are larger, namely 2272, 1280, 
1536, 1280, and for » < 2272 least factors to 10° are listed. Each prime found was 
tested as a possible factor of any Fermat number F,, = 2?” + 1. 

In addition to listing the 14 new factors of Fermat numbers previously re- 
ported by Robinson [1], the table lists 95-2 + 1 as a factor of Fs. This brings 
to 30 the number of known composite F,,. There are 118 primes greater than 278° 
in the table, but only the three already mentioned [1] are greater than 25°. 
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There is a list of the 42 Cullen numbers n-2* + 1 with n < 10* for which 
Cunningham [2] did not find a small factor. These were all tested, and only the 
smallest, 141-2"! + 1, is a prime. 

Robinson has presented this wealth of information in a very careful, pleasing, 
and useful format which is described in the four-page introduction. This was 
achieved despite the poor type font and limitations of the IBM 402. 


J. L. SELFRIDGE 
University of California 
Los Angeles, California 


1. R. M. Rosrinson, “Factors of Fermat numbers,”” MTAC, v. 11, 1957, p. 21-22. 


2. A. CunnincHAM & H. J. Woopa.t, “Factorisation of Q = (2¢ g) and (g.2¢ + 1),” 
Messenger Math., v. 47, 1917, p. 1-38. 


112[F ].—Emma Leumer, “On the location of Gauss sums,” MTAC, 1956, v. 


10, p. 194-202. at 
The generalized Gauss sum of order k is defined by Si = >> exp (2xim*/p), 
pb =kf + 1 and prime. ated 


It is known that —(k — 1)\/p < Si < (k — 1)\/d. The author investigates 
the distribution of S; for k = 3, 4, 5, 7 with respect to various partitions of the 
interval (—(k — 1)\/p, (k — 1)\/p). The research was undertaken to test 
Kummer’s conjecture which states that 5S; falls in the intervals (—2+/p, 
—V/?), (—V/b, Vd), (Vd, 20/p) with frequencies of 1 to 2 to 3. It was found 
that for the first 1000 primes of the form 6m + 1, these ratios became 3 to 4 to 5, 
tending to disprove the conjecture. 


Morris NEWMAN 
National Bureau of Standards 
Washington, D. C 


113[F].—R. J. Porter, Irregular Negative Determinants of Exponent 3n with 
their Critical Classes. Part II, from —D = 50,000 to 100,000. 203 typewritten 
pages, 25.5 X 10 cm., deposited in UMT Fre. 


This table is a supplement to Table of Irregular Negative Determinants of Ex- 
ponent 3n up to —D = 50,000, which were previously deposited in the UMT 
Fitz [MTAC, Review 3, v. 9, 1955, p. 26 and MTAC, Review 8&4, v. 9, 1955, 
Part I, p. 198]. 


R. J. PorTER 
266 Pickering Road 
Hull, England 


114[G, X].—D. H. Lesmer, “On certain character matrices, “Pacific J. Math., 

v. 6, 1956, p. 491-499. 

Two classes of matrices M, for which det M, \;(M), (M-"),; (if existing) 
and (M*);; can be given explicitly, are described. (Here \;(M) denotes the 
characteristic roots of M.) The matrices are of order p — 1, where # is an odd 
prime, and they are defined in terms of the Legendre symbol 


. 0 if p divides n 
x(n) = (*)- —1 if the congruence x*? = m (mod ) is impossible. 
4 +1 otherwise 
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The general element of a matrix in the first class is 


Mi; = a + bx(t) + cx(j) + dx(i) 
for any a, b, c, d. This is singular except for p = 3. The general element of a matrix 
in the second class is 
My =c+xlat+i+)) 


where c is any constant and a any integer. 
These matrices are, incidentally, useful for the evaluation of programs for 
inverting and finding the characteristic roots of matrices on computers. 


Joun Topp 


California Institute of Technology 
Pasadena, California 


115[L].—S. L. BrLousov, Tablitsy normirovannykh prisoedinennykh polinomov 
Lezhandra (Tables of normalized associated Legendre polynomials). Akad. 
Nauk SSSR, Energeticheskii Institut im G. M. Krzhizhanovskogo (Academy 
of Sciences of the USSR, Krzhizhanovskii Institute of Power). Moscow, Press 
of the Academy of Sciences, 1956. 380 p., 27 cm. Price 23.70 rubles. 


These are tables of the associated functions P,™ (cos @) defined by 


P,™ (cos 0) = /*F : “ + ar P,™ (cos @), 


the normalization resulting from this definition being such that 








ie [Pa™(x)Pdx = 1. © 


The functions tabulated are thus identical in every way with those whose behavior 
is exhibited graphically up to m = 8 in four diagrams in the second and subsequent 
editions of Jahnke and Emde [1 ]. Except for the difference of symbolism, they 
are also those denoted by F,,™ (cos @) in Section 16 of FMR Index. 

The arguments in the Belousov tables are m = 0(1)36, » = m(1)56, @ = 
0(2° .5)90°. Each pair of values of m, n belongs to one column of 37 entries corre- 
sponding to the 37 values of 6. The function values are to 6D, without dif- 
ferences. Nine figures were kept in the computations. . 

The author refers not only to several well-known tables, but also to tables by 
fa. M. Kheifets [2], which are stated to give values of associated functions for 

= 1(1)12, = m(1) 20, and each 5° of 6. On page 15 are given nine corrections 
(two of them taken from Kheifets) to numerical values contained in R. and L. 
Egersdorfer [3 ]. 

A. F. 

1. E. JAHNKE & F. Empe, Funktionentafeln mit Formeln und Kurven, second edition, Leipzig 

and Berlin, 1933, p. 178-179; also with various page numbers in later editions. 


2. fA. M. Kueirets, Tablitsy normirovannykh prisoeuinennykh polinomov Lezhandra (Tables 
of normalized associated Legendre polynomials), Moscow, Gidrometeoizdat, 1950. 

3. R. and L. EGERsp6ORFER, Formeln und Tabellen der zugeordneten Kugelfunktionen 1. Art 
von n = 1 bis n = 20. I. Teil: Formeln. Reichamt fiir Wetterdienst, Wéissenschaftliche 
Abhandlungen, Bd. !, Nt. 6, Berlin, 1936. 
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116[L].—Grorce C. CrarK & Stuart W. Cuurcu1, Tables of Legendre 
Polynomials P,, (cos @) for n = 0(1)80 and @ = 0°(1°)180°. Engineering 
Research Institute Publications, University of Michigan Press, Ann Arbor, 
Mich., 1957, ix + 92 p., 28 cm. Price $4.50. 


This is a straightforward table of the ordinary Legendre polynomials P,,(cos @). 
As the values Po (cos 8) = 1 and P, (cos0) = 1 are not listed, the tabular argu- 
ments are m = 1(1)80 and @ = 1°(1°)180°, the limit 180° rather than the usual 
90° having been chosen for convenience in application. The function values are 
to 6D, and are believed accurate to the last place given, since they were calculated 
to 11D on the Michigan Digital Automatic Computer (MIDAC). No differences 
are provided. Cosines required in the computations were taken from the NBS 
15-place tables. 

The table, whose special virtue is the high upper limit of 2, was computed in 
connection with a research program on the engineering applications of light- 
scattering which has been in progress in the Department of Chemical and Metal- 
lurgical Engineering of the University of Michigan for the past ten years. 


A. F. 
117[L].—Haro_p K. Crowper & GEorGE C. Francis, Tables of Spherical 
Bessel Functions and Ordinary Bessel Functions of Order Half an Odd Integer 
of the First and Second Kinds. Ballistic Research Laboratories, Memorandum 


Report No. 1027, Aberdeen Proving Ground, Maryland, 1956, 86 p., multi- 
lithed, 83’ X 11’’. Two copies deposited in UMT Fie. 


This gives tables of the functions 


inte) = 4/% Jule), male) = 4/Z Yous 


Ina4(X), Vn+4(x) 


for x = 1(1)50, 2 = 0(1)N, in all four cases, where N is such that 


and 


Vniy(x) < 10% < Vuys(x). 
Some typical values of N are 


x 1 10 20 30 40 50 
N 10 30 45 59 72 84. 


Values are given to 9 decimals for » < x and to 7 significant figures for m > x. 
These tables only partially overlap earlier tables from which the functions 
can be derived. We mention a few, extracted from material to appear in the pro- 
jected second edition of the FMR Index of Mathematical Tables [1]. 
Bessel Functions Jns4(x) 


12 decimals, m = w + 3(1)M + 4, x = 1(1)20, 
with | J,-4(x) | > 1, | Jnsa(x) | < 3-10-" 
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for each x, in the BAAS Report for 1925, page 221; this table was prepared by J. 
R. Airey [2]. 


6 decimals, 2 —64(1)344, x = 1(1)20 


—63}(1) + 6}, x = 21(1)30. 


nN 


This table is by E. C. J. von Lommel [3], and was partially reproduced (up to 
= 18}) in G. N. Watson’s Treatise on Bessel Functions [4]. It was similarly 
reproduced, to 4 decimals only, in the various editions of Funktionentafeln mit 
Formeln und Kurven by E. Jahnke and F. Emde [5 ]. 
A more recent tabulation by K. Reitz [6] gives 5 decimals or figures for 


n = —303(1) + 313, x = 0(0-2)20. 


Stokes’s Functions 


jn(x) 7 fx. Jura) n(x) = fz Ynia(x) a (—)*"*4y_,-1(x). 


The two volumes produced by the National Bureau of Standards Computation 
Laboratory [7] give 8- to 10-figure values, with central second or second and 
fourth differences for interpolation for m = —22(1) + 21 with 


= 0(.01)10(.1)25, 


in general, and for m = —31(1) — 23, 22(1)30 with x = 10(.1)25. The auxiliary 
function A,(x) = (4x)"m!J(x) is also tabulated for 
29 1 61 


61 
=S=—lU a _— — 1 —_— = ‘ 
m D (1) a? 9! ) a? J 0(.1)25 
with from 7 figures to 9 decimals, with central x-wise differences. 
Riccati Bessel Functions 


Xjn(x),  (—)*xyn(x) = —xj—-n-1(x). 


Here we mention only the book of tables by R. O. Gumprecht and C. M. 
Sliepcevich [8]. This gives values of the first function to 6 decimals and of the 
second to 6 or 7 figures, in each case with the first derivatives, for 


= 0(1)n + d, x = 1(1)6(2)10(5)100(10) 200(50) 400 


where d increases fairly steadily from 4 at x = 1 to 20 atx = 400. There is also 
another inconvenient table giving xj,(x) and its derivative only for very miscel- 
laneous arguments x from 1-2 to 640, and for m = 0(1)N, where N varies errat- 
ically from 6 near x = 2 to 436 at m = 640, usually §x < N < x + 10 approx- 
imately. 

J. C. P. MiL_er 
University Mathematical Laboratory 


Corn Exchange Street 
Cambridge, England 


1. A. FLetcHER, J. C. P. Miter, & L. RosENHEAD, An Index of Mathematical Tables, 
—— ¢ Computing Service Ltd., London, 1946. 
BAAS, Mathematical Tables ‘Report, “Bessel functions of half-odd integral order,” Cambridge 
University Press, 1925, p. 221-233. 
3. C. J. von LoMMEL, Munich Abhandlungen, v. 15, 1886, p. 529-664. 
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4. G. N. Watson, A Treatise on the Theory of Bessel Functions, Cambridge University Press, 
New York, 1922. 


5. EUGENE JAHNKE & Fritz Empe, Funktionentafeln mit Formeln und Kurven, Teubner, 

Leipzig and Berlin, 1909 and 1948. 
Reitz, Institut fur Praktische Mathematik, Darmstadt, Reports, ‘Tabellierung 

Besselscher Funktionen,” 4. Bericht, p. 2-9; 6. Bericht, p. 2-31; 9 Bericht, p. 2-37, 1945. 

7. NATIONAL BUREAU OF STANDARDS COMPUTATION LABorATORY, Tables of Spherical Bessel 
Functions, v. | and II, Columbia University Press, New York, 1947. 

8. R. O. Gumprecut & C. M. Siiepcevicu, Tables of Riccati Bessel Functions for Large 
Arguments and Orders, Engineering Research Institute, University of Michigan, Ann Arbor, 1951. 


118(L ].—F Rank E. Harris, “Tables of the exponential integral Ei(x),’”’ MTAC, 

v. 11, 1957, p. 9-16. 

This paper gives basic values of high precision for use in molecular structure 
calculations, with special reference to the region 4 < | x | < 50 in which | x | is 
inconveniently large for the Taylor series and too small for the asymptotic formula. 

Table 1. Values of —Ei(—x) to 18S and of —e*Ei(—x) to 19D for x = 
1(1)4(.4)8(1)50. 

Table 2. Values of Ei(x) to 18S and of e-*Ei(x) to 19D for the same values of 
x as in Table 1. 

Table 3. Tables of the functions 


h2 h" = 
R, (h) = nif -(1+4+E4 ve +e] 
2! n!\ 


which occur in a special interpolation formula. Values of R,(h) are listed for 
h = .01(.01).05, .1(.1).5, 1 to as many decimal places and up to as large values 
of m as are necessary. 


Table 4. Constants, namely Euler’s constant y and the modulus log we to 
24D, and e+" to 25D for the same values of / as in Table 3. 


A. F. 


119[L, S].—R. B. Drnctz, D. Arnot, & S. K. Roy, “The integrals A,(x) = 
(p!)— fe e(e + x) ede and B,(x) = (pl!) JP e(e + x)—e~de and their 
tabulation,” Appl. Sci. Research B, v. 6, 1956, p. 144-154. 


[L, SJR. B. Dince, D. Arnot, & S. K. Roy, “The integrals 


C(x) = (p= [ole + at) terde 
and 
Dj(x) = (p> [" o(e + 2)terede 
and their tabulation,” Appl. Sci. Research B., v. 6, 1956, p. 155-164. 
The integrals are tabulated as follows: U,(x) for p = —.5(.5)4, B,(x) for 
p = 0(.5)4, €,(x) for p = —.5(.5)5, D,(x) for p = 0(.5)6.5 all for 
x = 0(.1)1(.2)2(.5)10(1)20 


to 4S except for D, when x > 18, where 3S are given. It is pointed out that values 
of %, and %, for imaginary arguments can be obtained from the present tables 
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using relations such as %,(ix) = (p + 1)€,41(x) — ix€,(x). A table for A,(—x) 
is promised for later publication. Functions are given to facilitate interpolation 
between orders p for small x. No discussion of interpolation with respect to x is 
given. 

The second paper includes a table of the Fresnel integrals 


C(x) = +f J,(t)dt and S(x) = + Jy(t)dt to 12D for x = 0(1)20. 


This table was prepared using the representations of C(x) and S(x) as sums of 
Bessel functions of half-integral order which are tabulated in BAAS [1]. The 
authors point out an error in this table caused by transposition of digits: 


for Jss2(14) read 0.00437 04731 08. 


The reviewer has checked this value. 

The tables for A, and 8, were computed, using desk machines, by use of re- 
currence relations such as p%, + x%,1 = 1 and B, = x"{1 — (p + x)%,} 
using initial values obtained from W%o(x) = —e*Ei(—x) and 


W3(x) = (2/x)tF[(2x)*] 


where F(t) is the ratio of the tail area of the normal curve to its bounding ordinate 
at t, which has been tabulated (BAAS [2]). The tables for €, and D, can be 
obtained similarly from 


b(p — 1)6, + #6,» = 1 and Dp = 34 7{(D + 1) Gp — (6 — 1)G,} 


and expressions for ©» and G; in terms of Si(x) = f% ¢" sint dt and Ci(x) = 
Sz & cost dt and for €;(x) in terms of C(x) and S(x). The necessity for care in 
the use of the recurrence relations is pointed out. The authors state, in a letter to 
the reviewer, that the results were monitored by calculations of asymptotic repre- 
sentations of %, and ©, for large . 

The reviewer checked a substantial portion of the tables for %, and B, using 
SEAC. It is possible to compute %,, for integral p, by use of a Laguerre quadra- 
ture formula and, for half-integral p, by use of a Hermite quadrature formula after 
a trivial change of variable. We can handle %,, €,, and D, analogously using de- 
compositions into complex partial fractions. The error E in the use of these formu- 
las can be expressed in terms of the divided differences of functions such as 
y?/(y + x). For example, for U,(p > 0), one has 


| E| < tm)? /|T (+2) [+0 0O<t<@ 


where x; is the i-th zero of the Laguerre polynomial of degree m. Thus, e.g., for 
= 2,m = 15, p = 4acrude estimate of the error is 4 X 10-*. Asa sample of the 
accuracy actually obtained, agreement with the present tables was achieved for 
x > .7 in the cases p = Oand p = 1. 
In addition to the tables, the papers treat the following topics in connection 
with the functions %,, B,, €, and D,: physical applications (in particular, their 
origin in the theory of elemental semiconductors), recurrence relations, related 
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differential equations, relationships to other special functions, and asymptotic 
expansions. 
Howarp H. WIcKE 


Sandia Corporation 
Albuquerque, New Mexico 


1. BAAS, Mathematical Tables Report, Cambridge University Press, 1925, p. 221-233. 
2. BAAS, Mathematical Tables, v. 7, The Probability Integral, W. F. Sheppard, Cambridge 
University Press, 1939. 


120(L, S].—M. Dank & S. W. Barer, “The specific heat function for a two- 
dimensional continuum. Numerical values of 


CG. _ 6 [? dt ae 


C @iet@—-i #-i 
[MTAC, v. 9, 1955, p. 191-194] 





This function was computed by using a power series for the range 0 < x < 2.0 
and by using a different expansion for the range 2.0 < x < 16.0. The maximum 
error using seven terms of the power series is less than 0.5 X 10-* and for the 
second expansion it is approximately 2 X 10~-*. The table lists values x = (0.1)9.4 
and 10(.5)16, 5D. 


& Eh TF. 


121[M].—F. OBERHETTINGER, Tabellen zur Fourier Transformation, Springer- 
Verlag, Berlin, 1957, x + 213 p., 24 cm. Price DM 35.80. 


This collection of Fourier transforms contains 109 pages of Fourier cosine 
transforms, 91 pages of Foureir sine transforms, and 6 pages of exponential Fourier 
transforms. Also given are concise definitions of the function symbols used, a 
table of errata, and a short bibliography. The book appears to be the largest 
collection of Fourier transforms published to date. The author states in the preface 
that a considerable portion of the approximately 1800 correspondences collected 
in the volume is new. Fourier transforms of elementary functions occupy about 
two-fifths of the space; the remainder is taken up by transforms of higher trans- 
cendental functions, with the Bessel functions and their associates controlling a 
strong minority. 

Glancing over the book one cannot help being impressed by the enormous 
wealth of formulas stacked up in the book and by the astounding formal dexterity 
which Professor Oberhettinger must have commanded in deriving some of them. 
On the other hand one also notices some shortcomings. There is no indication of 
the sources of individual formulas, conditions for validity are given in many cases 
for real values of the parameters only; there is no system of numbering the formu- 
las that could be used for reference purposes; and the bibliography refers only to 
a few well-known titles. However, one must not forget that all these matters are 
of secondary interest to most of the actual users of such tables and that the ab- 
sence of a complicated scholarly apparatus makes for easier reading. 

The existence of some trivial printing errors (on p. 210 and 213) does not in- 
crease the confidence in the absence of non-trivial errors. With this reservation, 
the printing of the book is excellent; by its neat and clean type it is distinguished 
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from the ugly appearance of the Varitype process used in a recent American 
counterpart. Both author and publisher are to be congratulated on this excellent 
work. 


PETER HENRICI 
University of California 
Los Angeles, California 


122[M ].—C. F. LinpMan, Examen des Nouvelles Tables d’Intégrales Définies de 
M. Bierens de Haan, G. E. Stechert and Co., New York, 1944, 231 p., 27 cm. 
Price $5.00. 


According to the preface, Mr. Lindman has inspected Bierens des Haan’s 
famous tables page by page, checking the integrals against suck. original sources 
as were available to him and recomputing a considerable number of the integrals. 
He collected the results of his inspection in a weighty memoir, which was sub- 
mitted to the Royal Swedish Academy of Sciences and published in 1891. The 
volume under review is an unaltered reproduction of this memoir. Although many 
of Lindman’s comments consist in obvious restrictions on the parameters of an 
integral to insure convergence, there is a surprisingly large number of non-trivial 
corrections. It appears that the majority of the 486 tables which make up Bierens 
de Haan’s collection contain faulty results. 

In our review of the 1957 Hafner reprint of the Nouvelles Tables d’ Intégrales 
Définies [1] we pointed out that the extent and nature of the corrections made in 
that edition (as well as in an earlier 1939 reprint) were not evident. Several spot 
checks now reveal that the results of Lindman’s Examen have not been in- 
corporated. Thus the work under review is an indispensable companion volume 
also to the modern editions of Bierens de Haan’s Tables (An earlier review of 
Lindman’s Examen appeared in MTAC [2)]). 7 

PETER HENRICI 


University of California 
Los Angeles, California 


1. D. BrerNEs DE Haan, Nouvelles Tables f Tenterelee Définies, Hafner Publishing Co., 
New York, 1957. [Review 60, MTAC, v. 11, 1957, p. 11 

2 CHRISTIAN FREDRIK LINDMAN (1816-1901), tak des Nouvelles Tables d’Intégrales 
Définies de M. Bierens de Haan, Amsterdam, [sic] 1867. (K. Svenska Vetenskaps Akad., Hand- 
lingar, v. 24, no. 5, Stockholm, 1891.) New York, G. E. Stechert & Co., 1944. [Review 167, 
MTAC, v. 1, 1944, p. 321-322.] 


123[P].—LeEonarD Pope, Tables for Computing the Equilibrium Configuration 
of a Flexible Cable in a Uniform Stream, Report 687, The David W. Taylor 
Model Basin, Washington, D. C., 1951, 31 + 192 p. of tables, 27 cm. This 
report will be available from the Photoduplication Service, Library of Con- 
gress, Washington 25, D. C. 


[P.J—Leonarp Pope & Lours RosENTHAL, Cable Function Tables for 
Small Critical Angles, Supplement to Report 687, The David W. Taylor 
Model Basin, Washington, D. C., 1955, 105 p., 27 cm. 


(A) The equations describing the tension in and the configuration of a flexible 
cable immersed in a uniform steady stream and lying entirely in a plane are de- 
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rived. They are 














cos ¢ : 
decatiad fey t RRS, 
” q(¢) . 
? T 
"= Jeale) 
T COS ¢ 
c= [7 q(¢) de, 
?-rsing 
=| q(¢) de, 


where 
—sin g| sing | + weos¢. 


q(¢) 


Here ¢» is the value of ¢ at the point chosen as the origin of the coordinate system, F 
is the drag per unit length of cable when the cable is parallel to the stream, R is 
the drag per unit length of the cable when the cable is normal to the stream, f = 
F/R, T> is the tension in the cable at the point chosen as origin of the coordinate 
system, T is the tension in the cable at the point x, y,r = T/T», = Rx/To, and 
g= Ry/To. 

These functions called the cable functions are tabulated to four decimal places 
in this report. 

The values of the parameter f for which the cable functions have been evalu- 
ated are 0.01, 0.02 and 0.03. Another important parameter entering into the evalu- 
ation of these functions is ¢., the critical angle of the cable, the value of the angle 
¢ obtained when the cable is freely trailed in the stream and the angle for g(¢) = 
0. The values ¢- = 0°(5) to 85° are covered in the tables. These were computed 
by use of Simpson’s rule for the evaluation of integrals. The relative errors in the 
functions tabulated are said to be less than 0.001 percent. 

(B) The cable function tables for small critical angles reported herein are 
supplementary to those described above. The supplementary tables provide 
tables for critical angles in the range from 0 to 10 degrees in increments of one 
degree. More closely spaced intervals of the independent variable in the vicinity 
of the critical angle ¢. are also provided. Numerical integrations were made using 
the formula 


zy h h 
/ om ydx = 3 [yn + Ayn-1 + Yn—2 ] oy 90 [yn = 4Vn-1 + 6Yn—2 — 4yn-3 + Vaal. 
The authors state that “It is believed that the maximum error in any of the tabu- 
lated values is never greater than one unit in the least significant figure.” 

x. i. Fs: 


124[P, K, X, Z].—E. F. Breckenpacu, Editor, Modern Mathematics for the 
Engineer, McGraw-Hill, New York, 1956, xx + 514 p., 23 cm. Price $7.50. 


This collection of essays is based on a series of lectures organized in the Ex- 
tension Division of University of California at Los Angeles and given at that 
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University and repeated elsewhere. Those who did not hear the speakers will be 
grateful to them, and the editor, E. F. Beckenbach, for providing this volume. 
Engineers are some of the greatest customers of the automatic computers and 
much of the present volume is relevant in the field of this journal, and all our 
readers can profit from reading those essays. 

The essays are divided into three parts: (1) Mathematical Models, (II) 
Probabilistic Problems, and (III) Computational Considerations. The contribu- 
tors and titles are: R. Weller, Introduction; S. Lefschetz, Linear and nonlinear 
oscillations; R. Bellman, Equilibrium analysis: the stability theory of Poincaré 
and Liapunov; J. W. Green, Exterior ballistics; M. R. Hestenes, Elements of the 
calculus of variations; R. Courant, Hyperbolic partial differential equations and 
applications; M. M. Schiffer, Boundary-value problems in elliptic partial differ- 
ential equations; I. S. Sokolnikoff, The elastostatic boundary-value problem; N. 
Wiener, The theory of prediction; H. F. Bohnenblust, The theory of games; G. W. 
King, Applied mathematics in operations research; R. Bellman, The theory of 
dynamic programming; G. W. Brown, Monte Carlo methods; L. A. Pipes, Matri- 
ces in engineering; J. L. Barnes, Functional transformations for engineering de- 
sign; E. F. Beckenbach, Conformal mapping methods; C. B. Morrey, Jr., Non- 
linear methods; G. E. Forsythe, What are relaxation methods?; C. B. Tompkins, 
Methods of steep descent; D. H. Lehmer, High speed computing devices and their 
applications. 

It is not possible to discuss all the contributions; we shall only comment 
briefly on a few which are specially relevant to the field of this journal. 

In his essay in the first part, R. Courant indicates the importance of a double 
attack on problems, by modern computing techniques and by penetrating analy- 
sis. In the second part, the essays by H. F. Bohnenblust and G. W. Brown give an 
account of subjects which, although they have largely developed in the last few 
decades, are becoming more and more used in the programming and operational 
aspects of engineering. 

Some of the articles in the third part are not very closely related with practical 
computational considerations. However, those of Morrey and Forsythe do get 
down to numerical cases. Morrey discusses the solution of functional equations, 
mainly arising in calculus of variation problems. Newton’s method is discussed 
in the finite dimensional case, normed linear spaces are introduced and the Ray- 
leigh-Ritz method is discussed. Examples of the solution of a pair of equations in 
two variables, a first order differential equation and of a calculus of variations 
problem are discussed in detail. 

Forsythe gives a valuable introduction to the subject of relaxation methods, 
indicating the scope of the method and its extensions, giving illuminating exam- 
ples, and discussing the basic questions of convergence which arise. There is an 
ample bibliography. 

Tompkins discusses the method of steep descent in a general framework: 
the solution of many engineering and physical problems can be represented as 
minimizing some function of several variables, or some integral. There are ac- 
counts of the solution of a system of linear equations by the projection method 
of Kaczmarz and the conjugate gradient method of Hestenes and Stiefel. After 
appropriate generalizations, a representative problem of the calculus of variations 








‘sy? 


SS wee OTF |O 





REVIEWS AND DESCRIPTIONS OF TABLES AND BOOKS 285 


is discussed and some account is given of recent work in this field which is relevant 
in connection with computational problems. 

The concluding essay by D. H. Lehmer is a pleasant introduction to automatic 
computers and includes many significant remarks. After noting the inadequacy of 
current machines for handling differential equations involving three space varia- 
bles and time he adds, “‘It is only fair to state that there are also inadequacies in 
the numerical analysis of such problems.”’ He encourages engineers to do their 
own programming and coding—a list of the machine instructions is almost enough 
—and most of the tricks of coding will occur to an engineer “‘who uses his normal 
budget of native cunning.” 

Joun Topp 


California Institute of Technology 
Pasadena, California 


125[(S].—A. Rauman, Tables of Integrals A,(a) = ff d*e-*dd and F,(a) = 
SP Qo(A)A*e-“da, Extrait des Annales de la Societe Scientifique, T. LXIX, 
Serie I, 1955, p. 123-128. 


This paper gives the functions mentioned in the title (wien Qo(A) = $ ln A+ 
to ten significant figures for m = 0(1)10, a = .1(.1)3; 10S. The computations were 
carried out to 12 significant figures. Interpolation formulae are given. The func- 
tions here tabulated are used in calculating the wave functions for the hydrogen 
molecule. 


A. H. T. 


126(S].—_K. M. Case, F. p—E Horrman, & G. PLaczeK, Introduction to the 
Theory of Neutron Diffusion, v. 1, 1953, Los Alamos Scientific Lab., Los Ala- 
mos, New Mexico, viii + 174 p., 26 cm. U. S. Gov. Printing Office, Washing- 
ton, D. C., Price $1.25. 


This book discusses the determination of the distribution of neutrons in space 
and time in terms of the geometrical configuration and physical properties of the 
medium they are in under the restrictive assumption that the magnitude of the 
neutron velocity is unchanged on collision. A very detailed well written treatment 
of this restricted one-velocity theory is given. This theory is basic to the solution 
of more general and more realistic problems. Indeed, some of its results have im- 
mediate application in such problems as the slowing down of neutrons by elastic 
collisions and in other problems. 

The book is organized into three parts: Part A, Introduction; Part B, Propaga- 
tion in the absence of scattering collisions; and Part C, One-velocity theory of 
neutron differsion. 

Part B discusses the solutions of the continuity equation for the angular demity 
function, ¥(r, 2, #) the number of neutrons per unit volume and unit solid angle 
moving in the direction of the unit vector 2; r and ¢ denote the space coordinates 
and time respectively. The continuity equation is discussed for three types of 
streaming: with no sources in vacuum, with sources described by a function 
q(r, 2), and with sources g and absorbers. 
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Included in this part are four place tables for the collision probability for a 
slab of half-thickness ¢ = a/2, a sphere of radius a, and an infinite cylinder of 
radius a. These tables are listed in accordance with the values of b/l, a/l and a/1 
respectively where / is the mean free path. The variable a// (b/1) ranges from 0.00 
to 5.00 in steps of 0.01. 

Part C starts with a discussion of the transport equation and some of its 
properties, discusses this equation for a uniform medium with isotropic scattering 
and finally shows how the results obtained for the uniform infinite medium may 
be applied to the solution of finite problems. In this part there are fairly complete 
tables and graphs of functions entering in the solution of the transport equation. 

The Mathematical Tables Project is credited for the material in Tables 28a 
and 28b. The remaining computations reported in this work were carried out by 
the Los Alamos computing group under the direction of Bengt Carlson and Max 
Goldstein. 

A. Hi. TF. 


127[S}—M. Ferentz & N. RosENzwetc, “Table of F Coefficients,’ U. S. 
Atomic Energy Commission Report ANL-5324, Argonne National Labora- 
tory, Lemont, Illinois, 1955, 293 p. Price $6.30 (photostat), $3.00 (microfilm) . 


This report contains an extensive tabulation of a function F,(L, L’, j’, 7) of 
five arguments known as the F coefficient, which occurs in the formulas for the 
angular correlation between successive nuclear radiations. The arguments of 
the function are restricted in accordance with various triangular and other condi- 
tions. The table is in two parts. Part I contains integral spins (7,189 entries on 
151 pages) and Part II contains the half integral spins (6,413 entries on 133 
pages). The values of F are given to eight decimal places. The computations were 
carried out on the UNIVAC at the AEC computing facility at New York Uni- 
versity and the table was reproduced photographically from the original UNIVAC 
Output. The computations were checked by comparing the results with other 
computations and by computing F(L, L’, j’, 7) (=1). On the basis of these 
checks it is asserted that the entries in the table are correctly given to the sixth 
decimal place. 

A. Et. TF: 


128[S].—A. RauMAN, “Two centre integrals arising out of 2s and 2p atomic func” 
tions,”’ Acad. r. de Belgique, Cl. d. Sciences, Memoires, v. 14, €asc. 2, publ. 
No. 1660, 1955, 13 p., 29 cm. Price 10 Belgian francs. 


In working with the wave functions for diatomic molecules from 2s and 2p 
atomic wave functions centered on the atomic nuclei certain linear combinations 
of integrals previously tabulated by Kopineck [1, 2, 3] turn out to be very useful. 
These combinations of Kopineck’s integrals are tabulated herein to five places. 


A. H. T. 


i. HERMANN-JosEF Kopineck, ‘‘Austausch- und andere Zweizentrenintegrale mit 2s- und 
2p-Funktionen,”’ Z. Naturforschung, B. Sa, 1950, p. 420-431. 

2. HERMANN-JOSEF KoPpiINEcK, ‘‘Zweizentrenintegrale mit 2s- und 2-Funktionen II,” 
Z. Naturforschung, B. 6a, 1951, p. 177-183. 

3. HERMANN-JOSEF Kopineck, ‘‘Zweizentrenwechselwirkungsintegrale III,”’ Z. Naturfor- 
schung, B. 7a, 1952, p. 785-801. 
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129[S, V].—SverrRE PETTERSSEN, Weather Analysis and Forecasting, v. 1, Mo- 
tion and Motion Systems, Second Edition, McGraw-Hill Book Co., Inc., New 
York, 1956, xix + 428 p., 23 cm. Price $8.50. 


As the author states explicitly, this book is designed for use as a general text 
on weather forecasting. Accordingly, its subject matter and style of presentation 
are strongly slanted toward the interests and needs of the practicing meterologist, 
and will be of little direct value to the student of mathematics and numerical 
analysis. 

Two chapters, however, deal with computing methods that are interesting, 
if only for their novelty. In chapter 3, the author outlines a method for extrapolat- 
ing the position and intensity of certain definite features of atmospheric pressure 
patterns (e.g., “high,’’ “lows,” “fronts,” etc.) from one moment to the next. 
This method is based entirely on the differential geometry of the pressure field 
and its instantaneous time-derivatives (assumed known). Thus, since it contains 
essentially no element of physics, it is equally applicable to other time-dependent 
fields. 

Chapter 19 describes a graphical method, originally due to Fjortoft, for solv- 
ing an equation of the form 


Vy: + VN, — VY. = 0 


with y constant along a fixed closed curve in the (x, y) plane at all times ¢, and 
with y given everywhere within that curve at ¢ = 0. Interpreted geometrically, 
this equation states that the instantaneous local velocity of the V*y-field is a vec- 
tor equal in magnitude, but normal to Vy. This provides the basis for a graphical 
method of extrapolating the V*y-field at a slightly later time ¢, by moving the 
initial V*y-field for a short time ¢ with the initial local velocity. The remaining 
problem, that of recovering the y-field at time ¢ from the extrapolated V*y-field, 
is solved by a graphical method of repeated averaging—a procedure that is roughly 
equivalent to Richardson’s method of relaxation. Although such graphical meth- 
ods are doubtless very ingenious, it is the reviewer’s opinion that they are a step 
backward against the current trend toward the more easily automatized digital 
type of computation. 

Chapter 18, on Numerical Forecasting, was specially prepared for this volume 
by Dr. Arnt Eliassen of the Oslo University. This section of the book is a clear, 
concise, and remarkably complete account of the application of numerical methods 
and high-speed automatic computing techniques to the problem of weather pre- 
diction, and is probably of more general interest than the sections discussed above. 

On the whole, Professor Petterssen’s new book is clearly written, well illus- 
trated, and represents a truly heroic effort to bring the practicing weather fore- 
caster up to date with new and improved tools of his profession. 


Puitie D. THOMPSON 


United States Weather Bureau 
Washington 25, D. C. 
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130[S, V].—E. L. Harris & G. N. PATTERSON, The Boltzmann H- Function A pplied 
to the Shock Transition, Institute of Aerophysics, University of Toronto, 
UTIA REPORT No. 40, 1956, iv + 17 p. + 8 p. of figures + 6 p. of tables, 
28 cm. Available from Institute of Aerophysics on an exchange basis only. 


The authors evaluate an approximation to the Boltzmann H-function for a 
gas with a shockwave present, where 


H(z) =n | $(&, x) log nfat, 


n is the number of molecules per unit volume, and nf (x, t)dxdé is the number of 
molecules with velocity components in the range £; to £; + dé; in the volume con- 
taining the point between x; and x; + dx; (t = 1, 2, 3). The function f(é, x) is 
approximated by using the solutions of the hydrodynamical equations of conserva- 
tion of energy, means and momentum. Numerical solutions of the latter equations 
are involved in the calculations. The results of these calculations are listed for 
Mach numbers 1.5, 2.00, 2.3238, 3.00, 3.3764 and 4.00. 
A. H. T. 


131[X].—-Cornetius Lanczos, Applied Analysis, Prentice Hall, Inc., Eaglewood 
Cliffs, N. J., 1956, xx + 539 p., 21 cm. Price $9.00. 


This book is at first an abbreviated account of the author’s work in different 
fields of numerical analysis. Therefore the reader may take it as a guide to the 
papers published by Lanczos in various periodicals. The specific contributions of 
Lanczos to numerical analysis can be listed as follows. First of all he discovered 
the usefulness of Chebyshev polynomials in many fields where they had not been 
used before. The reader will find much valuable material in this connection, in 
particular for solving large linear systems of equations by relaxation, for solving 
algebraic equations, for telescoping power-series and for the approximation of 
solutions of linear differential equations (‘‘r-method”’). The author’s method of 
computing eigenvalues of matrices by ‘‘minimized iterations’’ is outlined a little 
sketchily in chapter III and his more recent invention of “spectroscopic eigenvalue 
analysis” is described. 

But the book gives moreover a general view of numerical analysis. It is very 
interesting for the specialist who will find even in its more elementary parts 
Lanczos’ personal attitudes and his sometimes very original ideas. It‘may however 
be a little hard to read for the beginner. Chapter I deals with algebraic equations 
and much weight is given to Bernoulli’s method called the ‘‘method of moments.” 
Stability is discussed by mapping the left half-plane onto the unit circle. One 
should observe that by this procedure the nearest root to the imaginary axis does 
not correspond necessarily to the nearest root to the unit circle. Chapter II begins 
with a theoretical introduction in linear algebra and proceeds to the numerical 
methods of solving linear equations and eigenvalue problems. Iterative methods 
are given in Chapter III. 

Chapter IV is devoted to harmonic analysis. One of the merits of the book is 
the extensive discussion of the numerical methods connected with network analy- 
sis and Laplace transforms. In particular the author develops four methods for 
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numerical inversion of Laplace transforms. The chapter is a little heavy because 
some results are discussed twice in the languages of Fourier and Laplace 
transforms. 

The more classical problems of data analysis (interpolation and smoothing) 
and quadrature are treated in Chapters V and VI. The reader will find some new 
and interesting points of view, for instance the improvement of Simpson’s rule by 
end corrections and a very elegant treatment of Hermite’s quadrature formulas 
based on Legendre polynomials. 

Error estimates are sometimes a little vague. As an example the estimate 
6-5.4 for Simpson’s rule can be made exact by multiplying the right side by 5/4 
and assuming that the fourth derivative of the integrand does not vanish inside 
the interval of integration. 

At the end of chapter VI there is an interesting application of Hermite’s 
quadrature to eigenvalue problems connected with linear differential equations. 

The last chapter of the book gives various methods for the expansion into 
power series and analytical extension. 

It seems to the reviewer that the modern literature on numerical analysis has 
been taken into account up to 1952 approximately. Throughout the whole book 
almost every computational method is illustrated by numerical examples. 

E. STIEFEL 
Swiss Federal Institute of Technology 


Institute for Applied Mathematics 
Zurich, Switzerland 


132[X].—HEtnz RutisHAusER, Der Quotienten-Differenzen-Algorithmus, Birk- 
hauser Verlag, Basel, Switzerland, 1957, 74 p., 23 cm. Price DM 8.50. 


The Quotient-Difference Algorithm in its simplest form is a numerically ef- 
fective method for determining the poles of a rational (or meromorphic) function 
from one of its Taylor expansions. Based on ideas due to Hadamard, Aitken, 
Lanczos, and Stiefel, and on concepts of the theory of continued fractions, the 
algorithm was developed by Rutishauser in three articles in Z. angew. Math. 
Phystk [1, 2, 3] and shown to be an efficient tool for the solution of a number of 
fundamental problems of numerical analysis. 

Each of the mentioned articles forms the backbone of one of the chapters of 
the work under review, but there are many changes and additions. Chapter I 
furnishes the theoretical basis of the algorithm. Among the new material we men- 
tion the didactically valuable rhombus rules (§3) and an addition theorem for 
continued fractions (§10). The latter appears to be a contribution to the formal 
theory of continued fractions and has several applications in the two later chap- 
ters. Chapter II deals with applications of the Quotient-Difference Algorithm to 
the summation of slowly converging series, to the factorization of polynomials, 
and to the interpolation of functions by sums of exponentials. Much of the new 
material included here tends to bring the algorithm closer to the computer. One 
of the disadvantages of the original form of the algorithm is that it occasionally 
may require dividing by small numbers. It is now shown (§9) how this can be 
circumvented by a simple device. Also, the quadratically convergent modification 
of the algorithm given in the original article for determining the real roots of a 
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real polynomial is now extended so as to yield also conjugate complex roots (§8). 
The original method for interpolation is now replaced by a numerically more 
stable process (§10). The section dealing with the summation of series is perhaps 
the least impressive of the chapter. No attempt is made to define a class of series 
for which application of the algorithm may speed up convergence. Chapter III 
describes the application of the algorithm to the determination of the eigenvalues 
and eigenvectors of a not necessarily symmetric or hermitian matrix. Again the 
author shows his concern with actual computation by including detailed flow dia- 
grams and (as in the other chapters) numerical examples. In three appendices 
the author discusses several extensions of the algorithm. One of these, the so-called 
L-R-algorithm, is treated more fully by Rutishauser in an article to appear in 
volume 49 of the National Bureau of Standards Applied Mathematics Series. 

On the whole the presentation follows the pattern of the underlying papers in 
Z. angew. Math. Physik. A few proofs that were originally missing (e.g., of 
Theorem 1 of chapter III) are now carried out in full. The algorithm is still based 
on the analytic theory of continued fractions, and the reader’s task is not made 
easier by the fact that no references are made to Perron’s standard works on the 
subject. [In an article which will appear also in the above-mentioned volume 49 
of the NBS Applied Mathematics Series, the reviewer has given an introduction 
to the algorithm which is based solely on the elementary theory of complex varia- 
bles. ] The printing of the whole work is excellent, and the fraction bar has ample 
opportunity to demonstrate its much-neglected faculty to make complicated 
formulas easier to read. 

The world of computation owes to the author already several pioneering con- 
tributions on automatic coding and on the numerical integration of ordinary dif- 
ferential equations. By his impressive work on the Quotient-Difference Algorithm 
and its ramifications, Rutishauser has firmly established his position as one of 
the most resourceful and inventive numerical analysts of the present. 


PETER HENRICI 


University of California 
Los Angeles, California 


1. HEInz RuTIsHAUSER, “Der Quotienten-Differenzen-Algorithmus,” Z. angew. Math. Physik, 
v. 5, 1954, p. 233-251. 

2. HErnz RuTIsHAusER, ‘‘Anwendungen des Quotienten-Differenzen-Algorithmus,” Z. angew. 
Math. Physik, v. 5, 1954, p. 496-508. 

3. Heinz RutisHAusER, ‘‘Bestimmung der Eigenwerte und Eigenvektoren einer Matrix mit 
Hilfe des Quotienten-Differenzen-Algorithmus,” Z. angew. Math. Physik, v. 6, 1955, p. 387-401. 


133[X, P].—G. Tempe & W. G. BickLey, Rayleigh’s Principle and its A pplica- 
tions to Engineering. The Theory and Practice of the Energy Method for the 
Approximate Determination of Critical Loads and Speeds, Dover Publications, 
Inc., New York, 1956, vi + 157 p., paper bound. Price $1.50. 


The authors state: ‘‘The object of this book is to explain and justify a rapid 
method of determining the fundamental period of a vibrating system or the con- 
dition of stability of an elastic system, with the degree of accuracy usually de- 
manded in engineering problems. . . . The fundamental principle is due to the 
third Lord Rayleigh, and it applies not only to vibrating systems with a finite 
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number of degrees of freedom, but also to continuous systems such as a stretched 
string or metal reed.” 

The authors only mention that the Rayleigh Principle is intimately related to 
the principle of least action and state that whenever the differential equations of 
a problem, dynamical or otherwise, are equivalent to a variational principle, the 
problem is always soluble by methods analogous to those discussed in this book. 
However, this intimate relationship is nowhere discussed nor is the relation be- 
tween variational principle for the system under discussion and the Induction 
Function G(x, s) (the Green’s function of the equations describing the system 
mentioned). 

An outstanding feature of this book is the use of this function G(x, s) to 
generate sequences of functions f,(x) which approach the first (and higher) 
proper function of the equations describing the system and the use of these appro- 
ximate proper functions for the determination of approximations to the proper 
values of the system. The major portion of the discussion is devoted to a method 
for determining the first proper function and proper value. However lower and 
upper bounds are given for the latter quantity and the lower bound is shown to 
depend on the ratio of the first two proper values. A method for estimating this 
ratio is also given. 

The introductory chapter gives a general account of Rayleigh’s Principle, of 
the problems to which it can be applied, and of the results obtained in subsequent 
chapters. 

Chapters I, II, V, and VI contain many illustrative examples in which the 
methods derived in some detail in Chapters III and IV are applied. Chapters I 
and II are mainly concerned with problems in elastic stability and Chapter V 
deals with vibrating systems. 

Chapter VI deals briefly with theorems relating the fundamental frequencies 
of two systems which themselves are related by inequalities in their mass distribu- 
tions, or potential energies. 

Chapter VII entitled Numerical and Graphical methods discusses one problem: 
The determination of the fundamental frequency of a tapered Aeroplane strut. 

Chapter VIII gives a few examples of the determination of equilibrium con- 
figurations of elastic systems. 

The book is extremely well written and will be particularily useful to those who 
desire closed formulas for various approximations in engineering problems. Mod- 
ern numerical analysts and other mathematicians will find this a thought-pro- 
voking book and very suggestive of methods for solving engineering problems. 
However it is an open question as to how suitable the methods discussed in this 
book are if numerical solutions are to be obtained by the use of automatic digital 
computers. 

A. H. T. 


134[X].—J. Greenstapt, ‘‘A method for finding roots of arbitrary matrices,” 
MTAC, v. 9, 1955, p. 47-52. 
The author proposes to triangularize an arbitrary matrix by a sequence of 


transformations of the form 
A—-B = StAS 
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where S is a unitary matrix operating on two rows and columns and S* denotes 
the transposed matrix. No proof is given that the sum of the squares of the abso- 
lute values of the sub-diagonal elements decreases to zero. It is known that this 
sum does not decrease monotonically. A variety of matrices whose elements were 
picked at random were treated by this method on an IBM 701 and it was shown 
that for this class of matrices the method converged. 


A. Ti. o 


135[X].—E. BopEewic, Matrix Calculus, Interscience Publishers, Inc., New York, 
1956, xi + 334 p., 23 cm. Price $7.50. 


Experience shows that at least half the problems of a modern computing 
organization involve “linear algebra.’’ Computers, in particular, will therefore 
welcome Matrix Calculus, which represents the first attempt to cover the whole 
field in a single volume. 

The book is in four parts. Part I, in 85 pages, treats the theory needed in 
practical applications. It considers elementary properties of vectors and matrices; 
errors in A caused by small uncertainties in A; measures of magnitude of A; 
decomposition theorems and orthogonalization; properties of latent roots and 
latent vectors; bounds for eigenvalues and determinants; and some theory of 
elementary divisors. 

Part II, in 100 pages, considers practical methods for solving linear simultane- 
ous algebraic equations. Direct methods include the elimination and condensation 
methods of Gauss and Jordan, both for exact and approximate solution; applica- 
tions of matrix decomposition into triangles and its connections with elimination; 
the Gauss-Doolittle process, the methods of Banachiewicz and Cholesky, and Ait- 
ken’s ‘‘triple-product”’ method. Details are given of desk-computing techniques 
and their variations of advantage for more automatic equipment. There is also a 
discussion of ill-conditioned equations and the ‘error’ (accuracy) and “‘exactness” 
(precision) of solutions. 

Iterative methods discussed include those of Gauss and Seidel, the accelera- 
tion devices of Aitken, and the steepest-descent methods_of Hestenes, Stiefel, 
and others. The general theory of iteration, with theorems on convergence and its 
acceleration, is given at length. Methods for improving an approximate inverse, 
procedures which always converge and which are therefore considered suitable 
for automatic computing, and the practical advantages of various scaling opera- 
tions, are also included. 

The inversion of matrices, in about 50 pages, is the topic of Part III. In addi- 
tion to the obvious applications and extensions of the methods of Part II there is 
a discussion of partitioning and bordering of matrices, an old but little known 
method for reducing a matrix to diagonal form, and a lengthy account of special 
treatment of the matrices arising in geodetic surveying. 

The last 100 pages belong to Part IV, the determination of latent roots and 
vectors. Iterative methods are used to determine the dominant root and associated 
vector, and perhaps one or two sub-dominant solutions if convergence is slow, in 
this respect an advantage. For accelerating convergence, or causing convergence 
to a wanted solution, methods of “vector-deflation” and ‘“‘matrix-deflation”’ are 
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fully explained: The construction of the characteristic equation is also discussed, 
together with its ill-conditioned dangers. There is a mention of Aitken’s sequences 
for “higher” eigenvalues, and a further discussion of elementary divisors. Then 
come the orthogonal transformation methods of Jacobi and Givens, perturbation 
methods for improving accuracy, the gradient methods of Hestenes and Karush 
and the use of Rayleigh’s quotient. 

Some direct methods are also mentioned, associated with the names of Le- 
verrier, Krylov, Duncan, Hessenburg, Samuelson, and others. These methods ef- 
fectively form the characteristic equation, and find favor only for matrices of low 
order. 

There is a bibliography of 111 references and a somewhat terse index. 

Numerous examples illustrate the text (even in Part I), the amount of com- 
putation in each method is carefully assessed, and there are many computing 
hints, close attention to checks, and a willingness, refreshing if not always believed, 
to give an opinion on the relative merits of the different methods, according to the 
nature of the matrix, the amount of information required about it, and the 
computing equipment available. 

Few omissions come readily to mind. Perhaps the “escalator’’ method was 
worthy of passing mention, however, and the absence of the name of Lanczos, 
from the section on latent roots, is somewhat surprising. 

This reviewer did not find the book easy to read. The emphasis on a “new 
matrix calculus,” with stress on rows and columns, avoidance of the use of in- 
dividual elements, and accompanying new ideas and notations, together with the 
vast amount of material contained in relatively few pages, cannot readily be as- 
similated from the depths of an arm chair. But there is nothing essentially difficult, 
and the student who really makes the effort with this book and succeeds, can claim, 
with its author, to be an expert in this field. 


L. Fox 
National Physical Laboratory 


Teddington, Middlesex 
England 


136[X, I].—H. Fisuman, “Numerical integration constants,”” MTAC, v. 11, 
1957, p. 1-9. 


Tables of the 6; and x; in the approximate quadrature formula 
1 m 
[water =  bielx,) 
= 


are given for m = 0(1)5, m = 1(1)8 to 12D. They were prepared on an IBM 650 
and are an extension of tables of P. C. Hammer, O. J. Marlowe, and A. H. Stroud 
[1] which covered the range m = 1, 2, m = 1(1)5 and m = 3, m = 1(1)4, but 
give 18D. The x; are the zeros of certain orthogonal polynomials g,,(x) defined by 





qm(x) = W/(m + 2m + 1) x (—1)*(™ ome 2 ‘) (7) = 
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which are given explicitly for the same range of m, n. The 6; are defined by 
m—1 —] 
b; = {y [aux F} . 


Joun Topp 
California Institute of Technology 
Pasadena, California 


1. Preston C. HAMMER, O. J. MARLoweE, & A. H. Stroup, ‘Numerical integration over 
simplexes and cones,”” MTAC, v. 10, 1956, p. 130-137. 


137[X, Z].—F. J. WeyL, Report on a Survey of Training and Research in Applied 
Mathematics in the United States, Monograph No. 1, Society for Industrial 
and Applied Mathematics, Philadelphia, Pennsylvania, 1956, vi + 58 p. 
Price $2.00. 


Durch die Erfindung der Rechenautomaten und die schnelle Entwicklung der 
numerischen Analysis ist der vorher ruhig dahinfliessende Strom der mathema- 
tischen Forschung in den letzten Jahren etwas turbulent geworden und wir stehen 
heute vor manchen Problemen grundsatzlicher Natur hinsichtlich des Verhalt- 
nisses der Mathematik zu Naturwissenschaft und Technik und hinsichtlich des 
mathematischen Unterrichts. 

Das dringende Bediirfnis nach Neufestlegung der Standorte und abklarender 
Aussprache hat in den vereinigten Staaten zum vorliegenden Bericht gefiihrt, der 
auf Grund von zwei Konferenzen und einer allgemeinen Umfrage bei den Uni- 
versitaten, der Industrie und den militérischen Organisationen von der mathe- 
matischen Abteilung des national research council ausgearbeitet wurde. 

Die Aufgabe war erstens, Charakter und Umfang der gegenwartigen Forschung 
in angewandte Mathematik festzustellen und nach potentiellen Méglichkeiten 
zur Vertiefung zu suchen und zweitens, Anregungen fiir die Ausbildung ange- 
wandter Mathematiker zu geben. 

Der Bericht umgrenzt zunachst den Begriff ‘‘angewandte Mathematik,” der 
ja sehr schwankend gebraucht wird, und erklart sehr zutreffend, dass angewandte 
Mathematik weniger ein Sachgebiet kennzeichnet als eine geistige Haltung, 
charakterisiert durch die bewusste Bereitschaft zur Zusammenarbeit mit Schwes- 
terwissenschaften, mit dem grossen Ziel, unsere Umgebung und uns selbst besser 
zu verstehen. Mit Recht wird dariiber geklagt, dass viele der fiihrenden Mathe- 
matiker, welche die Atmosphare des Hochschulunterrichts besttmmen, nicht 
imstande sind, Substanz und Ziel dieser Denkweise zu begreifen. 

Er stellt weiter fest, dass die nachste Zukunft durch eine niedagewesene 
Mathematisierung gekennzeichnet sein wird, die sich nicht nur auf Physik und 
Technik beschrankt, sondern auch alle unsere Lebensgewohnheiten beeinflussen 
wird. Als Griinde dafiir werden angegeben einmal die wachsende Kraft mathe- 
matischer Methoden, die Notwendigkeit feinerer Berechnungen in allen Zweigen 
der Technik und ferner das Erscheinen der Automaten. 

Es werden die vielen verschiedenen Wege beschrieben, die von den einzelnen 
Universitaten und Amtsstellen beschritten worden sind, um der steigenden 
Nachfrage nach angewandter Mathematik gerecht zu werden. Vor allem ist hier 
die neue Erscheinung des mathematischen Instituts zu nennen, dessen Erfolge 
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mehr auf der Aktivitat einer Arbeitsgruppe beruhen, als auf spektakuldren 
Leistungen Einzelner. Es wird auch durchaus hervorgehoben und als positiv 
gewertet, dass die Regierungsstellen und die militarischen Instanzen durch ihre 
Kontrakte oft diese Institute erst lebensfahig machen. Wiederum beklagt sich 
der Bericht dariiber, dass viele Mathematiker diese nationalen Notwendigkeiten 
nicht einsehen und Kritik iiben statt am Aufbau mitzuwirken. Es wird auch auf 
die drohende Gefahr hingewiesen, dass die mathematischen Wissenschaften in 
zwei Teile—den abstrakten und den angewandten—auseinanderfallen. 

Das untersuchende Komitee findet aber, dass die Entwicklung zu langsam 
und zu unausgeglichen vor sich geht. Vor allem fehlt es an qualifizierten Stu- 
dierenden, am Unterricht und an einer klaren Konzeption desselben. Die beun- 
ruhigenden Nachwuchsprobleme werden analysiert und auf verschiedene Quellen 
zuriickgefiihrt. 

Aehnliche Unzulanglichkeiten werden hinsichtlich der Rolle des Mathemati- 
kers in der Industrie festgestellt. Obwohl! die Industrie in hohem Mass Mathe- 
matik konsumiert, stellt sie selten Mathematiker an, sondern iibergibt mathe- 
matische Aufgaben ihren Ingenieuren. 


Der Bericht gipfelt in folgenden Empfehlungen: 


1) Es soll ein standiges Komitee gebildet werden, das die Beziehungen 
zwischen der Mathematik und ihren Anwendungen férdert. 

2) Die Universitaten sollen Forschungsprogramme in angewandter Mathe- 
matik fiir fortgeschrittene Studierende durchfiihren. 

3) Alle Mathematikstudenten sollen gezwungen werden, sich auch mit an- 
deren Wissensgebieten bekannt zu machen, in denen die Mathematik 
gebraucht wird. 

4) Mathematiker aus der Industrie oder aus staatlichen Organisationen sollen 
zu Gastvorlesungen an Universitaten eingeladen werden. 

5) Alle Dozenten der Mathematik sollen angehalten werden, Kurse in ange- 
wandter Richtung zu geben, sofern sie dies imstande sind. 


Ein zweiter Teil des Berichtes befasst sich mit einigen typischen neueren 
mathematischen Entwicklungen. Es wird zunachst die axiomatische Methode 
geschildert, sodann die Rolle der Wahrscheinlichkeit und Statistik in den Theorien 
der mathematischen Physik auseinandergesetzt und die Bedeutung einer flexiblen 
und anpassungsfahigen Anwendung mathematischer Methoden in der Technik 
besprochen. Die ganz jungen Gebiete des “operational research,’’ der mathe- 
matischen Theorien biologischer und soziologischer Phanomene werden gestreift 
und die numerische Analysis als Ganzes ist untersucht. 

Als Anhang zum Bericht findet man eine Analyse der Antworten auf den vom 
Komitee versandten Fragebogen. 

Es wurden hier absichtlich diejenigen Teile des Berichtes besonders gewiirdigt, 
die unabhangig von den spezifisch amerikanischen Verhdltnissen weltweite 
Giiltigkeit haben. Alles Gesagte trifft wértlich auch auf europdische Verhaltnisse 
zu und zwar in verscharftem Mass, da in manchen europadischen Staaten auch 
noch das aktive Interesse der Behdrden fehit und somit die finanziellen Mittel oft 
knapp sind. Ferner veranlassen die besseren Arbeitsverhaltnisse in Amerika sehr 
viele und oft die besten jungen Mathematiker aus Europa zur Auswanderung. 
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Der Wunsch nach einer griindlichen Neuorientierung des mathematischen 
Unterrichts hat einige europdische Lander veranlasst, internationale Besprechun- 
gen durchzufiihren und die Fragen gemeinsam zu regeln. Eine erste Konferenz 
fand im Juni 1957 in Miinchen statt am mathematischen Institut der dortigen 
technischen Hochschule. Vertreter von Deutschland, Frankreich, Oesterreich 
und aus der Schweiz arbeiteten eine Resolution aus, die sich weitgehend mit den 
Forderungen des amerikanischen Berichtes deckt und in mehreren Fachzeit- 
schriften (z.B. Zeitschrift fiir angewandte Mathematik und Physik, Heft VIII/5) 
verdffentlicht wird. Die an der Tagung gehaltenen Vortrage sollen ausserdem vom 
Miinchner Institut herausgegeben werden. 

Allen Fachleuten ist es wohl klar, dass angewandte Mathematik heute haupt- 
sachlich ein Erziehungsproblem ist. Es handelt sich weniger darum, immer schnel- 
lere Automaten zu bauen, als das Personal heranzubilden, das imstande ist, 
mathematische Probleme schnell genug fiir den Automaten vorzubereiten. 


E. STIEFEL 
Eidg. Technische Hochschule 
Institut fiir angewandte Mathematik 
Ziirich, Switzerland 


138[Z ].—Mitton H. Aronson, Editor, The Computer Handbook, First Edition, 
55 p. The Computer Handbook, Second Edition, 1956, 71 p., The Instruments 
Publishing Co., Pittsburgh, Pennsylvania, 28 cm. 


The first edition of this booklet contains six papers presented at the First 
Automation Exposition (New York, Nov. 29-—Dec. 2, 1954) and published in 
the magazine Instruments and Automation. The second edition reprints three of 
these and adds eight others. In both cases about half of the papers are on analog 
computers and applications, the other half being on digital computers and appli- 
cations. 

The articles are written for the general reader, rather than being highly tech- 
nical. They provide a limited amount of information about a number of specific 
computers and their business or industrial applications. The articles will not be 
useful to any person desiring detailed technical information about the computers 
discussed. 


Dovucias L. HoGAn 
12512 Epping Court 
Silver Spring, Maryland 


139[Z ].—Proceedings of the Symposium (Automatic Coding), 1957, Monograph 
No. 3, J. of the Franklin Institute, Philadelphia, Pa., 1957, vii + 118 p. Copies 
may be obtained from AUTO CODE, Franklin Institute, Philadelphia 3, Pa. 
Price $3.00. 


The monograph contains the following papers: 


AUTOMATIC CODING AT G. E. by Richard M. Petersen 

SYSTEMS OF DEBUGGING AUTOMATIC CODING by Charles Katz 

PRINT 1—AN AUTOMATIC CODING SYSTEM FOR THE IBM 705 
by Robert W. Bemer 
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THE PROCEDURE TRANSLATOR—A SYSTEM OF AUTOMATIC 
PROGRAMMING by Henry Kinzler and Perry M. Moskowitz 

OMNICODE—A COMMON LANGUAGE PROGRAMMING SYSTEM 
by Russell C. McGee 

A MATRIX COMPILER FOR UNIVAC by Laurence C. McGinn 

A MATHEMATICAL LANGUAGE COMPILER by Alan J. Perlis and 
Joseph W. Smith 

A MECHANIZED APPROACH TO AUTOMATIC CODING by E. C. 
Yowell 


There is a short introduction by John S. Burlew of the Franklin Institute, a 
transcription of the discussion of each paper, and a list of those attending the 
Symposium. 

It seems to be a sad fact that workers in the field of automatic coding, who in 
general are scientific or research personnel, bring a most unscientific and un- 
researchlike attitude to these symposia. The papers presented at this symposium 
are, with two and possibly three exceptions, simply descriptions of individual and 
special-purpose coding systems which are of interest to users of the particular 
computer for which they were written, and to nobody else. Most of the information 
contained in these papers could more easily and efficiently have been obtained by 
interested parties through a users’ organization; by this time such organizations 
have grown up about most large computers. Evidence of the inefficiency of com- 
munication is the turn generally taken in the discussions; questions are often 
requests for information already given in the paper, and comments generally 
consist of destructive or petty and obvious criticism. Constructive criticism—the 
contribution or generation of constructive ideas—is so rare as to be negligible. 

Some of the papers deserve at least individual mention; one such is the Katz 
paper. This contains a short discussion of debugging problems that arise in auto- 
matically coded programs, which are difficult of solution because, by the very 
nature of automatic coding, the programmer is insulated from (and may be com- 
pletely unfamiliar with) the computer. Mr. Katz describes an attempt to produce 
debugging routines for some compilers for UNIVAC. Automatic debugging is a 
most interesting and difficult logical problem, and this paper is one of the few 
in which more rather than less detail would have been welcome. 

The Perlis-Smith paper is another refreshing exception. It is, in the final analy- 
sis, a description of an automatic coding scheme for the IBM 650, but the language 
and organization are a pleasure to see; further, the approach to the problem is at 
once scientific and comprehensible. There is a clear, concise description in logical 
terms of just what the automatic coding problem is, and a clear and convincing 
justification of the methods used in solution. It may or may not be significant that 
this is the only paper arising from work at a university. 

The last paper, and the last deserving mention, is the Yowell paper. This is 
quite different from the others in that it describes an attempt (by The National 
Cash Register Company) to design a computer (the 304) with macro-commands; 
this is automatic programming by hardware rather than routines. The paper is 
at once wordy and short, but the novelty of the content makes it interesting. 
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The tenor of papers and comments makes it clear that the overwhelming 
majority of those attending the symposium have a distorted idea of just who and 
what automatic coding is for. Time and again the emphasis is on simplifying pro- 
gramming so that cheaper talent can be used to write programs. This is a case of 
being right for the wrong reasons. The saving in time and money in making do 
with low-priced programmers is as nothing when compared with the saving that 
can be effected if high-priced talent—scientists, engineers, executives—can have 
access to a computer without having to explain problems or learn programming. 
The only statement or comment made at the Symposium in support of this 
philosophy comes from Dr. Perlis and is, in part: ‘‘It is the man who proposes the 
problem who is important, and his ideas are the ones that count.” 


E. P. NEUBURG 


Route 2 
Silver Spring, Md. 


140[Z].—D. D. McCracken, Digital Camputer Programming, John Wiley and 
Sons, Inc., New York, 1957, vii + 253 p., 24 cm. Price $6.50. 


Digital Computer Programming offers for the first time an excellent introduc- 
tion to the problems of coding for a high speed stored program computer. Such a 
work has long been needed in this field; McCracken’s book will probably remain 
as the standard of comparison for some time. Since this book will undoubtedly be 
widely used, not only as a textbook in colleges but also as a training manual within 
industry, it is felt that a fairly detailed review of the work is indicated. In what 
follows, the plan of attack will be to cover the book more or less chapter by chap- 
ter, and, in general, point up some of its shortcomings. Most of these arise from 
viewing the book from the point of view of a neophyte in coding, and will possibly 
seem obvious to the old hands in the business, who can read Digital Computer 
Programming aimost as if it were a novel, exclaiming ‘‘Yes, yes!’’ at every page. 
Following this survey will come a few remarks on the general philosophy of the 
book. 

In chapter 1 the newcomer is given a brief introduction to stored program 
computers. It becomes obvious almost immediately that the work is intended as 
a textbook, as new terms are introduced with fair rapidity. Most of these terms are 
indicated in italics, and defined contextually, which makes for easy reading, but 
these terms must be mastered before going on, and will probably require several 
readings on the part of the student. In order to have a concrete machine with 
which to work, McCracken invents TYDAC, the TYpical Digital Auto- 
matic Computer, a decimalized version of a modern binary scientific machine. 
A brief history of computing follows, which is fairly complete with one exception. 
Most of the major users of large scale machines for scientific computing started 
out by applying standard punched-card machine techniques to engineering cal- 
culations, and many of the methods, traditions, and customs which the newcomer 
to scientific computing encounters have their roots in this early work. Some credit 
should be paid to these pioneers. 

A fairly comprehensive outline of the steps in coding skips mention of one 
of the most important steps—that of making sure the proposed problem is really 
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the one which should be solved. Often a great deal of labor may be saved by in- 
corporating several small problems into one large systems analysis. On the other 
hand, it may be too early in the investigation of a problem to tackle more than a 
small piece at a time. The large amount of time spent recoding problems points up 
a need for more time spent in this phase of solution. Also, no mention is made of 
the large amount of feedback that ordinarily occurs between the scope of the 
problem under solution and the method employed. 

There is a good survey here of typical applications of digital computers; this 
might also have included some examples of problems which are not suited to 
high-speed stored program machines. 

A lucid, though conventional, explanation of coding fundamentals is presented 
in chapter 2. The going is still fairly fast here, and much reference will be made 
back to the terms previously defined in chapter 1, and to the appendix which con- 
tains a complete summary of TYDAC. (A worthy addition to the book would 
have been a small card, loosely inserted, containing a summary of the instruction 
codes, etc. Practically every coder for a real machine carries one of these supplied 
by the manufacturer; TYDAC should not be an exception.) Cogent points which 
must be remembered are often hidden in the middle of a paragraph; the book will 
probably be much underlined. 

A good glossary of computing terms would be welcome here. Since this text 
will find such widespread usage, it is to be regretted that the author did not in- 
clude what might have become a definitive dictionary of computing terms. 

There is a good description of the execution of a three step program which 
adds two numbers, but the section on multiplication and division processes within 
the machine is probably overdone. A better use of space here would have been to 
present the arithmetic registers as they appear before and after the execution of 
the instructions, as is done in the section on scaling. It is pleasant to see an intro- 
ductory work which mentions the necessity of clearing the MQ register before 
division takes place, although more emphasis might have been placed here— 
especially in connection with small-integer arithmetic. 

In an example on p. 26 several opportunities have been missed to point up the 
possibility of subtle errors in even such a simple problem as (10-A + B)C. Here, as 
happens throughout the text, the examples seem fairly simple when the correct 
answer is given. Only three examples of coding in the entire text are pointed out 
as containing errors—certainly a far smaller proportion than maintains in real- 
life coding. A very worthwhile addition to the explanations would give examples 
which do not work, and the reasons why. 

The exercises appended here, as throughout the work, are excellent. Perhaps the 
most frequently neglected part of teaching coding is in this realm, and it is re- 
freshing to see examples more typical than A + B = C. 

Chapter 3 leads us painlessly through the mysteries of binary and octal num- 
ber systems. One addition to this chapter that could be made would explain the 
use of mixed-base number systems for compacting information within a word, as 
is often done to minimize input-output time or conserve storage space. 

The Decimal Point Location Systems, as described in chapter 4, are exceed- 
ingly well done. McCracken is to be congratulated on his fearless approach to the 
so-called ‘‘mysteries” of fixed-point coding. There are many experienced coders 
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who could profitably study this explanation. The first method, that of fixing the 
point in the middle of the word, seems rather awkward and artificial, and should 
probably have been cut a bit. Here would be a good place to explain the necessity 
of clearing the MQ before division, but this is not done. 

The “graphic’’ method has long been the favorite of this reviewer, who finds 
here the first really good explanation of the method. The examples of multiplica- 
tion miss the fact that in a large proportion of cases a left shift is necessary after 
the multiply, since both factors may not reach their maximum value at the same 
time. Also neglected is the problem of finding a full ten digit rounded quotient, 
which is never presented in the book. The scaling of the remainder after a division 
is never presented; this should be included. Some simplification could have been 
made by considering only the number of digits before the point, as the number 
after is redundant. To be complete, this method should include a set of rules 
similar to the 5 step rule given under the scale factor method, and would probably 
give the ‘“graphic’”” method more widespread use. 

The commonly used scale factor method is accompanied by a good set of rules 
for its employment, and should dispel much fear on the part of the coder approach- 
ing fixed point problems. 

In the section on scaling, the author has missed an opportunity to present the 
coder with a check on the programming of the problem that is often valuable. In 
coding fixed point arithmetic, the values of the shifts required should normally 
be small; e.g., 3 or 4 places at the most. If a shift of eight or nine is arrived at in, 
say, an addition, an insignificant quantity is probably being added. The state- 
ment of the problem may be checked at this point, and many errors have been 
detected in such a fashion. 

The coder should be warned, in conjunction with the problem of overflow, 
that in a large majority of cases the overflow indicator will be turned on by some 
part of the problem which does logical work, or by subroutines, and that it should 
be turned off prior to embarking on a sequence of fixed point arithmetical calcula- 
tions and checked immediately following the sequence. 

In this section, as elsewhere, a good deal more space could have been profitably 
devoted to pointing out possible pitfalls. The examples are generally good, al- 
though more requiring “tight” scaling to retain significance could have been in- 
cluded. 

Chapters 5 and 6, which deal with address computation and loops, lead the 
coder quite easily into the real power of stored program machines. Particularly 
to be commended is the decision to separate the address computation from the 
testing computation; introductory examples here are often much too com- 
plicated. Two objections were noted here. In several examples some of the neces- 
sary information is buried in the text, and it would have been better to write out 
cells containing instructional constants, etc., following the examples. Where this 
is not done, the somewhat confusing notation “‘Loc 1400” is used rather than the 
more common “‘L(1400)”’ or “Loc. of 1400.” 

This section would be a good place to point up the advantages of standardiza- 
tion when writing loops. There are many ways to write a given loop; the important 
point is for the coder to pick a single method to use and stick with it—the tempta- 
tion to be clever is seldom resisted by the neophyte, and usually leads to disaster. 
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The excellent coder usually develops his own style of coding, and this becomes a 
major factor in reducing his own mistakes. 

There is an excellent paragraph at the bottom of page 77 which bears repeating 
here. 

‘Depending on the nature of the loop and of the test, it is possible to make a 

truly remarkable variety of mistakes in testing. If the loop should be carried 

out m times, it is quite easy to make mistakes which will result in doing it: 

(1) not at all; (2) 2 — 1 times; (3) 2 + 1 times; (4) 2” times; (5) until the 

power fails or the machine breaks down. It is fairly safe to say that loops, 

although one of our most powerful tools, are also a very large source of errors. 

Whatever other prechecking systems may be used, it is always advisable to 

go back and check the loop-testing parts of the program. As indicated above, 

one simple way is to analyze what would happen if the loop were to be exe- 
cuted only once.” 

This succinctly conveys to the coder the main source of errors in coding. 
McCracken might have pointed out that once a coder learns to “count to one”’ 
(not only in writing loops, but in many other phases of a problem) he is well on 
the way towards eliminating mistakes in coding. 

On page 80 there is a discussion of the time vs. space problem which always 
seems to bother the new coder. Perhaps this would be a good time to inject the 
point that very often the best loop is not the shortest, not the fastest, but the one 
which works the first time. The amount of time spent in debugging a slightly 
faster or shorter loop often far overshadows the time saved in running the problem. 
Many coders adopt the policy of writing all loops in a standard fashion, checking 
the program out, and then going back and “cleaning up the code,” one loop at a 
time, when time considerations are really important. 

There is also a gem of wisdom at the beginning of chapter 6 which might have 
been emphasized more—in coding a loop, start from the inside and work out. This 
philosophy can be carried to all parts of a problem, and bears repeating many 
times, although the student may quickly find this out for himself in working out the 
fine examples at the end of this section. 

The chapter on flow charting is well introduced and presented; unfortunately 
no attempt has been made to give general rules for problem organization. More- 
over, no symbolism has been adopted for indicating the execution of a closed sub- 
routine, whose detailed flow chart may be given elsewhere. Some attention to 
this point might have helped to keep flow charts from growing almost beyond 
comprehension. 

TYDAC contains, as do most modern scientific computers, a set of index 
registers. A rather arbitrary, though adequate, set of instructions is provided for 
operating on them, and an extremely good explanation is made of the distinction 
between indexable and non-indexable instructions. Another example of a program 
with an error occurs here, and again the opportunity to point up the advantages of 
standardized coding is missed. 

In chapter 9 the student is introduced to one of the most powerful tools at his 
disposal—subroutines. The various methods for entering and leaving a subroutine 
are well presented, but two methods for communicating information to the sub- 
routine that are commonly used are not emphasized. (These are transmittal via 
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the machine’s arithmetic registers and via standard cells in memory.) Several 
other points are missed that are vital to the discussion. A generally used subrou- 
tine must somehow take care of all possibilities of error in inputs, with appropriate 
error warnings, and must also conform to some prescribed set of conventions re- 
garding the use of overflow indicators, index registers, etc. An example of a sub- 
routine with an error exit is given in one of the examples in a later chapter, but the 
point is not emphasized. It is in the attention to details such as these that a good 
subroutine is distinguished from a poor one. 

The statement is made, on p. 114, that it is not clear whether or not the use of 
index registers for linking to subroutines will grow. On the contrary, the reviewer 
feels that one of the most important advantages of the newer machines lies in their 
built-in features, such as the set index jump instruction and indirect addressing 
features, which permit subroutines to be constructed with great ease while coding. 
The closed nature of these subroutines greatly facilitates the check-out phase, 
and they will undoubtedly become a more and more important part of program- 
ming. It is at this point in the book that the coder might have been introduced to 
good coding practices, and it is unfortunate that the point was not made. 

The introduction to floating point operation, as given in chapter 10, is generally 
good. His warning against the injudicious use of floating point is not too strongly 
made, in the face of the general trend toward hardware incorporating this feature. 
The new coder may be lulled into a false sense of security if he is not strongly con- 
fronted with the pitfalls that can arise. One of the most time-consuming problems 
for a programmer using built-in floating point operations in a large problem con- 
cerns the treatment of underflows and overflows; the routines necessary are not 
mentioned at all. It should also have been explained that built-in floating point 
operations are not always exactly analogous to the corresponding fixed point 
operations; e.g., TYDAC’s MQ is cleared in floating point addition. This is often 
a big source of errors for the new coder. 

Input-output methods can be treated only briefly in a book such as this, as 
McCracken points out, since these are strongly dependent on the actual machine 
used. This is unfortunate, as perhaps the best way for a coder to learn the dif- 
ferences between numbers, addresses, locations, operations, etc., and at the same 
time meet non-numeric and numeric data and conversion and scaling problems is 
to tackle an input or output routine. This is especially true on a binary machine. 
The statement that a coder need never know input-output will be contested by 
many, including this reviewer. The inclusion of a console on TYDAC leaves some- 
thing to be desired in view of the general trend away from having the programmer 
touch the console, or indeed, operate the machine at all. There might have been a 
section included, more pertinently, on the problems of organizing the production 
phase of a problem so that it may easily be run by a professional operator who 
may have no knowledge at all of the coding or the problem itself. 

A load program is presented which uses an initial address and a word count on 
each card to control loading and the advantages are detailed. The disadvantages 
should have been presented, too. Suppose the deck has had corrections added, 
and is somehow out of order? There is no indication that this has happened, and 
the program might run through tocompletion, giving wrong answers, at a large 
cost in machine time. Similarly, no indication is given if a card is missing; much 
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running time can elapse before this is detected. Many problems arise in connection 
with input-output procedures, and the programmer should be made aware of 
these early in his career. 

The presentation of magnetic tape procedures is again necessarily brief. The 
inclusion of check-summing as a machine checking feature is good, especially in 
view of the fact that TYDAC does have a redundancy check indication. The 
necessity of having programmed checks on the machine operation, even when 
checking circuits are built in, could probably have been emphasized elsewhere 
with equal advantage. Here again the coder should be warned that his problem 
will be run by professional operators unfamiliar with his routine and that he should 
try to simplify tape handling procedures as much as possible. 

Chapter 13 presents a basic introduction to program checkout techniques. 
Since the average coder spends the largest amount of time in this phase of problem 
preparation, a much more extensive treatment of the subject would be indicated. 
Unfortunately, no general set of rules to guide the programmer has as yet been 
presented. McCracken fails to distinguish the two phases of program checkout: 
debugging and testing. By debugging is meant the process of making sure that the 
program does what the coder meant it to do, and by testing is meant the process of 
making sure that the program solves the problem it is intended to solve. The dis- 
tinction is strong, and the testing phase is often woefully neglected, even by ex- 
perienced programmers. 

The debugging process is fairly well covered, including the standard techniques 
of tracing, memory dumping, and breakpoint printing. Also welcome is a descrip- 
tion of the “‘decoding”’ process of checking a code before it is tried out on the ma- 
chine, and various methods are presented. An excellent example points up the 
necessity for picking realistic and non-degenerate test case data, and should 
decrease the number of times the cry, “But the test case checked perfectly!” is 
heard. 

No rule is given for the procedure to be followed when the newly coded pro- 
gram collapses; a simple one is to try to get the program to execute all of the in- 
structions from beginning to end without regard to the correctness of answers, and 
then to tackle the arithmetic portion of the code. All such rules are, of course, 
rather dependent upon the programming system used, and will vary from ma- 
chine to machine. Another fruitful area for making mistakes is ignored—the 
procedure for correcting mistakes once they are found. Nowhere is this problem 
covered, and as the neophyte programmer does not have symbolic techniques at 
his disposal at this point in the text, he may reach the conclusion that it is neces- 
sary to do a great deal of rewriting each time an error is discovered. A few words as 
to the necessity for leaving room for patching would be in order here. Also, were 
symbolic methods available to the coder at this point, the ease with which these 
permit the inclusion of personalized diagnostics could be presented. 

It is difficult, in a short’chapter, to present relative (or symbolic) programming 
techniques and give a good indication of their power. The advantages of these 
techniques are most strongly felt on large problems, and a student is not inclined 
to go through a large problem example in a new and unfamiliar notation. It is 
significant, however, that once a coder has mastered relative and symbolic coding 
systems he usually becomes much more skilled in programming actual machine 
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language. The second exercise in this chapter requires the student to find a 
“‘small’’ error in a routine which has been presented in the text, and points up the 
fact that even advanced methods lead not to perfect programs, but to different 
types of errors, with which the coder must learn to cope. 

The final sections on interpretive routines, double precision routines, and mis- 
cellaneous programming techniques are well written and the examples which ac- 
company these chapters provide the opportunity for the student to apply the prin- 
ciples learned in the early part of the book. 

In spite of the general level of excellence of this text, there are several areas 
which have been rather neglected if programming, rather than coding, is to be 
taught. The first of these is that of program checkout and testing, and the second is 
that of problem organization. In the majority of classes using this as a text, it is 
difficult to see how the student can be made aware of the many mistakes he will 
undoubtedly make in coding the examples given. Certainly few instructors will 
have time to discover, by the process of decoding, all of the subtle mistakes which 
an inexperienced coder can make in even the simplest examples. One solution 
which will undoubtedly be adopted at a few schools, is to code an interpretive 
routine for a high speed computer to simulate the features of TYDAC, and to 
provide laboratory sessions in which to try out the problems. 

If this is not practical, or if a high speed machine is not available, the only 
alternative would seem to be to standardize coding techniques and penalize the 
student who uses fancy routines much as he would be penalized in actually trying 
to check them out on an operating machine. Certainly, as has been indicated in 
individual cases above, there are many opportunities in the text to point out the 
various mistakes that can be made, and it would be a great help if a greater por- 
tion of the text, and particularly the examples, was devoted to finding errors in 
routines. There is a trend, particularly among the users of the largest scientific 
machines, towards finding mistakes without the use of elaborate machine debug- 
ging programs. Examples might be given which would contain the original coding 
for a routine, and the results of running this coding on TYDAC. (This routine 
stopped at xxxx. Why?) If necessary, the relevant portions of a memory or console 
dump could be supplied. 

Perhaps the greatest difference between the performance of a mediocre coder 
and a top-notch programmer lies in the way in which they organize a problem. 
It is a far cry from writing a set of twenty or fifty instructions to completing a 
problem containing many thousand orders, and many coders fail to bridge this 
gap. The good programmer thinks at all times about the possibility of errors, and 
so arranges his problem to eliminate as many mistakes before they occur as he can. 
One technique that is widely used is to break the problem into small, self-contained 
subroutines, trying at all times to isolate the various sections of coding as much as 
possible. At first glance this appears extremely wasteful of instructions—many of 
them are communicating information between routines that could well communi- 
cate directly with one another. It is in this direct communication, however, that 
many mistakes are made, and by isolating the various portions of the code as much 
as possible the problem is reduced to many much smaller ones. The truth of this 
seems very obvious to experienced coders, yet it is hard to put across to the new- 
comer. An extremely valuable addition to the McCracken text would be a rather 
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large problem, perhaps requiring as many as several thousand instructions, broken 
down into small logical blocks, with descriptions of these blocks. Similar problems 
could be presented for analysis by the student. This sort of work also has the ad- 
vantage of being quite easily checkable by the instructor. 

The experienced prograinmer, on reading the text, will probably wonder, as 
did the reviewer, why symbolic coding was not introduced much earlier rather 
than being relegated to a few pages at the end of the book. Certainly everyone who 
is introduced to computer programming in industry is given symbolic coding 
shortly after he learns the basic operation of the machine. The aim of symbolic 
coding is to help alleviate some of the burdens of coding, help him in checking 
and testing his program, and it is probably most useful to start the would-be coder 
off with as much help as possible, leaving him to concentrate on the essential 
problems of programming. Many excellent programmers would be lost without 
such an aid, and the student is probably no exception. 

One final criticism of the book will be raised by many. Why, it will be asked, 
bother to teach many students to program for a mythical machine which they 
will never actually use, and bears only a confusing relationship to the machine 
they will ultimately program, and, additionally, one with so many instructions? 
This reviewer feels, on the contrary, that it is probably worthwhile for the coder 
to be familiar with more than one machine, even if one of them is mythical. 
Machines are continually being manufactured which are larger, faster, and with 
a larger set of instructions, and the switch from TYDAC to an actual machine will 
adequately prepare the coder for future machine changes. 

In conclusion, it should be stated once more that Digital Computer Program- 
ming will probably be the standard text on coding for some time to come, and 
rightfully so. 

C. L. BAKER 


The RAND Corporation 
Santa Monica, California 


141[Z ].—W. Stuckin, Mind and Machines, Penguin Books, Inc., Baltimore, 
Maryland, 1954, 221 p., 18 X 11 cm. Price $0.50. 


This book attempts a comparison of human brains with digital and analogue 
computing machines in the endeavor to answer some questions raised by the exist- 
ence of these machines. These questions: of similarity of mind and machine, of 
how the human brain works, of what constitutes “thinking,” of whether machines 
“think,” are all open fields for controversy. The author presents a many-sided 
discussion, leaving it to the reader to make his own decisions, and without always 
disclosing where he himself stands on some of these problems. 

Cybernetics, and the concept of negative feedback are discussed in an attempt 
to determine how the human brain functions. Machines, particularly electronic 
machines, use the negative feedback principle extensively. How valid are our at- 
tempts to introduce this concept into our discussions of the nervous system? The 
author shows clearly the value and the limitations of the principle of negative 
feedback in this application. 

In one sense, however, he persistently ignores a fact which to this reviewer 
should be more often used in these arguments, particularly on the question of 
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whether what a computing machine does constitutes real thinking. This fact is 
that man has built the machine and has put into it, by virtue of his designs, a 
predictable pattern of behavior. True, this may be a very complicated pattern, 
it may involve decisions simulating human thought processes, it may even simu- 
late “‘learning.’’ But in every case, in the final analysis, a machine has been built 
which can only do what the designer says it can do. 

Thus it seems a little odd to take this simulator, this computing machine, and 
use it as a model of what we suppose the brain to be like. It would be very surpris- 
ing if it did not resemble the brain in some sense, since we have built it with our 
own logical thought processes as a model. This seems to be a weakness of the argu- 
ments for a mechanistic interpretation of the human brain. 

One minor criticism concerning usage should be made. The author states, on 
page 176, for example, that digital and analogue computing machines “think” 
in that they perform mathematical calculations. This is not correct in a strict sense. 
These machines perform arithmetical calculations. 

This is an informative book, and while it avoids dogmatic answers to the 
problems it states, it does give a wide basis for consideration of validity of the 
many hypotheses which have been suggested. 

W. Sluckin combines interests in two fields ideally situated to help in the dis- 
cussion of these problems: he is both a psychologist and an electrical engineer. 

Each chapter concludes with a pertinent bibliography for additional reading. 


F. H. HOLLANDER 


University of California 
Los Angeles, California 


142[Z].—Elektronische Rechnenmaschinen und Informationsverabeitung (Elec- 
tronic Computers and Information Processing), Nachrichtentechnische Fach- 
berichte (Technical Information Report), Friedr. Vieweg and Sohn, Braun- 
schweig, Germany, v. 4, 1956, viii + 299 p., 29.5 cm. Price DM 26. 


This volume reports the proceeding’s of an international meeting held in Darm- 
stadt, Germany, from October 25th to 29th, 1955 under the sponsorship of the 
German Association for Applied Mathematics and Mechanics (GAMM), the 
German Society of Electrical Engineers (VDE), the German Mathematical 
Society (DMV), and the Association of German Physical Societies (VDPG). 

Among the seven introductory papers, three were presented by U. S. repre- 
sentatives with H. H. Goldstine evaluating computer speeds required to solve 
significant differential and integral equations by existing methods; A. S. House- 
holder discussing numerical mathematics; and H. H. Aiken projecting the future 
of automatic computing machinery. A. D. Booth of England discussed input- 
output systems. The problems of data processing, programming and switching 
and memory techniques were respectively considered by R. Piloty of the Tech- 
nische Hochschule Miinchen, H. Rutishauser of the Technische Hochschule, 
Ziirich, and H. Billing of the Max Planck Institute, Gottingen. 

The largest group of papers concerned themselves with descriptions of Euro- 
pean computer developments and will have value as source material not readily 
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found elsewhere. The computers reported on are tabulated below: 


COUNTRY 
Austria 


Belgium 


Czechoslovakia 


Germany 


Netherlands 


Sweden 
Switzerland 


U.S. A. 
USSR 


ORGANIZATION 


Institute for Low Frequency 
Techniques, Technische 
Hochschule, Vienna 

Institute for the Encourage- 
ment of Scientific Research in 
Industry and Agriculture and 
National Fund for Scientific 
Research, Antwerp 
ARITMA, Prague 

Institute of Mathematical 
Machines, Academy of Sciences, 
Prague 


Institute of Applied Math. 
Tech. Hochschule, Dresden 
Institute of Electrical 
Information Handling 
Tech. Hochschule, Munich 
Institute for Experimental 
Mathematics, 

Tech. Hochschule, Darmstadt 
Leitz Optical Works, 
Wetzlar 

Max Planck Institute, 
Gottingen 


PTT Laboratories, Hague 
Mathematische Centrum 
Amsterdam 


Matematikmaskinnamnden 
Stockholm 

Institute for Applied Math. 
Tech. Hochschule, Zurich 
IBM 

Academy of Sciences, Moscow 


MACHINE 
Logical Function Computer 


IRSIA-FNRS 


Calculating Punch 
SAPO 


Digital Plotter 

Linear Equation Solver 
Relay Network Synthesizer 
Di 


PERM 


DERA 


Z5 
Gt 


G2 

G3 
PTERA 
ARRA 


ARMAC 
BESK 


ERMETH 
705 


BESM 
URAL 


Eleven papers describe approaches to programming in Germany and the 
Netherlands, with consideration given both to changes in machine design to 
facilitate programming and to compiling and interpretive methods. 

In the area of numerical analysis a dozen papers cover problems in linear alge- 
bra, iterative processes, interpolation, hyperbolic partial differential equations, 











308 REVIEWS AND DESCRIPTIONS OF TABLES AND BOOKS 


quadratures, nonlinear differential equations, hydrodynamical equations, linear 
programming, aircraft design computation, weather forecasting, matrix inversion, 
and astronomy. 

The concluding miscellaneous section includes a survey of logical function 
computers, a description of contact grids as an aid to relay network synthesis, a 
discussion of statistical programs in industry and description of an equivalence 
algebra for representation of digital machines. 

Aside from the content of the papers, interested researchers will find a num- 
ber of new references associated with several of the articles. 

Moreover, the publisher has been quite considerate in providing a list of ad- 
dresses of the authors as well as short translated summaries of all the papers. 


G. EstRIn 


University of California 
Los Angeles, California 


TABLE ERRATA 


Reviews and papers in this issue mention errata in the following works: 


A. J. C. CUNNINGHAM & H. J. WooDALL, Factorization of y* + 1, y = 2, 3, 5, 6, 
7, 10, 11, 12 up to high powers (n), Francis Hodgson, London, 1925. (See 
RAPHAEL M. ROBINSON note, ‘Some factorizations of numbers of the form 
2" + 1,” p. 265.) 


R. B. Dincie, D. Arnot, & S. K. Roy, “The integrals 


A,(x) = (p) [ &(e + x)—'e-de and B,(x) = (p+ | P(e + x)—~e-de 
0 0 
and their tabulation,’’ Review 119, p. 279. 


D. N. Leumer, List of Prime Numbers from 1 to 10006721, Review 107, p. 272. 


257. —C. A. Coutson & W. E. Duncanson, “Some new values for the exponen- 

tial integral,”’ Phil. Mag., v. 33, 1942, p. 754-761. 

A comparison of tables of Ei(x) and —Ei(—x) appearing in this paper with 
more elaborate tables by Harris [1] reveals a total of nine rounding errors and two 
more serious errors in the former. 

Rounding errors appear in the 10-figure values given by Coulson and Duncan- 
son for Ei(x) when x = 20, 31, 46, and 47, and for —Ei(—x) when x = 16, 
20, 23, 29, and 39. 

Their claim of accuracy to within 2 units in the last figure is refuted, however, 
by the following emendations. Corresponding to x = 44 and 45 their approxima- 
tions to Ei(x) should be increased by nearly 8 units in the last place; that is, for 
2.990444711, read 2.990444719, and for 7.943916028, read 7.943916036, respec- 
tively. 

The accuracy of Harris’s table of the interpolation coefficients R,(1) was con- 
firmed by me by an independent calculation, and these data were then used to 
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verify the consistency of his 18-figure values for Ei(44) and Ei(45) with the ad- 
jacent entries in his table. 


Joun W. WRENCH, JR. 
Applied Mathematics Laboratory 
David Taylor Model Basin 
Washington 7, D. C. 


1. Franx E. Harris, “Tables of the exponential integral Ei(x),”” MTAC, v. 11, 1957, p. 9-16. 


258.—CaRL-ERIK FROBERG, Hexadecimal Conversion Tables, C. W. K. Gleerup, 
Lund, Sweden, 1957, (MTAC, Review 82, v. 11, 1957, p. 208.) 
The following erratum has been found: 
Page 11, for 0.2 65 00D41 DF3B6 45A1D 
read 0.2 65 00D4F DF3B6 45A1D. 


CARL-ERIK FROBERG 


University of Lund 
Lund, SWEDEN 


NOTES 
Acknowledgments to Referees 


The editors of Mathematical Tables and Other Aids to Computation were happy 
to take advantage of the action of the Council of the American Mathematical 
Society at their December 1956 meeting in Rochester, New York, and furnish to 
them a list of persons who have refereed papers for MTAC during 1955-1956. 
This list was incorporated with names furnished by various other journals and 
published on pages 27 and 28 of AMERICAN MATHEMATICAL SOCIETY 
NOTICES, April 1957. The editors of MTAC deeply appreciate the services 
rendered by its referees. 

G. we Ee 


The Illinois Journal of Mathematics 


The Illinois Journal of Mathematics is devoted to the publication of original 
research papers in pure and applied mathematics. The first two issues of this 
journal have appeared with the tables of contents listed below. 


Contents, Vol. 1, No. 1 March, 1957 


A fundamental inequality in the theory of valuations, by I. S. Cohen & 
Oscar Zariski 

Kompakte projektive Ebenen, by Hans Freudenhthal 

Homology of Noetherian rings and local rings, by John Tate 

Finite dimensionality of certain transformation groups, by Deane Mont- 
gomery 

On modules of trivial cohomology over a finite group, by Tadasi Nakayama 

Markoff processes and potentials I, by G. A. Hunt 

Inequalities for asymmetric entire functions, by R. P. Boas, Jr. 

On a class of linear differential equations with periodic coefficients, by 
Jack K. Hale 
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On a problem of Picard concerning symmetric compositums of function- 
fields, by Jun-ichi Igusa 
On matrix classes corresponding to an ideal and its inverse, by Olga Taussky 


Contents, Vol. 1, No. 2 June, 1957 


Classes of finite groups and their properties, by Reinhold Baer 

The irreducible representations of a semi group related to the symmetric 
group, by Edwin Hewitt & Herbert S. Zuckerman 

On Ingham’s trigonometric inequality, by L. J. Mordell 

Fonctions aléatoires 4 correlation lineaire, by Paul Levy 

Sur la corissance radiale d’une fonction méromorphe, by Paul Malliavan 


On some applications of dynamic programming to matrix theory, by 
Richard Bellman 


Manuscripts intended for publication in this journal may be written in 
English, French, German, or Italian and should be addressed to the Illinois 
Journal of Mathematics or one of the following editors: 


Reinhold Baer, Mathematisches Seminar der Universitat, Frankfurt am 
Main, Schumannstrasse 58, Germany. 

J. L. Doob, Department of Mathematics, University of Illinois, Urbana, 
Illinois. 

A. H. Taub, Department of Mathematics, University of Illinois, Urbana, 
Illinois. 

George W. Whitehead, Massachusetts Institute of Technology, Cambridge 
39, Massachusetts. 

Oscar Zariski, Department of Mathematics, Haryard University, 2 Divinity 
Avenue, Cambridge 38, Massachusetts. 


A. H. Taus 


Policy Committee for Mathematics, Technical Advisory Committee to the 
National Bureau of Standards, Annual Report, May, 1957 


1. The Bureau of Standards Technica! Advisory Committee for Mathematics 
has reviewed the work of the Applied Mathematics Division of the Bureau 
(Division 11) and finds that in spite of the difficulty this division is having in 
maintaining its present staff and acquiring new staff, it is providing very excellent 
mathematical, computational, and statistical services to other agencies of the 
government as well as to other Bureau divisions. The four sections of Division 11; 
Numerical Analysis Section (11.01), the Computation Laboratory (11.02), the 
Statistical Engineering Laboratory (11.03) and the Mathematical Physics Sec- 
tion (11.04), have interesting, useful and imaginative research and training 
programs which make them effective organizations for the discharge of their 
responsibilities. 

2. The Numerical Analysis section (11.01) has continued its outstanding 
work in modern Numerical Analysis and the use of computing machines in pur- 
suing research problems in mathematics. In addition it has inaugurated a train- 
ing program in numerical analysis for senior university staff. 
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The purpose of this program is to give regular university staff heretofore 
specializing in fields different from numerical analysis, a training in that field 
which will enable them to direct the operation of a university computing center, 
and to organize training and research in numerical analysis on their return to their 
own institutions. Thirteen people began work on this program on February 11, 
1957 and continued to be in residence at the Bureau of Standards until June 7, 
1957. The funds for stipends to the participants and other expenses in connection 
with this training program were supplied by the National Science Foundation. 

Trainees were taught programming for SEAC and were given a very inten- 
sive course of lectures on classical numerical analysis, linear equations and matrix 
inversion, quadrature and ordinary differential equations, characteristic values 
of matrices, partial differential equations, integral equations, linear programming, 
and non-linear equations. In addition lectures were given covering a variety of 
other topics including recursive functions and Turing machines, statistical cal- 
culations, discrete variable problems, and some experimental mathematical cal- 
culations. 

The lecturers included, in addition to members of the staff of the Applied 
Mathematics Division of the Bureau, mathematicians of note presently engaged 
in research on one or more of the topics listed above. 

This training program should greatly stimulate teaching and research in 
numerical analysis in universities. It complements the work being done at some 
universities and does not conflict with it since the universities are in the main 
concerned with the training of students and this program has as its trainees es- 
tablished mathematicians. 

3. The Computation Laboratory (Section 11.02) like other sections of the 
Division, is having difficulty in maintaining its staff. If this situation is viewed 
as an educational service to the community in that it furnishes to industry and to 
other agencies of the government young people with good training in computing 
procedures, the Bureau may well take pride in the role it is playing. 

However, the problem is a difficult one, particularly as the Section experiences 
a need for the expansion of its staff to provide for the increased work it must un- 
dertake with the delivery of the IBM 704 and the obligation to provide some pro- 
gramming help to DOFL. 

The Handbook of Mathematical Tables is progressing well. The Section is to 
be particularly commended for drawing into the work experts in the special phases 
of mathematics that are represented either by securing the temporary services of 
these experts at the Bureau, or by contracting out the preparation of specific parts 
of the book. 

4. Statistical Engineering Laboratory (Section 11.03) has as a major activity 
the consulting with scientists in other divisions of the Bureau on statistical aspects 
of their problems, including design of experiment and analysis of results. It is 
performing this function well, and its services are in great demand. The SEL 
continues to be a center of basic research in experimental design; ten research pa- 
pers in experimental design are in process of publication in various technical 
journals. The reliability study recently undertaken by SEL has made a good be- 
ginning and promises to produce valuable results. 

5. Mathematical Physics Section (11.04) carries on fundamental research on 
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problems of its own, and has the function of providing high level specialist con- 
sulting services to other parts of the Bureau in the field of mathematical physics. 
The committee finds that a good program of work is going forward within the 
limitations of the small staff available. 

The idea of having a central consulting group on difficult problems of mathe- 
matical physi¢s to serve the rest of the Bureau is sound and commendable. But 
it can only function effectively if its staff includes creative scientists of the highest 
caliber. By its nature the work cannot be carried on effectively by persons in the 
lower ranks. It follows necessarily that this section must develop in a way that its 
staff will have considerably higher than average grade ratings as compared with 
the rest of the Bureau of other comparable scientific laboratories. 

The development thus far has been mostly in the direction of fundamental 
work on theoretical mechanics. In the years to come corresponding staff strength- 
ening should occur to permit the rendering of effective consulting service in 
thermodynamics, statistical mechanics, chemical physics, electrodynamics, and 
so on. 

During the past year important contributions were made by the section in 
exact non-linear water wave theory, theoretical elasticity, the theory of elastic 
surface waves, the theory of diffraction and reflection of electromagnetic waves, 
and the study of the Fourier transformations for distribution functions occurring 
in applied statistics. 

The strength of the section in carrying on this program of highly original work 
has been greatly enhanced because the Bureau has been following a wise and highly 
praiseworthy policy of having on its staff a succession of some of the most re- 
nowned mathematical physicists in the world for varying periods as visitors. 
These have included Stoneley of Cambridge, England, Synge of Dublin, Eire, 
Schultz-Grunow of Aachen, West Germany, Burgers of Delft, Holland, and Walz 
of Géttingen, West Germany. 

The field of mathematical physics is one that gains in creative productivity in 
unusual degree from stimulating contacts with great men in the field from other 
lands. The Bureau is indeed to be congratulated, in the opinion of the committee, 
on having developed such a first-rate stimulating atmosphere that men of such 
high distinction find it worth while to spend some time there as visiting research 
workers, and to give the Bureau the benefit of sharing their latest research ideas. 
This phase of the program represents a very genuine step forward which deserves 
every encouragement. ” 

Since further recruiting will be required to maintain the research of this 
section at its present level of excellence, and to extend their consulting service to 
other divisions of the Bureau, it may be useful to review the philosophy behind 
the establishment of Section 11.04, to see how important its maintenance might 
be. Many Divisions of the Bureau deal with various aspects of classical (non- 
atomic) physics and, in order to improve their measurement techniques and 
broaden their service coverage, they require detailed theoretical computations in 
the fields of classical physics. Many of these computations are difficult enough to 
be classed as important research projects in theoretical physics, requiring theoret- 
ical physicists of exceptional ability to solve them successfully. 

One way to get such problems solved is for each Division of the Bureau to hire 
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theoretical specialists in their own field, but this has several disadvantages. The 
number of challenging problems arising in any one Division may not be sufficient 
to attract the services of a first-class man and the general attitude of the Division 
may be such as to discourage theoretical research in new directions. Moreover the 
specialist, being isolated in one Division, may find it difficult to obtain the advice 
and encouragement from other theoretical workers, in other specialties, which 
often can make the difference between failure and success in some difficult 
problem. 

Section 11.04 represents an alternate way of providing such specialized con- 
sulting services. The theoretical physicists are kept as a group, supporting each 
other in their work and able to divert manpower, as needed, to various important 
problems as they arise in the various Divisions. This arrangement utilizes scarce 
manpower more efficiently and makes it easier to attract good men, but it places a 
considerable responsibility on the head of the Section to keep in touch with the 
various Divisions of the Bureau so that the Section can be of real help to them 
(even, at times, to the extent of persuading the Divisions to get interested in an 
area which they had neglected before) . 

If the Bureau is convinced that the maintenance of Section 11.04 is the correct 
way to provide consulting services for the Divisions on difficult theoretical 
problems, it must be prepared to make special efforts to recruit the experienced 
personnel needed by the Section. A Section with these responsibilities, staffed 
by mediocre people, would defeat its purpose. If such recruiting involves the 
hiring of specialists from Europe, special arrangements should be made to make 
this possible. 

Respectfully submitted, 


A. H. Tavs for the 

Technical Advisory Committee 
for Mathematics 

Davip BLACKWELL 

E. U. ConDON 

MarkK Kac 

Puitie M. Morse 

Mina REES 

A. H. Taus, Chairman 


Index 


In this volume, MTAC has continued its policy of trying to review all pub- 
lished tables, other than the most elementary ones, which depend on or contribute 
to computation. In addition, there have been reviews of works, essentially book 
length in size, contributing in some way to numerical analysis or its applications. 
A classified index attempts to furnish the reader with as much information as is 
conveniently summarized about these works. The index includes references to 
tables from number theory, algebra, analysis, statistics and physics, and to works 
or tables relating to numerical analysis, to the social and biological sciences, and 
to computing equipment and its use. 

In the present volume, the author and the nature of the function tabulated 
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are put close together. This is a result of a suggestion by Dr. J. C. P. Miller and 
seconded by others. In the section devoted to analysis, the classification of the 
Fletcher, Miller, and Rosenhead Index has been retained without radical change. 
In the section devoted to statistics, an outline prepared by Professor H. O. Hart- 
ley has been used with considerable augmentation. Also, it should be noted that, 
while analysis tables are assigned fairly consistently to the class to which the 
tabulated function belongs, the statistics index frequently assigns a table to the 
section devoted to the use of the function. Some ambiguity in classification is un- 
avoidable, and multiple entries are common. 

In the tables from number theory, the classification of the Lehmer Index is 
followed without essential change (except for the addition of a section devoted to 
tables pertaining to analytic number theory not suitably assigned elsewhere). 

Detail in the classification of the works on numerical analysis seemed justified 
by the amount of correspondence and conversation which exists concerning ade- 
quate presentations of various types of material ; this has led to voluminous multi- 
ple listings. 

The indexing was done largely by the Chairman of the Editorial Committee, 
but he gratefully acknowledges the help of W. J. Dixon, Rudolph Hiisser, and 
J. D. Swift—the first two in connection with the classification of tables from statis- 
tics, and the last in connection with number theory. 

Cc oe 


CORRIGENDA 


D. R. Morrison, “A method for computing certain inverse functions,” 
MTAC, v. 10, 1956. 


for * read 

p. 205, line -4, 1/2(2* — 1) logs (1 + 2e) 

p. 206, line 14,f ff we 

p. 206, line-11,f" fo | 

p. 207, line -4 5 arc cos (1/y) 5 — arc cos (1/y) 

° te i ° ° € log? 

p. 205, lines —3, —2, omit It is less than 4 we 
H. J. HAvUER 
E. A. Fay 


U. S. Naval Ordnance Test Station 
China Lake, California 


Review 24, MTAC, v. 11, 1957, p. 31, line 3 from bottom, 
Exponent should read r not t. 
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AUTHOR INDEX 


Papers and Technical Notes 
[Short papers and notes are marked (N) in this index] 


Author Title 
BELFORD, G., L. Jack- Table Pertaining to Solutions of a Hill Equation (N) 
son Lastett, & J. N. 
SNYDER 
BoLey, Bruno A. A Method for the Numerical Evaluation of Certain Infinite 
ee ee oe ee ee 
CAMPBELL, Epwin S. A Method for Integrating a Set of Ordinary Differential Equa- 


tions Subject to a Type of Numerical Indeterminacy.... . 
CHAMBERS, LL. G. Note upon the Numerical Evaluation of Limits of Sequences 

a a ee ee a eee Poe 
Cook, J. M. Rational Formulae for the Production of a Spherically Sym- 

metric Probability Distribution (N).................... 
Diyxstra, E. W. A Method to Investigate Primality (N)................... 
DoyLe, THomas C. Inversion of Symmetric Coefficient Matrix of Positive-Defi- 


Rk Aci Ph nes dcdodnencnabeehs 
Farrincton, C. C., On the Numerical Solution of Sturm-Liouville Differential 
R. T. Grecory, & Equations 


A. H. Taus 
FisHMAN, HERBERT Numerical Integration Constants......................... 
Goon, I. J. On the Numerical Solution of Integral Equations (N)....... 


GREENSPAN, DONALD On the Numerical Evaluation of the Stokes’ Stream Function. 
HamMeER, Preston C., Numerical Evaluation of Multiple Integrals I 
& A. WAYNE WYMORE 


Harris, FRANK E. Table of the Exponential Integral Ei(x)................... 
Hart, RoGer G. A Formula for the Approximation of Definite Integrals of the 
Normal Distribution Function (N)....................- 
Heaty, M. J. R. A Rotation Method for Computing Canonical Correlations (N) 
Hitcuock, A. J. M. Polynomial Approximations to Bessel Functions of Order 
Zero and One to Related Functions. (N)................ 
Hotiapay, Joun C. A Smoothest Curve Approximation....................... 
Linc, Car-Binc Tables of Values of 16 Integrals of Algebraic-Hyperbolic Type. 
Lonoman, I. M. Tables for the Rapid and Accurate Numerical Evaluation of 
Certain Infinite Integrals Involving Bessel Functions... . . 
LuKE, YUDELL L. On the Computation of log Z and arc tan Z............... 
MENDELSoBN, N. S. An Iterative Method for the Solution of Linear Equations 
Based on the Power Method for Proper Vectors (N)...... 
MENDELSOBN, N. S. The Computation of Complex Proper Values and Vectors of a 


Real Matrix with Application to Polynomials (N)........ 
Mever, H. I. & B. J. A Method of Inverting Large Matrices of Special Form (N).. 
HOLLINGSWORTH 


Miter, J. C. P. Note on the General Solution of the Confluent Hypergeo- 
ET LE eS | 
Perce, WILLIAM H. Numerical Integration over the Spherical Shell............. 


RoBInson, RAPHAEL M. Factors of Fermat Numbers (N)..............---.+-.-+-- 

Rosinson, RAPHAEL M. Some Factorizations of Numbers of the Form 2* + 1 (N)... 

SaLzeR, HERBERT E. Equally-Weighted Quadrature Formulas for Inversion Inte-- 
a sien ch-<4 ok RinebS 5s 0 Shhh ss 00.00.0509 Mo phaee’ 


Page 
79-81 


. 229-233 


19-21 


81-82 
195-196 


55-58 
131-150 
1-9 
82-83 
150-160 
59-67 

9-16 


265 
83-86 


86-88 
233-243 
160-166 


166-180 
16-18 


91-94 


94-97 


97-99 


21-22 
265 
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Author 
SHELDON, J. W., 
B. ZonpEx, & M. 
FRIEDMAN 
SmitH, HERSCHEL F. 
SouTHARD, THomaAs H. 


Stecun, IRENE A. & 
Mitton ABRAMOWITZ 
Stroup, A. H. 


THACHER, HEnry C.,JR. 


VISWANATHAN, R. V. 
WaGnNER, HARVEY M. 
WARNER, F. J. 

Wir, HERBERT S. 


Wirkes, M.V. & D.J. 
WHEELER 


Title 
On the Time-Step to be Used for the Computation of Orbits 
Dey DI io. oon onc i oc cece ccccvccevces 


On a Generalization of the Prime Pair Problem............ 
Approximation and Table of the Weierstrass Function in the 

Equianharmonic Case for Real Argument (N)............ 
Generation of Bessel Functions on High Speed Computers. . 


Remarks on the Disposition of Points in Numerical Integra- 
a a a: eo ee ee ere eer 
Optimum Quadrature Formulas in s Dimensions. . 
Solution of Poisson’s Equation by Relaxation Method—Nor- 
mal Gradient Specified on Curved Boundaries. . . 
A Lemma for Automatic Optimum Pusgnimming o on : the IBM 
RR ee ia) ee AG ie | 
On the Solution of “Jury” Problems with Many Degrees of 
PON Se OT eer oe FOE EIU Ac lis sono dewdile 
An Open Formula for the Numerical Integration of First 
Order Differential Equations (N)......................- 
Note on “A Method for Computing Certain Inverse Func- 
GR PPE CURA. CoP Ee eas Fes eee eees ccosccgdbatves 


Page 


181-189 


249-254 


99-100 


. 255-257 


257-261 


- 189-194 


67-78 


101-104 


268-271 


201-203 


INDEX OF REVIEWS BY AUTHOR OF WORK REVIEWED 


Author 


Allen, E. E. 
Anderson, Einar 
Arndt, D. 


Armour Research Fdn. Proc. 94 


Aronson, Milton H. 


Banerjee, D. P. 
Barber, S. W. 
Bateman, Harry 
Battin, Richard H. 
Beckenbach, E. F. 
Belevitch, V. 
Belford, G. 

Bell Aircraft Corp. 
Belousov, S. L. 
Bennett, Albert A. 
Bickley, W. G. 
Bishop, R. E. D. 
Bodewig, E 
Brombacher, W. G. 


Case, K. M. 
Chisnall, G. A. 
Churchill, Stuart W. 
Clark, George C. 


Couffignal, L. 
Crandall, S. H. 
Crank, J. 


Crowder, Harold K. 





Review Review 
Number Page Author Number 
36 38 Dank, M. 120 
90 212-214 | Davis, P. 84 
119 279-281 | Davison, W. R. 30 

215-216 | de Haan, D. Bierens 60 
138 296 de Hoffman, F. 126 
Dingle, R. B. 119 
is 62 Dodge, H. F. 32 
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69 117-118 Drew, Thomas B. 40 
63 112 Durand, David 17 
124 283-285 
78 205-206 | Eckert, W. J. 101 
59 110 Elektronische Rechenma- 
86 210 schinen und Informations- 
115 276 verabeitung , 142 
69 117-118 | Emersleben, O. 56 
133 290-291 | Emersleben, O. 65 
89 212 Emersleben, O. 66 
135 292-293 | Erdélyi, A. 67 
= = Fedorova, R. M. 49 
126 285-286 | Ferentz, M. 88 
6 24 Ferentz, M. 127 
116 277 Fieller, E. C. 12 
116 277 Fishman, H. 136 
70 118-119 | Forsythe, George E. 47 
99 219-221 | Francis, George C. 117 
71 119-120 | Franklin Institute Proc. 139 
117 277-279 | Fraser, D. A. S. 16 
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Author 


Fréberg, C. E: 
Fréberg, Carl-Erick 


Garcia, Juan 
Gellman, H. S. 
Gloden, Albert 
Glowatzki, Ernst 
Godwin, H. J. 
Greenstadt, J. 
Greenwood, J. Arthur 
Guttman, Irwin 


Haan, D. Bierens de 
Hammer, P. C. 
Hannan, E. J. 
Harris, E. L. 
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Partial Differential Equations 

44 Papers in proceedings symposium on numerical 
analysis 

53 Tables for approximate solution of Laplace’s equation 

71 The Mathematics of Diffusion 

73 Papers in report of convention on digital computer 
applications 

97 Methods and applications in book on mathematical 
physics : 

98 Elementary discussion in advanced calculus text 

99 Elementary numerical treatment of equilibrium and 
propagation problems in text on engineering 
analysis 

124 Hyperbolic equations and applications, boundary 
value problems in elliptic equations, elastostatic 
boundary value problems, essays in book on mathe- 
matics in engineering 

142 Hyperbolic equations, hydrodynamics, in proceed- 
ings of symposium on electronic computers 
Integral Equations 

43 Monte Carlo methods in proceedings symposium on 
Monte Carlo 

97 Some applications in book on mathematical physics 

124 Prediction theory, essay in book on mathematics in 
engineering 
Special Functions 

67 Higher Transcendental Functions 


Page 


217 


218 


219 


283 


306 


46 


108 
119 
121 


217 


218 


219 


283 


306 


44 


283 


114 
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Review 
1ge Author/s Number Description Page 
Sir Harold Jeffreys & 97 Bessel, hypergeometric, Legendre, I, and elliptic 217 
Bertha Swirles functions in book on mathematical physics 
117 F. G. Hildebrand 98 Elementary discussion, Bessel, Legendre, f, hyper- 218 
geometric functions in advanced calculus book 
121 Asymptotic Expansions 
Sir Harold Jeffreys & 97 Chapter in book on mathematical physics 217 
217 Bertha Swirles 
218 Definite Integrals 
D. Bierens de Haan 60 Table of definite integrals 111 
219 C. F. Lindman 122 Critique of de Haan’s tables 282 
Integral Transforms 
283 Sir Harold Jeffreys & 97  Sometransform applications in book on mathematical 217 
Bertha Swirles physics 
306 F. B. Hildebrand 98 Forty-one Laplace transforms and elementary dis- 218 
cussion in advanced calculus book 
F. Oberhettinger 121 Table of Fourier transforms 281 
E. F. Beckenbach, Editor 124 Functional transformations for engineering design, 283 
essay in book on mathematics in engineering 
46 Cornelius Lanczos 131 Harmonic analysis and related topics sections of a 
book on applied analysis 
108 Sampling, Generating Random Variates, Monte Carlo 
119 NBS AMS No. 12 42 Proceedings of a symposium on Monte Carlo 43 
121 H. A. Meyer, Editor 43 Proceedings of a symposium on Monte Carlo 44 
Proc., Second Annual 94 Machine program for Monte Carlo 215 
217 Computer App. Symp. 
E. F. Beckenbach, Editor 124 Monte Carlo methods, essay in book on mathematics 283 
218 in engineering 
219 
Applications to Chemical Sciences 
Thomas B. Drew & 40 Advances in Chemical Engineering 41 
283 John W. Hoopes, Jr., 
Editors 
Bibliographic Works 
306 A. S. Householder 46 “Bibliography on numerical analysis” 49 - 
George E. Forsythe 47 ‘Selected references on use of high-speed computers 49 
for scientific computation” 
A. V. Lebedev & 49 Guide to tables 104 
R. M. Fedorova 
44 
TABLES FROM NUMBER THEORY 
217 
Euler’s Totient Function and its Inverse, Sum and Generalizations 
283 E. H. Neville 3: Denominators of Farey Series of order 105 23 
Ordinary Factor Tables 
D. N. Lehmer 107. Factors of N,0 < N < 100 17000, 2, 3, 5,7 KN 272 
114 
Tables of Factors of Numbers of Special Form 
Albert Gloden 109 Factors of N* + 1, N = 3001(1)6000 274 














326 TABLES FROM ALGEBRA 


Review 
Author/s Number Description 

R. M. Robinson 111 Least factor f of N = k-2" + 1,k = 1, < 2272, 

f < 105; k = 3," < 1280;k = 5,” < 1536; 

k = 7, < 1280;k = 9(2)99, n < 512:f < 104 

Lists of Consecutive Primes 

D. N. Lehmer 107. _—~ Primes p, 1 S p S 100 06721 
John Leech 110 Groups of 4, 5, 6, 7 primes » with closest possible 

spacing 50 < p < 100 17000 

Lists of Primes of Special Form 

R. M. Robinson 111 Primality of N = k-2" + 1,k = 1, < 2272; 

k = 3," < 1280;k = 5," < 1536;k = 7, 

nm < 1280; k = 9(2)99, nm < 512 

Diophantine Equations of Degree > 2 
D. H. Lehmer 50 Solutions of x? + y3 + 2? = 1 
Non-linear Forms, their Classes, and Class Numbers 

R. J. Porter 113 Irregular negative determinants of exponent 3m with 

their critical classes, 50000 s —D S 1,00000 

Tables Related to Cyclotomy 
Emma Lehmer 112. Distributions of Gauss sums of orders 3, 4, 5, 7 
Tables Relating to Algebraic Number Theory 
H. J. Godwin 51 Quartic fields 
Other Tables Relating to Analytic Number Theory 
O. Emersleben 56 Some values of Epstein’s-zeta-function 6D 
D. H. Lehmer 108 25,000 zeros of ¢(s) 
Tables for Radix Conversion—Decimal to Radix 2" and Radix 
2* to Decimal 

Carl-Erik Fréberg 82 Hexadecimal conversions: integers 1 (1) 1024(16)4096 


and 10*(10*)10**+1, k = 2(1)12; fractions x = 
10-*(10-*) 10-**2, & = 2(2)16; m-10*, m = 1(2)9, 
k = 12(1)12; 10* and 10-*, k = 13(1)25; some 
constants, all in hexadecimal form. 
x = 167*(16-*) 16-*+1, k = 1(1)10 

in decimal form r 

A. van Wijngaarden 83  Decimal-octal conversion x = 0(8)99992; 2", ” = 
1(1)50 exact; 2-*, m = 1(1)50, 20D. Octal 10*, 
nm = 1(1)18 exact and 10, m = 1(1)18 to twenty 
octal digits 


TABLES FROM ALGEBRA 


Tables Involving Marks from Algebras or Arithmetics 
Other than the Usual Marker Systems 


H. J. Godwin 51 Real quartic fields 


Matrices of Real Numbers 


W. L. Wilson, Jr. 53 Inverses to Laplacian operators over triangular grid 
of n(n + 1)/2 points, m = 2(1)7 


Page 
274 


272 


274 
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275 
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208 
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107 


109 


208 


208 


107 
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Author/s 
Richard B. Smith 


D. H. Lehmer 
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Review 
Number Description Page 
95 A, and B,, nm = 2(1)15, where a; = 1,0; = 216 
(¢ +7 — 1)7% = 2(1)mand bi; = (¢ +7 — 1) 
114 Two classes of matrices with elements 0, 1 and 275 
formulas for inverse, powers, eigenvalues 


TABLES FROM ANALYSIS 


Multiplication Tables and Proportional Parts 


Walter Schmidt 


Franz Triebel 


V. Belevitch & F. Storrer 


H. Nagler 
V. Belevitch & F. Storrer 


Hans J. Maehly 
V. Belevitch & F. Storrer 


H. S. Uhler 


J. L{Hodges, Jr. 


S.'C. Nicholson & J. Jeenel 


75 m-n,n = 1(1)200,m = p + @where p = 0(1)100, 204 
6 = .1(.1).9, %, %, 36, %, %, 34, 1D 

76 mn, n = 1(1)100, m = 1(4%4)100, 2D; also m = 205 
101(1) 1000, » = 1(1) 100 exact; also m — 1(1)300, 
n= ¥, 4%, %, 2D 


Reciprocals 
78 1/x, 15S (polynomial approximation) 205 
Square Roots and Related Functions 


77. WN, N = 1(1)10000, 15D 205 
78 x, 15S (polynomial approximation) 205 


Fractional Powers of Real Numbers 


4 2, —4 sx S }, 10D (rational approximation) 24 
78 10° or 10* — 1, 15S (polynomial approximation ) 205 
Factorials and Sub-factorials 
1 996! and 1000! exact 22 
Binomial Coefficients 
22 4+#C,, A = .5(1)19.5, # = 1(1)18, 7S 31 


Constants Involving + 
2 #, 3089D 22 


Logarithms of Integers, Common Constants, etc. 


Frank E. Harris 


Frank E. Harris 


D. Bierens de Haan 
C. F. Lindman 


V. Belevitch & F. Storrer 
Juan Garcia 


118 log we, 24D 279 
Euler’s Constant 7, e*’, and Related Constants 
118 y, 24D 279 
Definite Integrals 
60 Definite integrals 111 
122 Critique of de Haan’s tables 282 
Common Logarithms 
78 log wx, 15S (polynomial approximation) 205 
106 log m, m = 1(1)9999, 5D with radical interpolation 272 
aids 


Logarithms to Non-decimal Integral Bases 


Hans J. Maehly 


4 log sx, 12 S x S 1, 10D (rational approximation) 24 








328 TABLES FROM ANALYSIS 


Review 
Author/s Number Description Page 
Natural Trigonometric Functions with Argument in Radians 
R. B. Horgan 5 Sin xand cos x, x = n-10*,n = 1(1)9,k = —3(1)1 24 
V. Belevitch & F. Storrer 78 Sin x, 15S (polynomial approximation ) 205 
Natural Trigonometric Functions with Argument in Time 
L. S. Khrenov 79 Sin x, cos x, tan x, cot x, sec x, cxc x 206 
x = 0'(4*)34, 5S, Aand P.P; cot x, csc x, 
x = 0(.1*)8"(1*)40™, 8S; sin? (4x), 
x = 0(4*)124, 5S 
Logarithms of Trigonometric Functions with Argument in 
Degrees, Minutes, Seconds 
Juan Garcia 106 log sin @, log cos @, log tan @, log cot @, 272 
0° (various) 10” (.1’’) 140’ (1’’) 16’40” (10’’) 10° 
(20’’) 80°, etc., 5D with radical interpolation aids 
Inverse Circular Functions 
Hans J. Maehly 4 arctanx, —1 Sx S 1, 9D (rational approximation) 24 
V. Belevitch & F. Storrer 78 arctan x, 15S (polynomial approximation) 205 


Tables Concerning Right-angled Plane Triangles 


E. H. Neville 3 arc tan y/x in degrees, x = 1(1)105, y = 1(1)x, 13D 23 
and in radians y = 1(1)105, x = 1(1)y, 15D, r = 
(x* + y*)4, x = 1(1)105, y = 1(1)x, 13D, and 
Inr, y = 1(1)105, x = 1(1)y, 15D 
Natural Values of e . 
Frank E. Harris 118 eh = .01(.01).05, .1(.1).5, 1, 25D 279 
Natural Values of e= 
Frank E. Harris 118 e*,h = .01(.01).05, .1(.1).5, 1, 25D 279 
Logarithms of Hyperbolic Functions 
Homer S. Powley 81 log cosh x, x = 10(.5)20(5)40, 3D 208 
Natural Logarithms of Integers 
E. H. Neville 3. Inx,x = 1(1)160, 15D; kin 10, k = 1(1)20, 15D 23 


Natural Logarithms of Reciprocals, Square Roots, and Other Functions 


M. R. Sampford 21 xinx/l — x, x = 0(.01).99, 4D 30 


Circular and Hyperbolic Functions of Complex Argument, 
Rectangular Values for Argument and Function 


Jgrgen Rybner 80 sinh (A + 4B) and cosh (A + iB),0 5 A <4, 207 
tanh (A + 7B),0 S A S 3 (alignment charts); 
with several auxiliary functions 


Exponential Integrals (Ei(x) 


Frank E. Harris 118 Eji(x), « = 1(1)4(.4)8(1)50, 18S and e*EI(x), 279 
x as above, 19D, with interpolation aids 








Frank 


NBS 


NBS 


Bell 


~ 


en) 
—" 


Ha 











TABLES FROM ANALYSIS 329 
Review 
‘age Author/s Number Description Page 
Exponential Integrals — Ei(— x) 
24 Frank‘ E. Harris 118 —Ei(—x)-x = 1(1)4(.4)8(1)50, 18S and 279 
205 —eEi(—<x), 
x as above, 19D, with interpolation aids 
206 The Gamma Function 
F. B. Hildebrand 98 I(x),x = 1(.01)1.99,4D 218 
The Error Integral, Tables of the Ordinate (the Integrand) 
NBS AMS No. 41 33 H' (x), x = 0(.0001)1(.001)5.6(various), 15D 36 
Tables of the Error Integral 
2 
" NBS AMS No. 41 33 H(x), x = 0(.0001)1(.001)5.6(various), 15D 1 — 36 
H (x), x = 4(.01)10, 8S 
Higher Integrals and Derivatives, etc., Involving e~** 
a Bell Aircraft Corp. 86 Circular normal integrals 210 
205 1 
P(R) = — exp — [(x — a)? + (y — b)?/2]dydx, 
2a z24y2<R?2 
- R = 0(.01)4.59, ~/a? + B = 0(.01)3 
Legendre Polynomials P,, (cos @) 
George C. Clark & 116 P, (cos 6), = 1(1)80, @ = 1°(1°)180°, 6D 277 
Stuart W. Churchill 
Associated Legendre Functions, P,,, (cos 0), P,.™ (cos 0), etc. 
79 S. L. Belousov 115 F,™ (cos @) or P,™ (cos 6) 276 
2n + 1 (n — m)! 
= ———— 2. (cos 6), 
79 / 2 (n+~m)! 
m = 0(1)36, m = m(1)56, @ = (2°.5)90°, 6D 
08 Bessel Functions: Unlimited Order 
Harold K. Crowder & 117 Jngg(x) and Vayg(x) = 1(1)50, m = 0(1) suchthat 277 
George C. Francis Vnsa(x) < 10! < Yuyg(x), 9D for nm < x, 7S for 
23 n2=x 
Expressions Involving J ,,(x) or Y,,(x). Spherical Bessel Functions 
30 and Riccati-Bessel Functions 
- _ 
Harold_K. Crowder & 117. jn(x) = Wz In44(x) and yn(x) = iz Vny3(x), 277 
x 20 
George C. Francis x = 1(1)50, 2 = 0(1) such that Vay3(x) < 10% < 
97 Y,439(x), 9D for n < x, 7Sforn = x 
Zeros and Turning Values of Bessel Functions 
F. B. Hildebrand 98 ay.” = 1(1)5, P = 0(1)5, 3D, where J,(a,,) = 0 218 
Modified Bessel Functions of the First Kind: Integral Order 
9 E. E. Allen 36 Io(x), —3.75 S x S 3.75, 8D (polynomial approxi- 38 
mation) 
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Author/s 


Number 





TABLES FROM ANALYSIS 


Review 


Description 


Auxiliary Functions for the Determination of Modified Bessel Functions 


E. E. Allen 


T. Pearcey 


J. Arthur Greenwood & 
David Durand 


M. Ferentz & C. Harrison 
T. H. Southard 


F. H. Hollander & 
C. B. Tompkins 


D. H. Lehmer 
D. H. Lehmer 


36 = In(x), —3.75 S x S 3.75, 8D; Io(x)xte*, 3.75 < 
x < o, 7D; (x)/x, —3.75 Ss x S 3.75, 8D; 
Ii(x)xte=, 3.75 S x < &, 7D; Ko(x) + 
log. (.5x)Io(x), O < x S 2, 7D; [Ki(x) —- 
log. (.5x) Ii (x) x, 0 < x S 2,7D; Ko(x)xte#, 2S 
x < », 7D; Ki(x)xte, 2 <5 x < ~,7D; (poly- 
nomial approximations) 


Fresnel Integrals 


87 C(V/2x/r) and S(+~V/2x/r), x = 0(.01)1, 7D, x = 
1(.01)50, 6D, & 


Integrals Involving Bessel Functions 


17 P(r,n) =r [ [Jo(x) Si (rx)dx, n = 6(1)24,7r = 
5(.5)12(1)2, 5D 

88 x7} [ Jo(y)dy, x = 0(.01)31, 4D 

‘Walesstontatinn Elliptic Functions 


57 PBs; 0, 1) — 1/22, g = 0(.1).8(.05) 1.55, 7D with 
Sn?. Also x*9 (x; 0, 1), 0 < x S 1.53, 7S (poly- 
nomial approximation ) 

58 x(x; 0, 1),0 <x S 1.53, 7S (polynomial approxi- 
mation) 


The Riemannian Zeta Functions 


52 10,000 zeros of ¢(s) 
108 25,000 zeros of ¢(s) 


Hypergeometric Functions, Confluent Hypergeometric Functions 


C. E. Fréberg 


35 Formulas for computation of Coulomb wave func- 
tions for untabulated arguments about 5S 


Debye Functions, Radiation Integrals, Fermi-Divac Functions, etc. 


D. K. C. MacDonald & 
Lois T. Towle 


2'dz 
9 #—1) (l-e4) 


r = 2,3, 4,6," = .1, .25, .5, 1, 1.2, 1.5, 2(1)6, 8, 
10, 13, 20, ~, 4 to 6S 





37. J, (x) = 


Miscellaneous Integrals Involving Circular and Exponential Functions 


S. Johnston 


R. B. Dingle, D. Arndt, & 
S. K. Roy 


z 
34  Sievert’s integral, | e-Avec®Jg A = 0(.5)10, 
0 


x = 0(1°)90°, about 5S 


119 Integrals won | (2 + x)~e-*de, 
0 


x = 0(.1)1(.2)2(.5)10(1)20;¢ = 6 = 1, 

p = —.5(.5)4;@ = 1,6 = 2, p = 0(.5)4; a = 2, 
b= 1,p = —5(.5)5;4 = b = 2, p = 0(.5)6.5; 
mostly 4S 


Page 


38 


210 


211 


110 


110 


107 


273 


37 


38 


36 


279 
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Review 
"ge Author/s Number Description Page 
6 {? 2 
M. Dank & S. W. Barber 120 The function of — =. —_ = . 281 
38 x? 0 ef —1 e-—i 
x = 0(.1)9.4 and 10(.5)16, 5D 
Miscellaneous Functions 
0 
O. Emersleben 56 The Epstein Zeta Function (p)Z H (>), ? = 109 
00 
1(1)4(2)8, 6D; °° (2), (p, g) = (0, 4), 
(34, 4), (4, 34), (0, 4), C6, 46), 26, 34), C6, 6), 
™ <p; 2/00 (1), 2 44,60; 2|¢| (0) 6D 
’ Prp ’ , ’ ’ ; ’ 
F. H. Hollander & 58 Solutions of y’ + (A + B cos 2? + Ccos4#+ 110 
C. B. Tompkins D cos 6t)¥Y = 0, for A = —.5(.1) —.2(.05) — 
28 -15(.03).15(.05).2(.1).5, B = 0(.2).8(.1) 1.6(.2)2; 
C=0,D = —.05(.05).10 and C = —.5(.1).5, 
D = 0, about 6D 
211 G. Belford, 59 Solutions of Hill equations 110 
L. Jackson Laslett, & 
J. N. Snyder 
110 Ernst Glowatzki 61 dm = (sin 6)#sn(mK/n; sin 6), 6D, A = a,‘Ila?,. 111 
with « = —1 for m odd and e = 0 for m even, 6S, 
—In A, 3D, 6 = 0(1°)90°, m = 1(1)n,m = 1(1)12 
= PN 1\* 
110 A. Rahman 125 Integrals / (3 In : + ) Me-“7dr, « = 0, 1,n = 285 
. oad 
0(1)10, @ = .1(.1)3, 10S 
107 Interpolation Coefficients and Formulas 
273 G. A. Chisnall 6 Modified j-th differences, 7 = 2(2)10 and multipliers 24 
for Chebyshev-Everett interpolation 
NBS AMS No. 35 54 Lagrangian interpolation coefficients, 3 to 6 points, 108 
37 @ = 0(1’)60’, 8D 
H. E. Salzer 55 Complex osculatory interpolation in a square, 2to4 109 
points, ~, g = 0(.1)1, exact 
K. A. Karpov 85 Lagrangian coefficients A;(t), i = —1(1)2,¢ = 209 
—1(.001)2 andi = —2(1)2,# = —2(.001)2, 6D 
38 Frank E. Harris 118 Ra(k) = nI[1 — e*2te A*/k!], h = .01(.01).05, 279 
.1(.1).5, 1, values of m and precision as needed for 
interpolating Ei(x)4 < |x| < 50, 18S 
Extrapolations 
Herbert E. Salzer 9 Am(n)-m = 4(1)10, 25 
36 nm = 11(1)50(5)100(10)200(50)500(100) 1000, 
exact, where S,, the sum of » terms of a sequence, 
is estimated by 2, An (n)S, 
R. B. Horgan 10 Am(m), m = 4(1)10, 25 
79 n = 11(1)50(5)100(10) 200(50) 500(100) 1000, 
15D 
Coefficients for Integration and Summation Using Ordinates 
P. Davis & P. Rabinowitz 84  Abscissas and weights for Gaussian quadrature with 209 
2, 4, 8, 16, 20, 24, 32, 40, and 48 points, 20D 














332 TABLES FROM STATISTICS 


Review 
Author/s Number Description 
H. Fishman 136 Best values for b; and x; in 


1 
/ x"g(x)dx = DR 1bjg(x;), 
0 


nm = 1(1)5,” = 1(1)8, 12D 


Coefficients for Integration and Summation using Differences 


2 \2 
K. S. Kunz 8 2S, — S.-,] (") and (*").2 = 0(1)6,” = 1(1)p, 
exact; where S, = 2," 7-1(2n + 1)!/(n!), n = 
1(1)6, exact 


Numerical Integration, Multiple Integrals 


P. C. Hammer, 7 Polynomials gn" orthogonal with respect to x” on 
O. J. Marlowe, & (D, 1), m = 1(1)5, 2 = 1, 2 and m = 1(1)4, 
A. H. Stroud n = 3; zeros of gn” and weights and arguments 


for numerical quadrature over a triangle, 18S 


Bibliographic Material Relating to Tables from Analysis 


A. V. Lebedev & 49 Guide to tables 
R. M. Fedorova 


TABLES FROM STATISTICS 
Normal Distribution. Area, log Area and Ordinates as Function of x/./2 


2 - . 2 
NBS AMS No. 41 33 H(x) = = e-*"da and H'(x) = wa es, 
tT Jo T 


x = 0(.0001)1(.001)5.6(various), 15D and 1 — 
H(x), x = 4(.01)10, 8S 


Normal Distribution. Measures of Location and Dispersion 
other than the Standard Deviation 


G. E. Noether 14 For subsample of size N, factors to estimate o from 
range, 3S and to give 90% and 98% confidence 
limits, 2D, N = 2(1)100 


D. P. Banerjee 15 80, 90, 95, and 99% points 
D. A. S. Fraser & 16 Tolerance region factors, k variates k = 1(1)4, 
Irwin Guttman n — k = 1(1)30, 40, 60, 120, ~, 8 = .75, .9, .95, 


-975, .99, .995, 4S 


Normal Bivariate Distribution 


Peter Ihm 18 A 3D table for 95% and 99% confidence that a sam- 
ple from bivariate normal distribution indicates 
1 — 02 =p =0 


1 
Bell Aircraft Corp. 86 Circular normal integrals P(R) = an 
T 


I. ee C(x — a)? + (y — b)*/2dydx, 
z*+y 


R = 0(.01)4.59, 
Va + @ = 0(.01)3, 5D 
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TABLES FROM STATISTICS 
Review 
Author/s Number Description 
Normal Distribution. Miscellanea—Quality Control 
R. L. Storrer & 30 k, 2D and F, 3D, for quality levels of .065, .1, .15, 
W. R. Davison .25, .4, .65, 1, 1.5, 2.5, 4, 6.5% and sample sizes 


5(5)40, 50, 70, 130, 200 for inspection by average 
range and sample standard deviation, one and two 
sided specifications U and L suchas X + k&R S U, 
or Max allowable average range = F(U — L) 


G. J. Lieberman & 31 ‘Multi-level continuous sampling plans” 
Herbert Solomon 
H. F. Dodge 32 “Chain sampling inspection plan” 


Auxiliary Tables for Estimates or Tests Involving f Variates 
Sigeiti Moriguti 25 99.9% confidence limits for a variance component 
connected with x? distributions with degrees of 
freedom f = 6(2)12, 15, 20, 30, 60 and f; = 

1(1)6(2)12, 15, 20, 30, 60, ©, 45 or 3D 


Cumulative Binomial Series, Tail Area of I_(a, b) 


Harvard University 11 Cumulative binomial distributions E(n, r, p) = 
n 


2:-,(7) ea — p)"",r = 0(1)n, p = .01(.01).5, 


Me, Ma, , ¥, He, Ke, 4, 34, Ke, Ke, = 
1(1)50(2) 100(10) 200 (20) 500(50) 1000, 5D 


Percentage Points of I. (a, b). Confidence Belts for Binomial Distribution 


Motosabura Mauyama 20 5% and 1% control limits for limits for r-th order 
observations, N = 3(1)6(2)10,r = 1(1)N 


The Beta Function, Moments of Binomial Distribution, Moments of I.(a, b) 
J. L. Hodges, Jr. 22 «4ttC,, A = .5(1)19.5, ¢ = 1(1)18, 7S 
Shoji Ura 26 vw =([(fit 1)/fede, fi = 1(1)10, 12, 15, 20, 24, 30, 
60, 120, ~, fe = 2(2)20, 24, 30, 40, 60, 120, @~, 
2D, where g is Tang-Lehmer power function, with 
significance level .05 error probability. 


Miscellaneous Functions Based on the Binomial Distribution 
J. L. Hodges, Jr. 22 =AttC, A = .5(1)19.5, ¢ = 1(1)18, 7S 


Tables of F-Distributions and Multivariate Variance Tests 
D. P. Banerjee 15 80, 90, 95, and 99% points of sample standard devia- 
tions, normal bivariate universe, N = 3(1)30, 
p = 0(.1).9, where p is the universe correlation 
coefficient 


Miscellaneous Continuous Distributions 
L. H. Miller 23 such that a = Prob (D, 2 e), a = .005, .01, .025, 


.05, .1, N = 1(1)100, where S, is observed, F true 
distributions, D, = max {S,(x) — F(x)} 


Geometric Probability Distributions 
J. Arthur Greenwood & 17 Distributions of R, = [(E>—, cos ty)? + (Z>_1 sin 
David Durand ty)2}}, where & is uniform variate (0, 27), 2 = 
6(1)24, Ra S .5(.5) 12(1), 5D. Also 95% and 99% 
points, 3D, etc. 
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334 TABLES FROM STATISTICS 





Review 
Author/s Number Description 
Distribution of Serial Correlations and Related Statistics 
E. J. Hannon 13. The following functions for testing correlation be- 
: . (1 + Pipe) (1 — p?) 
tween time series. and 
(1 — pipe) (1 — pi*) 


(1 — pips?) (1 + pips)? 
2(1 — pit) (1 — pot) 
and 
(1 + Pips — 5)*[1 + pid) — 5(p:* + Pips) ] 
Cl — pid) — b(pip2 — pr*) 
2(1 — pit) (1 — ps?) (1 — b) [1 + pe? — B(1 — D2?) ] 
b= —.6(.2).6, pi, po = .6, .8,and pi = pe = A, 2D 
J. F. Scott & V. J. Small 28 Tables for regression coefficient between trend re- 
duced time series approximately Markoff 





, fi, P2 = —.8(.2).8, 2D, 





Multiple Correlation Coefficients 


F. M. Lord 19 Riss = {r2, + 72, — 2niFi9723)/1 — 72,)}4 (nomo- 
graph) 
Non-parametric Two Sample and Analysis of Variance Randomization Tests 
E. S. Page 29 5 and 1% points for m = max{S,— min S;}, 
0S tSn 0St<r 


S, = Z{ sgn (x; — 0), x; symmetric about @, 27 
values of m 21 S m S 118. Power comparisons with 
sign test 


Runs 


Michio Takashima 27. ‘tg and (tg), a = .01, .05, m, m = 1(1)25, exact, 
where ¢, and (¢,) are smallest integers giving prob- 
ability Sa of run of length ¢ classes of either or 
selected one of two classes of objects 


Tables Relating to Efficiency and Goodness of Fit 


M. R. Sampford 21 Efficiencies of moment method for fitting truncated 
negative binomial distribution, Rk = .5, 1(1)5, 
p = .5, 1, 2, 5, 3S 


Miscellaneous Tables Relating to Fitting of Frequency 
Distributions, Moments of Moments and Cumulants, 
Means from Non-normal Parents, etc. 
C. W. Topp & F. C. Leone 24 a;*andé,b = 1,7 = .01, .02, .05, .08, .09, .1(.05).95, 
a = .05(.05)1, where a; is 3rd moment, 
8 = (2a, — 3az? — 6) (a, + 3)! of F(x) 


a x 
= 5 (2b% — x*)* + (1 — a) 5? *¢ 00,1] 


Random Deviates from Mechanical Devices—Quasi Random Numbers 


NBS AMS No. 12 42 Methods of generation and test 
H. A. Meyer 43 Generation and test 
Random Normal Deviates 
E. C. Fieller, T. Lewis, & 12 %ay, a = 0(.1).9,¢ = 1(1)3,000 where x,; ¢ N(0, 1) 
E. S. Pearson and a ~ 0 is correlation coefficient between x»; 


and xq;, 2D 
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TABLES 


FROM PHYSICAL SCIENCES 


Description 


Bibliographic Material Relating to Statistical Tables 


Guide to Mathematical Tables 


TABLES FROM PHYSICAL SCIENCES 


Review 
Author /s Number 
A. V. Lebedev & 49 
R. M. Fedorova 
G. Belford, 59 
L. Jackson Laslett, & 
J. N. Snyder 
Sir Harold Jeffreys & 97 
Bertha Swirles 
E. F. Beckenbach, Editor 124 
A. Rahman 125 
F. H. Hollander & 58 
C. B. Tompkins 
G. Belford, 59 
L. Jackson Laslett, & 
J. N. Snyder 
Marcel van Laethem 72 
Speer Products Co. 91 
Homer S. Powley 92 
E. F. Beckenbach, Editor 124 


R. E. D. Bishop & 
D. C. Johnson 
Zygmunt Dowgird 


Leonard Pode and 
Leonard Pode & 
Louis Rosenthal 

E. L. Harris & 

G. N. Patterson 


Einar Anderson, 
Torben Krarup, & 
Bjarner Svejgaard 

Sverre Petterssen 


General References 
Solutions of Hill equations 


Methods of Mathematical Physics 


Modern Mathematics for the Engineer 
Tables of Integrals A,(a) = / A*e- “dd and F, (a) = 


1 
| Qo(A)A* e~ “dr 
1 
Particle Mechanics 
Solutions of y’ + (A + B cos 2¢ + C cos4#+D 
cos 64)¥Y = 0, for A = —.5(.1) —.2(.05) — 
-15(.03).15(.05).2(.1).5, B = 0(.2).8(.1) 1.6(.2)2; 
C =0,D = —.05(.05).10 and C = —.5(.1).5, 


D = 0, about 6D 
Solutions of Hill equations 


Ion trajectories in cylindrical condensor 

Ballistic slide rule 

Ballistic functions A V?/700*, log V/u, N for Z = 
20000(1000)45000 and 50000, 2D, 4D, and 2D 
respectively 

Exterior ballistics (essay) 


Mechanics of Rigid and Deformable Bodies 


89 


96 


123 


130 


90 


129 


Vibration Analysis Tables 


Cracovians 
Fluid Mechanics 
Tables concerning equilibrium configuration of a 


flexible cable in a steady stream 


Boltzmann’s H for gas with a shock wave, Mach 
1.5, 2, 2.3238, 3, 3.3764, 4 


Geophysics 
Geodetic Tables, International Ellipsoid 


Weather Analysis and Forecasting 


335 


Page 


104 


110 


217 


285 


110 


110 


120 
214 
214 


283 


212 


216 


282 


288 


212 


287 











336 TABLES FROM PHYSICAL SCIENCES 
Review 
Author/s Number Description 
Thermodynamics, Statistical Mechanics, etc. 
6 ff? #dé 2x 





C 
M. Dank & S. W. Barber 120 Specific heat function — = 


= eJoe@—-1 @—1' 


x = 0(.1)9.4 and 10(.5)16, 5D 
Exchanges, Reaction Kinetics, Reaction Rates, etc. 


J. Crank 71 The Mathematics of Diffusion 
Electricity and Magnetism; Electric Waves; Electric Filters 
Ernst Glowatzki 61 Cauer parameters @ = 0(1°)90°, m = 1(1)n,” = 
1(1)12; am to 6D, A to 6S, —In A to 3D 
T. Sasaki 62 Coastal bending of radio waves 
Jg@rgen Rybner 80 sinh (A + 2B) and cosh (A + 7B),0 S A S 4; 


tanh (A + ¢B),0 S A S 3 (alignment charts); 
several auxiliary functions for filter design 


R. B. Dingle, D. Arndt, & 119 Integrals (p!)-} [ (ee + x)—e-de, 
0 


S. K. Roy x = 0(.1)1(.2)2(.5)10(1)200;@ =b=1,p = 
—.5(.5)4;¢@ = 1,6 = 2,p = 0(.5)4;¢ = 2,56 = 1, 
pb = —.5(.5)5;@ = b = 2, p = 0(.5)6.5, mostly 


4S for application to semi-conductors 


Crystal Structure, X-ray Diffraction, Electron Microscopy 


R. A. Sparks, R. J. Prosen, 38 Summary description of automatic computer codes 
F. H. Kruse, & for crystal structure calculations 
K. N. Trueblood 

O. Emersleben 65 Madelung constants ‘ 

O. Emersleben 66 Madelung constants 


Atomic and Molecular Physics and Chemistry, Quantum Mechanics 








. "d nie 
D. K. C. MacDonald & 37. Ir(x) = [ tr = 2,73, 4, 6, 
Lois T. Towle aMor Do ) x 
x = .1, .25, .5, 1, 1.2, 1.5, 2(1)6, 8, 10, 13,720, @, 
4 to 6S 
Frank E. Harris 118 Ei(x),4 <|x| < 50, 18S 
~ A+1\" 4, 
A. Rahman 125 Integrals 4 In > " Aj" “MD, e = 0, 1,8 = 
‘ ms 
0(1)10, 2 = .1(.1)3, 10S, for calculating wave func- 
tion of hydrogen molecules 
A. Rahman 128 Integrals from 2s and 2p atomic functions, 5D 
Nuclear Physics and Scattering 
(Including Photon Scattering and Compton Effect) 
C. E, Fréberg 35 Formulas for computation of Coulomb wave func- 
tions for untabulated arguments, about 5S 
H. A. Meyer 43 Monte Carlo methods 
J. L. Wolfson & 64 Electron kinetic energy for momentum 
H. S. Gellman B, = 100(1)2000(10) 10000(100) 20000 Gauss-cm, 


5S, PP 


Page 


119 


111 
207 


279 


39 


113 
1i4 


38 


279 


285 


286 


37 


44 
112 
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rT 


AUTOMATIC COMPUTATION 337 
Review 
Page Author/s Number Description Page 
K. M. Case, F.de Hoffman 126 Theory of neutron diffusion with some tables 285 
& G. Placzek 
281 M. Frentz & 127  F,(L, L’, 7’, 7), 8D, for 13,602 values of arguments 286 
N. Rosenzweig 
WORKS CONCERNING SOCIAL SCIENCES 
= Operations Research, Simulation 
E. F. Beckenbach, Editor 124 Applied mathematics in operations research, Monte 283 
lll Carlo methods (essays) 
Econometrics 
111 
207 J. Halcomb Laning, Jr. & 63 Random Processes in Automatic Control 112 
Richard H. Battin 
S. Vajda 93 The Theory of Games and Linear Programming 215 
Proc., Second Annual 94 Linear programming 215 
279 Computer App. Symp. 
E. F. Beckenbach, Editor 124 Dynamic programming (essay) 283 
Management Science (for Programming see Econometrics) 
Proc., Second Annual 94 Symposium proceedings 215 
Computer App. Symp. 
National Phys. Lab. 105 Wage Accounting by Electronic Computer 225 
39 Learning Theory, Psychometrics, etc. 
W. Sluckin 141 Mind and Machines 305 
113 Competition, Attention, Warfare, etc. 
1i4 S. Vajda 93 The Theory of Games and Linear Programming 215 
E. F. Beckenbach, Editor 124 Theory of games (essay) 283 
Translation, Pattern Recognition 
38 Proc., Inst. of Elec. 73 Conference report, several papers on mechanical 121 
Engineers translation and on character recognition 
” Biological Applications 
9 
S."Johnston 34  Sievert’s Integral, | e~4 20949, A = 0(.5)10, x = 36 
285 ° 0 
0(1°)90°, about 5S 
286 AUTOMATIC COMPUTATION 
General Works on Automatic Computation 
M. V. Wilkes 41 Automatic Digital Computers 42 
Clarence L. Johnson 45 Analog Computer Techniques 48 
37 Proc., Inst. Elec. 73 Report of convention on digital computer techniques 121 
Engineers im | 
44 S. A. Lebedev 74 Electronic Computing Machines,and Information 124 
112 Processes, pamphlet 
George R.’Stibitz & 100 Mathematics and Computers 221 
Jules A. Larrivee 











Author/s 
W. J. Eckert & 
Rebecca Jones 
E. F. Beckenbach, Editor 


Heinz Rutishauser 
F. J. Weyl 


Milton H. Aronson 

L Lf 

Elektronische Rechnen- 
maschinen und Informa- 
tionsverabeitung 


Number 





AUTOMATIC COMPUTATION 


Review 
Description 
101 Faster, Faster 


124 High speed computing devices and their applications, 
essay in book on mathematics in engineering 

132 Applications of the quotient-difference-algorithm 

137 Report on a Survey of Training and Research in 
Applied Mathematics in the United States 

138 The Computer Handbook, general articles on digital 
and analogue computing 

142 Report of symposium on electronic computers 


Radix Conversion—Decimal to Radix 2* and Radix 2* to Decimal 


Carl-Erik Fréberg 


A. van Wijngaarden 


Proc., Inst. Elec. 
Engineers 


Proc., Inst. Elec. 
Engineers 

Proc., Symposium of 
Franklin Institute 

D. D. McCracken 

Elektronische Rechnen- 
maschinen und Informa- 
tionsverabeitung 


Hans J. Maehly 


E. E. Allen 


T. H. Southard 


F. H. Hollander & 
C. B. Tompkins 


82 Hexadecimal conversions: integers 1 (1) 1024(16) 4096 
and 10*(10*)10**!, k = 2(1)12; fractions x = 
10-*(10-*) 10-**2, k = 2(2)16; m-10*, nm = 1(2)9, 
k = 12(1)12; 10* and 10-*, k = 13(1)25; some 
constants, all in hexadecimal form. 

x = 16-*(16-*) 16-41, 
k = 1(1)10 in decimal form 

83 Decimal octal conversion x = 0(8)99992; 2%,” = 
1(1)50 exact; 2-*, m = 1(1)50, 20D. Octal 10*, 
nm = 1(1)18 exact and 10, m = 1(1)18 to twenty 
octal digits 


Analogue-Digital Conversion 


73 Papers in report of convention on digital computer 
techniques f 


Logic of Coding 
73 General coding strategy 


139 Proceedings of a symposium on automatic coding 


140 Digital Computer Programming, a text book 
142 Parts of report of symposium on electronic computers 


Rational Approximations 


4 loge x, 4 S x S 1, 10D; 2, —4 S x S 3, 10D; 
arc tanx, -1 Sx*x31,9D 

36 Bessel functions: Jo(x), —3.75 S x S 3.75, 8D; 
Ih(x)/x, -3.75 S< x S 3.75, 8D; In(x)xte=, 
3.75 Sx < @, 7D; Ko(x) +log. (.5x)Io(x), 0 < 
x S 2, 7D; [Ki(x) —log, (.5x) li(x) ]}x, 0<x S 2, 

7D; Ko(x)xtet, 2 S x < @, 7D: Ki(x)xtet,2 Sx < 

©, 7D, polynomial 

57. Weierstrass elliptic function, ep (z,0,1),0<28 
1.53 7S 

58 x*@B(x, 0, 1),0 <x S 1.53, 7S 
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296 


306 


208 


121 


121 


296 


298 
306 


24 


38 
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208 


24 


38 


110 


110 





Author/s 
V. Belevitch & F: Storrer 


R. K. Richards 

J. Halcomb Laning, Jr. & 
Richard H. Battin 

Proc., Inst. Elec. 
Engineers 

W. Sluckin 

Elektronische Rechnen- 
maschinen und Informa- 
tionsverabeitung 


AUTOMATIC COMPUTATION 


Review 
Number 


78 


Description 
x71, «74, sin x, arc tan x, 10%, 10° — 1, logis x, 15S, 
polynomial 


General Aspects of Machine Design 


48 
63 


73 


141 
142 


Arithmetic Operations in Digital Computers 
Random Processes in Automatic Control 


Papers in report of convention on digital computing 
techniques 

Mind and Machines 

Papers in proceedings of symposium on digital com- 
puters 


Machine Descriptions—Specific Electronic Machines 


Proc., Inst. Elec. 
Engineers 


S. A. Lebedev 

W. J. Eckert & 
Rebecca Jones 

Elektronische Rechnen- 
maschinen und Informa- 
tionsverabeitung 


73 


101 


142 


DEUCE, MERCURY, PEGASUS, ELLIOTT 400 
series, HEC, IBM 650, 704, 705, MANCHESTER 
MARK II, NICHOLAS, EDSAC II, IMP, ACE, 
BESM, in report of convention on digital computer 
techniques 


Vienna Logical Function Comp., IRSIA-FNRS, 
SAPO, Dresden D1, PERM, DERA, WETZLER 
Z5, Géttingen G1, G2, G3, PTERA, ARRA, 
ARMAC, BESK, ENMETH, IBM 705, BESM, 
URAL, in report of symposium on electronic 
computers 


Machine Programs and Codes for Specific Problems 


R. A. Sparks, R. J. Prosen, 


F. H. Kruse, & 
K. N. Trueblood 

T. R. Norton & A. Opler 

Num. Anal., Proc., Sixth 
Symp. in App. Math. 
AMS 

Proc., Inst. Elec. 
Engineers 

Proc., Second Ann. 
Computer App. Symp. 

D. H. Lehmer 

L. Pease 

Sverre Petterssen 


George E. Forsythe 


W. G. Brombacher, 
Julian F. Smith, & 
Lyman M. van der Pyl 


38 


39 


44 


73 


94 


102 


103 
129 


104 


Crystallographic calculations ° 


Mechanized search for organic compounds 
Generation of permutations 


Several problems considered in convention on digital 
computer techniques 
Several papers in computer applications symposium 


“Sorting cards with respect to a modulus” 
Control board for IBM 602A calculating punch 
Numerical weather forecasting 


Computing Bibliography 


“Selected references on use of high-speed computers 
for scientific computation” 
Guide to Instrumentation Literature 
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Page 
205 


50 
112 


121 


305 


121 


124 
222 


39 


40 
46 
121 
215 


224 
225 
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225 








340 ERRATA, NOTES, AND CORRIGENDA 


TABLE ERRATA 
No. Author Title Page 
252 A. A. ABRAMOV Tablitsy \n T[z] v Kompleknoi Oblasti................. 125 
253 AKADEMIIA Nauk SSSR Desiatiéna Tablitsy logarifmov kompleksnykh chisel i 
perekhoda ot dekartovykh koordinat k poliarnym. Tablitsy 


PIR saetttgarccrsars<+s+csscscbecss ss tees 125-126 
254 AKADEMiA Nauk SSSR Tablitsy Integralov Frenalia..............00.00cceeeee 126 
255 ISABELLE ARSHAM Chebyshev Coefficients for Chebyshev Polynomials of Orders 
12 and 24 under the General Linear Transformation (U) 
126-127 
256 GEORGE WELLINGTON Smithsonian Logarithmic Tables to Base e and Base 10.... 226 
SPENCELEY, RHEBA 
Murray SPENCELEY, & 
EvuGENE RHODES EPPERSON 
257 C. A. Coutson & “Some new values for the exponential integral”’........ 308 
W. E. Duncanson 
258 CaRL-Er1Ik FROBERG Hexadecimal Conversion Tables............02220eeeeee 309 
Also see the following reviews: 
27 Micuto TAKASHIMA “Tables for testing randomness by means of lengths of 
SN Eee eET Rhee eh a CS dN sc cacciecccsses 33 
37 D. K. C. MacDonaLp & “Integrals of interest in metallic conductivity”........ 38-39 
Lois T. TOWLE 
77 H. NAGLER Table of Square Roots of Integers. ...........0000cceeee 205 
80 J¢RGEN RYBNER Nomogrammer over komplekse hyperbolske funktioner.... 207 
82 Car_-Erik FROBERG Hexadecimal Conversion Tables... .......20.ccccecsees 208 
87 T. PEARCEY Table of the Fresnel Integral to Six Decimal Places... 210-211 
107 D. N. LEHMER List of Prime Numbers from 1 to 10006721............. 272 
ao 
119 R. B. Dinezz, D. Arnpt, & “The integrals A,(x) = on | &(e + x)—e~de and 
0 
‘co 
S. K. Roy B,(x) = w= | I os ek evesecccceasis 279 
0 
% NOTES 
~seesete Sealine per @ promreaso delle Gciense. oo... oo... 65 os on svn ccc ccnc cece cwecenes 51-52 
ek cs ke bs ns 6 6 he 06. ss 49.505/0> See eee ee 53-54 
i ahi ch Es atin Black sie.e ann nee.s Tie opse-e aeecis ee 6-48 127-128 
ni cred add nines die. omaacee baanuees meni 129 
Wayne State University—Conference on Matrix Computations....................000005 129 
Handbock of Mathematical Tables—National Bureau of Standards.................. 226-227 
International Congress of Mathematicians—August 14-21, 1958........... lies sna pene 227-228 
ne I NS, PS, SEU TEA tess a ccc ccc sect onesie geceaeys 228 
ee lia ice ait 0.6 s ldklia 0 kes bch scheme eh oh esitnens 309 
IE EE TT A 309 
Policy Committee for Mathematics Technical Advisory Committee to the National Bureau of 
I SN, MN an 5 ds bh punto cindy < sis 64.00 66-0 che cdakarts.s 310 
Ree a SE RA a as alcidin «0 aaa Seip ya Peete ek Ral ama <P Es Mis gh edsloune o> Sgt dp kgs b's 4420 313 
CORRIGENDA 
TO PAPERS: 
Authors Title Page 
Mark LorKin (EEE E ET OTE, ae ere 130 


D. R. Morrison “A method for computing certain inverse functions’’....... 314 
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CORRIGENDA 341 


Title Page 
ee Me FEES GRE BUD Bio 6 kk nis de ds kc csscvscsnsees 129 
Proc. of the Sixth Symp. in App. Math., AMS............. 228 
Acht- und neunstellige Tabellen zu den elliptischen Funktionen 130 
“Bounds for the ratio of range to standard deivation”..... . 54 
“A family of J-shaped frequency functions”.............. 314 
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Arithmetical Tables. Mathematical Constants 
Powers 


Logarithms 


o 7.97 > 


Circular Functions 


Hyperbolic and Exponential Functions 
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Theory of Numbers 


Higher Algebra 


m © 


Numerical Solution of Equations 
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Finite Differences. Interpolation 
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Statistics 

Higher Mathematical Functions 
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Ze Pr Ro 


Interest and Investment 
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Actuarial Science 

Engineering 
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Navigation 

Aerodynamics, Hydrodynamics, Ballistics 
Economics and Social Sciences 


Numerical Analysis and Applied Mathematics 
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